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Preface 


This book is intended to be a realistic introduction to the theory 
of solid state physics. It is meant to fill the gap between the more 
empirical books such as Kittel’s Introduction to Solid State Physics 
or Dekker’s Solid State Physics, and the more advanced theoreti- 
cal texts such as Kittel’s Quantum Theory of Solids. No claim is 
made that the subject is covered comprehensively. This is 
clearly not possible within a book of this size. However, it is 
hoped that the included material will give a good basic back- 
ground to the subject. 

This text should be useful to research workers both in industry 
and university, and to final honours students. It should be 
especially appropriate to the several new M.Sc. courses in Solid 
State Physics. The reader is assumed to have taken a first course 
in Quantum Mechanics. 

Problems are included at the end of each chapter and a few 
of these have been used to extend the theory slightly. M.K.S. 
(Sommerfeld) units are used throughout and the charge on the 
electron is taken as —e = —1-602 X 10°C, 

H. CG. 
Tynemouth, March 1968 
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The One-electron Approximation 


§ 1.1 The quantum mechanics 


When it is noted that the interatomic distance in solids is of the 
order of a few angstrom units, it is not at all obvious how it is 
possible to separate the motion of the electrons in any sense. 
This monograph is primarily concerned with the results ob- 
tainable from the one-electron approximation and it is the 
purpose of the first chapter to indicate the relevant theoretical 
background. Initially it is necessary to explain how the motion 
of the electrons can be treated apart from the lattice of heavy 
ions and then to deal with the separation of the electrons them- 
selves. 

The results of the Schrédinger wave mechanics will be used 
freely throughout this book and perhaps it is as well to recount 
briefly the main results of quantum mechanics. For further in- 
formation the reader is referred to any of the standard texts such 
as Schiff or Mandl. 

In quantum theory, the state of any physical system is de- 
scribed as completely as possible by a wave function % which in 
general is complex and depends on the position co-ordinates of 
all the NV particles and on the time. Strictly 4 depends upon the 
spin co-ordinates as well. This will be emphasized where appro- 
priate. 

The real quantity %* dr is the probability that the co- 
ordinates 9; of the particles lie within the volume element 
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dr = dgidge ae dgsn 
at the time ¢. 

To be physically meaningful, the wave functions must be 
single-valued and in general they are continuous and quad- 
ratically integrable. The latter property makes it possible to 
normalize the wave function to unity, 


i.e. | ptibdr = 1 


The physical observables like energy and angular momentum 
which are present in classical theory in the form of dynamical 
variables are represented in quantum mechanics by operators 
which are both linear and Hermitian. A linear operator A has 
the property 


A(ari + dpe) = aAya + bAde (1.1) 


where #1, #2 are two state functions and a, 6 are constants. ‘The 
Hermitian condition demands that 


[tr*Adade = ( | o* Addr)” (1.2) 


There are some general rules to act as a guide in finding the 
form for these operators. In the Schrédinger representation they 
are: 


(i) the operator representing one of the position co- 
ordinates q; or the time ¢ is simply q; or f. 
(ii) the operator representing one of the components of 


e « h 0 
momentum f; conjugate to q; 1s 7 Oa, 
i 


The quantum operator representing a dynamical variable 
which is expressed in terms of the position co-ordinates, mo- 
mentum and time is obtained by substitution using the above 
rules. 

As an example consider the angular momentum of a point 
mass moving in some orbit. Classically the angular momentum 


2 


THE ONE-ELECTRON APPROXIMATION 


L of such a particle about an axis through the origin is given by 
the vector product 

L=rA\p 
where r is the radius vector and p the momentum of the particle. 


For a particle moving in the x—y plane the angular momentum 
vector lies along the z-axis and 


Lz = Xpy — Iba 


The quantum mechanical operator which represents the 
angular momentum is 


h/ a 0 
f= 5 Jan) 
If the state function ¢ is an eigenfunction of the operator A 
corresponding to a classical variable, 


i.e. Aydt = ays (1.3) 
where a is some constant, then in this state the variable has pre- 
cisely the value a.a is an eigenvalue, It can be shown that the 
eigenvalues of a Hermitian operator are real and that if ¥; and 
ys; are eigenfunctions belonging to different eigenvalues then 
they are orthogonal, 


i.e. | WtWdr =0 i ¥j 
Two operators A and B are said to commute if 
(AB) = (BA)y 


when applied to an arbitrary wave function %. Quantum 
mechanical operators do not always possess this property, The 
operators representing the position co-ordinate x and the con- 
jugate momentum fp, do not commute, 
hl d a h 
1.€. =| q-X — X57 => 

( Ox =) 2 id 
When two operators commute it can be shown that there exists 
a set of functions which are simultaneously eigenfunctions for 
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both operators. This result is important in the classification of 
states and is referred to again in Chapter 3. 

In classical mechanics the Hamiltonian H defines the total 
energy for a conservative system. The time variation of the 
state is determined through Hamilton’s equations. The most 
important quantum mechanical operator is the analogous 
Hamiltonian # obtained by substitution. The dependence of 
the state function % on the time is given by 


hop 
Hy = —=7e (1.4) 
This is the time-dependent Schrédinger equation. 

If y% is an eigenfunction of the Hamiltonian operator, the 
eigenvalue gives the energy E of the state and the wave function 
merely changes its phase in time, 

1.e; Hb = Eb (1.5) 
and Ye(E) = po(E)e teen (1.6) 
Equation (1.5) is the time-independent Schrédinger equation. 


In an eigenstate of the energy, the average value of any 
observable A, defined by 


| Yr Abd 


is a constant if A does not explicitly involve the time. The wave 
function is then said to define a stationary state. 

For those readers completely new to solid state theory it is 
possible to omit the remainder of this chapter on a first reading. 


§ 1.2 The Born-Oppenheimer approximation 


The classical Hamiltonian for a collection of N particles of 
masses m; and charges e& respectively, is given by 


2 1X4 B e1€} 
H= ‘cegeee I. 
Yom * 22,,2.4ne0 [ni — Fy (7) 
G#j) 
where r; is the position vector of the 7th particle. The first sum. 
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is taken over the 3N components of momenta and represents the 
kinetic energy of the system. The second term represents the 
electrostatic potential energy of the system and the sum is over 
all the possible pairs of particles. The corresponding quantum 
mechanical Hamiltonian is (omitting spin-dependent terms), 


—h? I eit; 
= 3am)" tak Fearon 

where V;2 operates on the spatial co-ordinates of the ith 
particle. The eigenfunctions of this operator must involve the 
3N position co-ordinates of the system (4N if spin is included). 

In a crystal there are two types of particles present, the 
nuclei and the electrons. For our purpose the structure of the 
nucleus may be ignored. 

The Schrédinger equation for a stationary state of a crystal 
may be written 


H (Ri, Re,...jr1,4e,...) FY (Ri, Re. . .571, 972, - -) 
= FE, ¥ (Ri, Re,...; 71, re,.-.) (1.9) 


where the wave function ‘¥ and the Hamiltonian are functions 
of all the nucleus and electron co-ordinates R; and r; respec- 
tively. The number of variables is reduced if the atomic 
electrons are separated into core and valence electrons, al- 
though this division is somewhat arbitrary. The core electrons 
have atomic wave functions which do not overlap to any appre- 
ciable extent at the observed interatomic spacing in the crystal 
and consequently may be considered to be unaffected by com- 
bining the atoms to form a crystal. Conversely, the wave 
functions of the valence electrons overlap considerably, and are 
modified by the process of crystal formation. 

The Hamiltonian of (1.9) then contains all the interactions of 
an assembly of ion cores and valence electrons and may be 
expressed as 


H( Ri, Ro. . .3 iy ¥25 +) 
= H# ion( Ri, Re, ee | +- # erectron(F1, Y2,.-. “) 
+ # erectron( Ri, Ro, ee e3 H1; F2;,-.- z) (1.10) 
ion 
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2 
where H ion = => aa Vr2z + ae LW (Ri — R;) 
t tH) 


H erectron = — ai Nae oP => > Z 
7 2m 29745 4rreo [ri — Fy 
and H# electron = >> v (ri — Rj) 
-ion tj 

The ionic Hamiltonian # ion describes the motion of the set of 
like ions, each composed of a nucleus and a shell of core elec- 
trons, which interact through a potential W(R; — Rj). #etectron 
corresponds to the valence electrons which interact through 
a Coulomb potential. The electron-ion Hamiltonian takes 
account of the interaction between the ionic cores and valence 
electrons which are assumed to act through a potential 
o(rs — R;). 

The total crystal potential may be written 


Ve = Vion + V(R1, R2, oo 05 1, F2,.- ) 
with 
Vion oe 4> >W(R: _ R;) 


t#j 
e 
(EEF gral + RP Be 


where Vion is the potential of interaction between the ion cores, 
and V is the potential of interaction of the electrons with the ion 
cores and with each other. 

Because the electronic masses are much less than the ionic 
masses, the electrons move much more rapidly than the ions. 
In the ‘Born—Oppenheimer’ approximation (1927), the 
Schrédinger equation is solved for the N electrons in the field 
of all the ions in some fixed configuration. The electronic wave 
function %(Ri, Re, ...; ri, re, ...) depends on the variable 
positions of the electrons r; and also on the nuclei-co-ordinates 
R; regarded as parameters. The full crystal wave function is 
written 


W(R1, Re, ...5 71, r2,...) = O(Ri, Re, ...). 
W(Ri, Re, ...3 "1, ¥2,.-.) (1.11) 
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where ® is the lattice wave function describing the ionic motion 
and is a function of the co-ordinates R; only. 
7 satisfies the equation, 


+ V(Ri, Ro,...3;m,4e2,.. )| (Ri, Re, ...5 71, r2,..-) 
= E(Ri, Re, ...) (Ri, Re, ...5 7, %,...) (1.12) 


where m is the electron mass and V,,2 operates on the electron 
co-ordinates. The electron energy eigenvalue E(R1, Re, ...) 
depends upon the particular ionic configuration considered. 
This electron energy then appears as part of the effective 
potential for the nuclear motion. 


2 
{-5 47 ¥8} 4 [Vion(Ri, Re, --.) 
4+ E(Ri, Re, .. vn} O(Ri, Re, ...) 
= E,®(Ri, Re,...) (1.13) 


The error produced by using the Born—Oppenheimer 
approximation is small when the characteristic frequencies of 
vibration of the ions are much smaller than the characteristic 
electronic frequencies. Generally this implies that the ionic mass 
may be much greater than the electronic mass which is, of 
course, correct. However, if the energy required to excite an 
electron to a higher level is small, implying a small electronic 
frequency, the approximation breaks down. This is the case for 
metals and then the ionic function ® depends on the electron 
co-ordinates as well as the ionic co-ordinates, 

Equation (1.12) describes the electronic motion in a static 
lattice. In practice it is solved assuming both the crystal sym- 
metry and the lattice constant. This implies that the electronic 
potential is a function of the electron co-ordinates only, and so 
the electron energy cannot be used in (1.13) to determine the 
nuclear motion. The equation becomes 
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|-2 EVE + V(r 1 «+ -) [W(t res) 
= Eib(ri, r2,...) (1.14) 


the atomic co-ordinates being given their observed mean 
values. The potential includes all the interactions involving the 
valence electrons and can be written 


Vin, re, ...) = Vari) + Velrs, ve, -- .) 


where Vz(r;) represents the interaction between the ith electron 
and the ion cores and where V, represents the total electron— 
electron interaction. 


§ 1.3 One-electron approximation 


It is extremely difficult to solve atomic or nuclear problems in- 
volving more than one electron. Even in the simplest cases 
complicated numerical procedures have proved necessary. The 
difficulty arises because the electron—electron potential contains 
terms of the form ¢/47re0 | ri — rj | involving the spacial co- 
ordinates of two electrons. In the one-electron formulations 
approximate solutions which treat the electrons separately are 
considered. The Hartree method employs a simple product 
wave function of the form 


b(n, ... rn) = drilri)da(re)...¢n(rw) (1-15) 


It is important that the one-electron functions be normalizable 
and that there be as many as electrons present. For the ground 
state, the variational theorem of quantum mechanics requires 
that the expectation value of the energy be a minimum. Con- 
sequently, the best Hartree solution is that which minimizes 
the integral 


fern, scp Pw), « «5 ¥N)dT (1.16) 
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where # is the Hamiltonian of equation (1.14), subject to the 
restrictive conditions 


[oe (ri) del) =I (1.17) 


This procedure may be carried out by the method of Lagrange 
multipliers and the one-electron wave functions which minimize 
(1.16) and satisfy (1.17) are given by 


KH idi(r1) = Exdi(ri) (1.18) 
where #; is the one-electron Hamiltonian. 


h2 


ns 
2m 


Vi? + Vi(n1) 
+ ay 3* (r2)dy(r2) a ae) 


47eq |r — re | 

The summation in the final term is over all the one-electron 
wave functions except the ith. Each term may be given a 
physical significance. The first term represents the kinetic 
energy of the first electron, the second term represents the inter- 
action between the electron and the ion cores and the final 
term represents the average Coulomb energy of the electron in 
the field produced by all the other electrons. 

The Hamiltonian is not identical for all electrons and the 
eigenfunctions are not orthogonal. 

The simple Hartree method does not entirely neglect the in- 
fluence of one electron upon another. Each electron moves in 
the average field of the remainder. However, there is a serious 
fault in the Hartree wave function. In quantum mechanics it 
is impossible to distinguish between identical particles. This 
means that if any two electrons are interchanged the initial and 
final states of the system must have identical physical pro- 
perties. Consequently the wave function can only change by an 
unimportant phase factor. If the operator Pi interchanges the 
co-ordinates (including spin) of the 7th and sth electrons then 


Prep = etd 
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By repeating this interchange the original wave function is 
restored. Hence it follows that ef = 1 or ef = +1, 


i.e. Puy = tye 


The wave function is either symmetric or antisymmetric with 
respect to permutations of the electron co-ordinates, when spin 
is included. Experimental evidence (Pauli principle) shows that 
electrons, and in fact all particles with half-integral spin, have 
antisymmetric wave functions. Such particles are called 
fermions. 

The Hartree-Fock approach is to recognize that all electrons 
are identical and to start with a determinantal wave function 
satisfying the Pauli principle. 


=k | di(n, 01), dr(re, a2)... dr(rn, on) 


dw(n, 01) éy(ry, on) (1.20) 


(r, «) denote the position and spin co-ordinates and & is a 
normalization constant (4 = 1/4/N! ifthe ¢’s are orthonormal). 
It should be noted that if ¢; = ¢,;, two rows are identical and 
the determinant vanishes. Only one electron may occupy a 
given spin-orbital. Similarly if (#:, 0) = (rj, oj) two columns 
are identical and the determinant again vanishes. Two elec- 
trons with the same spin are never at the same position in 
space. The one-electron functions are again chosen to minimize 
the expectation value of the Hamiltonian (1.14) subject to the 
condition that each is normalized. It is necessary that the 
functions be linearly independent so that it is possible to form 
an orthogonal set from them (Adams, 1961). 

The set of equations which the one-electron wave functions 
must satisfy can be written in the ‘standard form’ 


H bili) = Eide (ri) (1.21) 


where the one-electron Hamiltonians are given by 
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_ _# y* (r2) $;(re) 
i a ie + Vira) + 2. 4meo | ri — Fe dre 


ae _dt(n)b(n)bln)bln) a), 


= ameg| ri — ra | de® (ri) de(r1) (1.22) 


The first two terms correspond precisely to the first two terms 
in the Hartree Hamiltonian (1.19). The summation in the 
third term is over all the ground state electron wave functions 
including the ith and so represents the Coulomb potential 
energy of the electron in the average field of all the electrons 
including itself. The fourth term is summed only over the wave 
functions with spin parallel to the zth and is a new term brought 
about by the use of a determinantal wave function. 

The equations (1.21) have two types of solutions. First of all 
there are the N one-electron wave functions which are occupied 
in the ground state of the system. Apart from these solutions 
there are an infinite number of solutions corresponding to un- 
occupied excited states. All the Hartree-Fock solutions have 
the advantage of being orthogonal to each other, a property the 
Hartree solutions do not possess. 

The fourth exchange term in (1.22) must somehow compen- 
sate for the Coulomb self-interaction of an electron with itself 
which is physically not acceptable. If only the 7th member of 
the summation in the exchange term is used, then the Hartree— 
Fock and the Hartree Hamiltonians are identical. The Hartree 
equations are not unreasonable for free atoms, for they repre- 
sent an electron moving in the field of an ionized atom which is 
physically to be expected. However, if the atoms are assumed to 
form a crystal with a large enough lattice constant so that 
interatomic interactions are small, the Hartree equations in- 
stead of depicting electrons which, at any instant, are essentially 
moving in the field of an ionized atom, describe electrons which 
move in the field of a neutral atom. 

This arises because the term omitted in the Hartree 
Hamiltonian for the self-interaction of an electron is now 
negligible as the wave functions are spread throughout the 
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crystal and not concentrated about an atom. Only an extremely 
small part of the charge density of a crystal wave function is 
associated with a particular atom site. Therefore it may be 
reasonable to use the Hartree equations for free atoms but not 
for crystals. 

The Hartree-Fock equations correct for this fault of the 
simple Hartree solutions. An exchange charge density is defined 


by 


. Parallel 
or ae > $i eee ere ea) 


The Hartree-Fock equations (1.21) are ordinary Schrédinger 
equations of an electron moving in the field of the nuclei and of 
all the electrons including itself, diminished by the field due to 
the exchange charge (1.23). 

If the ith state is occupied then 


[e(rs, ri)dr2 = 1 (1.24) 


and so the ith electron moves in the Coulomb field of all the 
electrons with spin not parallel to it, but the charge of the 
electrons with parallel spins is diminished by one. This charge 
- may be shown to be removed from the vicinity of the electron 
and constitutes the exchange hole. This means that the occupied 
solutions of the Hartree-Fock equations represent electrons 
moving in the field of the ion core and in the Coulomb field of 
(N — 1) electrons where WN is the total number of valence 
electrons. 

If a state is not occupied in the ground state, it is clear that 
the exchange charge 


fore rijdrz = 0 


and so there is no exchange hole and the Coulomb potential 
arises from all the electrons. This is the type of potential an 
extra electron would be expected to experience if added to the 
N electrons already present and if its effect on the ground state 
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functions is ignored. Physically these results are sensible, and 
correct the results of the simple Hartree equations. 

In practice the Hartree-Fock equations are solved in a self- 
consistent manner. An initial set of one-electron wave functions 
is taken and the Hamiltonian #; obtained from equation (1.22). 
Then equations (1.21) are solved and a new set of wave func- 
tions obtained, and the whole process is repeated until the 
initial and final wave functions are in close enough agreement. 

The Hartree-Fock equations have been applied to certain 
atoms but little application to solids has been made as yet due 
to the difficulties involved, although some greatly simplified 
methods have been used. 

Koopmans (1933), has shown that the magnitude of the one- 
electron energies given by the equation (1.21), represent the 
ionization energies of the corresponding electrons when only 
one is excited. These energies are of great importance in the 
theory of electron bands. 

It should be noted that the expectation value of the many- 
electron Hamiltonian (1.14) is less than the sum of the one- 
particle energies FE; as this sum counts each electron—electron 
interaction twice. (For further details of the one-electron 
approximation see Seitz, 1953.) 

As an electron travels through the crystal, it repells other 
electrons and a ‘Coulomb hole’ is formed around it. There is 
correlation between the motion of the electrons. The use of an 
anti-symmetric wave function in the Hartree-Fock approach 
does take account of correlation between electrons of parallel 
spin, to a certain extent, but correlation between electrons of 
anti-parallel spin is ignored. Correlation reduces the total 
energy by keeping electrons apart and following the variational 
principle, the Hartree-Fock approach must be a better approxi- 
mation than the Hartree method. The difference in the many- 
electron energy obtained by the Hartree and Hartree-Fock 
methods is called the exchange energy. This should not be given 
too much physical significance as it is just the difference between 
two mathematical approximations to the truth. 
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The total Hartree-Fock energy Fur is defined by the expecta- 
tion value of the full Hamiltonian for a crystal with fixed ions, 
obtained from a determinantal wave function. This is the ex- 
pectation value of the Hamiltonian in (1.14) together with the 
inter-ionic repulsion terms 


> S>W(R: — Ri) 


t#] 


The ionic terms are added to avoid divergencies. For the ground 
state, the Hartree-Fock energy is an upper bound to the exact 
eigenvalue E. The energy error (E — Er) is called the ‘cor- 
relation energy’ and is a measure of the accuracy of the 
Hartree-Fock approach. The correlation energies for atoms 
are of the order of a few electron volts. 

_ Even though the interaction between two electrons may be 
expected to be of the same magnitude as the interaction be- 
tween an electron and an ion the same distance apart, the one- 
electron approximation is very useful as the following chapters 
will show. Bohm and Pines (1953) have indicated why the one- 
electron approximation gives such good results. They have 
shown that the screening effect of the mobile electrons limits 
the long-range Coulomb interaction between them, but this 
work is beyond the scope of this book. 

The one-electron approximation is not sufficient to fully ex- 
plain such phenomena as cohesive energy or ferromagnetism. 
The full set of solutions of the Hartree—Fock equations (1.21) 
form a complete set of one-electron functions. The totality of 
N x N determinantal wave functions, form a complete set of 
N-electron antisymmetric wave functions and so the exact 
solutions of the many-electron Schrédinger equation (1.14) may 
be expressed as a linear combination of these determinantal 
wave functions. This process is called ‘configuration inter- 
action’. 
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Problems for Chapter | 


1. Write down the Schrédinger quantum operators represent- 
ing the three Cartesian components of angular momentum. 
Verify the commutation relation 


LrLy — LyLy = hf, 


2. An electron is constrained to remain inside a linear box of 
length 2a by walls of infinite potential. Inside the box the 
potential function is zero. Write down the Hamiltonian for the 
particle and find the normalized eigenfunctions for the time- 
independent Schrédinger equation. 

Show that the allowed eigenvalues are given by 


_ neh 


" ~ 32ma® 


where n is restricted to integral values. 

3. Prove that the eigenvalues of a Hermitian operator are 
real. 

4. Prove that the eigenfunctions of a Hermitian operator, 
belonging to different eigenvalues, are orthogonal. 

5. The Hamiltonian operator that describes the motion of an 
electron in a potential V is 


# eve py 


Show that this Hamiltonian is Hermitian. The potential may 
be taken to be real and also to be a function of position only. 

6. The classical Lagrangian for an electron of mass m and 
charge —e in an electro-magnetic field is 


L = 3mV? + ep — e(V. A) 


V is the electron velocity and A and ¢ are the vector and scalar 
potentials of the field. Show that the component of momentum 
conjugate to x; is 
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Derive the classical Hamiltonian. Write down the correspond- 
ing quantum mechanical Hamiltonian. 

(The component of momentum conjugate to x; may be 

: oL 
defined by the Lagrange equation fj = =) 
‘ 

7. One approximation to a metal is to consider the motion of 
a gas of N electrons of charge —e in a volume Q. To ensure the 
neutrality of the system a uniform background of positive 
charge is added, having a charge density equal to the average 
charge density of the electrons. In the Hartree approximation 
the one-electron functions satisfy the set of equations 


f2 € 
ee Lak) ye Re ene ee 
| om i yreo[rn—r| rots 


ai py eet tee boar = Eyhi(ri) 


4meo | ri — rez 


The positive background charge density is p, = Ne/Q. 
Show that the set of plane waves 


$x(r) aa wie oo" 


form a self-consistent solution of these equations if terms of the 
order of ¢2/(Q)* are neglected. 

(This result shows that the energy of the electron gas in this 
approximation is purely kinetic.) 

8. In Problem a, if the origin is taken at one end of the box, 
the electron wave functions are (ignoring spin) 


bn = Jz sin (F x] N= 1,2,... 


Write down the spacial part of the determinantal wave function 
yb(x1, x2) describing two electrons, with identical spin, in the 
two lowest energy levels. If one electron is at the point a/2 in 
the linear box, the probability distribution of the other is 


w(2 xe) W(% xs) Confirm that this distribution is 
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Se Tin wee 
2a” 2a 4/2 a 

Sketch this function and explain how it illustrates the ‘Fermi- 
Hole’ concept. Compare this result with that obtained by the 
simple Hartree product wave function. 
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Free-electron Model 


§ 2.1 The Sommerfeld model 


In the one-electron approximation, the wave equation to be 
solved is 


2 vey + Vin) = Bp (2.1) 


The Hartree-Fock potential V takes account of the interactions 
of the electron with the ionic cores and with the averaged charge 
of the other electrons. Due to the positive point charges at each 
atom site, the potential must possess singularities and the 


ANNAN 


Fic. 2.1 Section through a line of atoms 


Zero | 
potential 


potential plot along a section through a line of atoms is illus- 
trated in Fig. 2.1. 

In 1928 Sommerfeld suggested a model in which the potential 
is taken as a constant inside a metal so that the force acting on 
an electron is zero. Essentially, the metal is replaced by a box 
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containing a non-interacting gas of electrons; the free-electron 
gas. The electrons are not assumed altogether independent as 
they are taken to obey the Pauli principle. The free-electron gas 
does bear some relation to the state of the valence electrons in 
the alkali metals. 

Consider a cube of side a with the origin of the co-ordinate 
system at a corner of the box. Inside the box the potential 
energy of an electron is negative (i.e. attractive) and equal to 
—W and outside the box the potential energy is zero. The 
Schrédinger equation to be solved is 


2/02 a% . 
C2 | 0% | O%p\ (2.2) 
and —= (st T He Oy? - aa) ee otherwise 


The wave functions and their derivatives must be continuous at 
the walls of the box and the bound state wave functions must 
also approach zero at infinity. As W ~ 5 eV for a metal then it 
can be shown that for a box of macroscopic dimensions 
(a ~1 cm) the wave functions are vanishingly small outside 
the box and consequently it is reasonable to use the simplified 
boundary conditions ¢ = 0 at the box walls. 
The equation to be solved is 
Op 


Or dp bas 
aa tae Be + Pre + k%s =o inside the box (2.3) 


where the electron kinetic energy 
Ey = E~— (—W) 
=p (2.4) 
This relation implies that k must be real. For the electron to be 


‘trapped’ in the box the total energy must be negative and so 
Ex < W. 
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The partial differential equation (2.3) is solved by the method 
of separation of variables. By writing 


B(x, I, 2) = X(x) ¥(y)Z(z) (2.5) 
equation (2.3) becomes 
1d?X , 1d?¥ , 1 d2Z 
Yar tyaypt+zaatP=o — (26) 
The first term does not involve _y or z, the second does not in- 
volve x or z and the third does not involve x or_y. This means 
that 


1 d2X 

¥ dee th? =0 

1 d2Y 

¥ ay t fyi =0 

1 d2Z 

7 dz2 + kz? =o (2.7) 
where kz, ky, kz are constants such that 

k2 = ky? +- ky? + kz? (2.8) 

Each of the equations (2.7) has a general solution of the form 
A(x) = A sin kx + B cos kex kz2 >0o (a) 
A(x) = Ax +B kz2 =o (6) (2.9) 


X(x) = Asinh|ke|x+Bcosh|kx|x kez <o (c) 
where A, B are arbitrary constants. 

The only solution which can satisfy the boundary conditions 
is (2.9(a2)) and the wave function (2.5) is then composed of 
terms of the type 
sin sin 
Kyy oe kez (2.10) 
This wave function must vanish at all points on the cube surface 
and the only product satisfying this condition is 


A sin kyx . sin kyy . sin kzz (2.11) 


with eel be ne ke a 
= a” a” a 
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where I, m, n are integers. The electron kinetic energy is 


Ey = peal + m2 + n?) (2.12) 


It is clear from (2.11) that all the linearly independent solutions 
are obtained if l, m, n take positive values only. Observe that it 
is the application of the boundary conditions that quantizes the 
energy. The energy levels lie very close together (~ 10-4 eV 
apart) and for most purposes may be taken as continuous; the 
energy levels are said to be quasi-continuous. The total electron 
energy is 


E= 


2 
mann +m +n) —W (2.19) 


The number of electron states (neglecting spin) with kinetic 
energy less than some given energy £;, is easily obtained. It is 
the number of integral points within the positive octant of the 
sphere of radius r where 

he 
Ey = 


Smae* 27 (2.14) 


For large values of r the number of the points is simply the 
volume of the octant 

1) 248 

(fer ) 


So the number of electron states with energy less than £, is 
8m 
M(Ex) =Q aa si ) Eg (2.15) 


where Q = a? is the volume of the box. If spin is considered the 
number of states is twice that given by equation (2.15). The 
number of states per unit energy range is called the aan of 
states N(E) and is related to M(E) by 


iz N(E)dE = M(Ex) (2.16) 
N(Ex) = (See (2.17) 
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Hence 


N(Ex) = Q “(Gr) Ee (2.18) 
Note that the number of states per unit energy per unit volume 
is independent of the dimensions of the box. Again, when spin 
is included the density of states is twice that given by (2.18). 
A more complete treatment shows that (2.18) is valid when 
non-cubic volumes are considered. 

The solution of the wave equation given by (2.11) is a con- 
sequence of the particular boundary conditions imposed. Sup- 
pose that instead of a box in which the electrons are trapped, 
there is an infinity of such boxes stacked together and all 
identical, The appropriate boundary conditions are the cyclic 
or Born von Karman conditions. 


p(x + 4, y, z) —. p(x, + 4, z) a (x, I, das a) 
= (x,y, 2) (2.19) 


The normalized solutions of the wave equation (2.3) are now 


h(x, Dds z) = Ta et (Kast + kyy + kez) (2.20) 


kz = (*2), ky = m(22), Ke = n( ==) (2.21) 


where /, m, n are integers. Here the integers take all values, both 
positive and negative, as the corresponding solutions are 
linearly independent. 

The electron kinetic energy is now 


2 
2ma2 


Ey = (12 + m2 + n?) (2.22) | 
The formulae (2.15) and (2.18) still apply. This may at first 
seem surprising because of the factor-four difference in the 
right-hand sides of equations (2.12) and (2.22). However, when 
it is remembered that both positive and negative values of 
l, m, n are now to be considered, this result quickly follows. 
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(Note — when the radius of a sphere is halved the volume is 
decreased by a factor of eight.) 

Suppose there are N ‘non-interacting’ electrons contained in 
the box at absolute zero temperature. If the electrons obeyed 
classical rules then all of them would be in the lowest level. But 
this is not the case as electrons must obey the Pauli principle 
which demands that no two electrons possess the same set of 
four quantum numbers. Three of these numbers are /, m, n of 
(2.22) (or (2.12)) and the fourth is spin. Each pair of states 
defined by /, m, n may be occupied by two electrons. At absolute 
zero the electrons occupy energy levels from zero to a maximum 


value Em defined by 
N aw{[8m\i_, ; 
3 = 2G Fr) Eat 


: _ (A*\/3\', 

i.e. En = (sz)(2) n (2.23) 
where n is the number of electrons per unit volume. The value 
of Em which corresponds to the valence electron density in 


metals is of the order of a few electron volts. The average 
kinetic energy of the electrons is 


p= [ E,N(Ex)d/ [ *" N(Ey)dE 


= SEm (2.24) 


This result will be referred to again in the discussion of the co- 
hesive energy of metals in Chapter 7. 

At absolute zero there is a bounding energy below which all 
the electron states are filled and above which all states are un- 
occupied. As the temperature is increased some electrons are 
excited into higher states and the boundary between occupied 
and unoccupied states becomes blurred. In this case it is appro- 
priate to talk about the probability of a state being filled and 
resort must be made to statistical mechanics. 
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§ 2.2 Fermi—Dirac statistics 


The fundamental result of statistical mechanics is that the prob- 
ability of the system being in an allowed state with energy £ is 
proportional to exp(—£/kT) where & is the Boltzmann constant 
and 7 is the absolute temperature. This relation holds in 
quantum mechanics as well as in classical mechanics. However, 
it is only useful in quantum theory if the correct many-electron 
eigenstates are known. In the one-electron approximation the 
system is described in terms of the single electron eigenstates. 
The many-electron wave function must change sign under the 
interchange of any two electrons and this imposes a restriction 
on the allowed eigenstates of the complete system. Only two 
electrons can occupy the lowest pair of one-electron states 
(specified by J, m, mn) and it is meaningless to talk about 
the probability of all the electrons (assumed more than 
two) occupying this lowest energy as that situation can not 
exist. 

The Fermi—Dirac distribution law gives the probability that 
a single-electron state is occupied, remembering that the many- 
electron wave function must be antisymmetric with respect to 
interchange of electrons. Clearly the Fermi—Dirac law is not as 
fundamental as the Boltzmann law as it depends on a Hartree— 
Fock type of approximate wave function, but even so it is im- 
portant and is derived below. (The reader is referred to Rush- 
brooke’s book (1949) for a background to statistical mechanics 
and to Schrédinger’s text (1946) for a more detailed study of 
quantum statistics.) 

It is assumed that the electrons form a gas of ‘independent’, 
structureless particles. The electron interactions are either 
ignored completely or approximately taken into account in an 
average way by the Hartree-Fock method. The electrons are 
Fermi particles; they are identical and at most one can occupy 
a completely specified single-electron energy state. If the energy 
levels are extremely close together and almost form a con- 
tinuum then it is possible to group the energy levels in bundles 
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and give each bundle a relevant average energy. This is cer- 
tainly the case for the levels in a box of non-microscopic di- 
mensions and is also true in crystals. 

Suppose the 7th bundle contains w; states with an average 
energy £;. It is assumed that w; and the number of electrons 7; 
in the bundle are both large. Clearly nj; < w; as there is at most 
one electron in each state. ‘The number of distinguishable ways 
of distributing the electrons among the states in the bundle is 


w,! 


ny! (wr =~ n)! (2.25) 


In deriving equation (2.25) it must be remembered that the 
electrons are identical and so distributions differing only by an 
interchange of electrons are equivalent. For a given set of 
distribution numbers m, the total number of distinguishable 
ways of arranging all the electrons among the states in all the 
bundles is 


t 
‘=TI Wt. 


; ny! (wy _ n;)! (2.26) 


There are restrictions placed on (2.26). First, the sum of the 
distribution numbers must add up to the total number of 
electrons and second, the sum of the electron energies must 
equal the total energy E. 


2m =N (a) 


Yuki =E (b) en 

t 
This latter condition requires further consideration. It has 
already been stated in Chapter 1 that the total energy of the 
system is not simply a sum over the one-electron Hartree-Fock 
energies. However, it can be shown that if the excited many- 
electron state is close to the ground state (i.e. only a few 
electrons are raised by a small amount) then the results ob- 
tained using the second of equations (2.27) are essentially 
correct. This is true for thermal excitations at ordinary tem- 
peratures (see Slater, 1956). 
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There are other sets of distribution numbers, apart from 
those in (2.26) which satisfy (2.27) and so give rise to other 
distinguishable arrangements of electrons. All such possible 
arrangements for a given total number of electrons and a given 
energy are taken to be equally probable. However, for suf- 
ficiently large N it can be shown that the observed averaged 
properties of the system are identical with those obtained from 
the most probable set of distribution numbers found by maxi- 
mizing (2.26) subject to the conditions (2.27). In fact it is more 
convenient to maximize In ¢. 


In ¢ = > [In (wi!) — In (mi!) — In {(wi — ma) 9] 
w, and m are both large and if (w; — m) is also large then 
Stirling’s approximation for the factorial may be used. 
In (x!) =x In x — x 


A rigorous treatment would require some justification of these 
approximations, but it will suffice to note here that the final 
result obtained may be verified by more powerful methods. 


In ¢ = > [wy In wy; — m4 In m — (wi — m) In (wi — )] 


The method of Lagrange multipliers is used to maximize In ¢ 
subject to the restrictions (2.27). Then 


io [inet (dm — N) +8 (3mBi — £)] =o forall 


where «, 8 are multipliers, 
i.e. —Inmg —1 +n (wj—-m)+1+a+ 8k; =o 
at (2.28) 


and mM = Sake Ly 


a and 8 are of course determined by (2.28) together with the 
restrictive conditions (2.27). The entropy S is a function of the 
internal energy E, the number of electrons N and the volume 2 
and is given by 

S = k In tmax. 
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where & is Boltzmann’s constant. From the thermodynamic 


relations 
(= ae (Sx) cee 
5E ya TT ON/z.qa I 


it is possible to identify B as —1/kT and « as £/kT were ¢ is the 
chemical potential per electron. In solid state theory ¢ is 
generally referred to as the Fermi energy and denoted by Er. 
The distribution function f = n;/w; gives the probability that a 
given state is occupied and 


I 
J 4) = gamer Ty (2.29) 


This is the Fermi—Dirac distribution law. At T=0, f=1 
for E; < Er and f =o for E; > Ey. All states with energies 
less than Er are filled whereas those with energies greater than 


F(Eé) 
1-0 


0°5 


Ee 3 
AT 


Fic. 2.2 Fermi distribution 


Er are empty. As the temperature is increased from absolute 
zero the distribution function departs from this step character. 
Some electrons in the states immediately below Er are excited 
into levels just above. The only region where f is appreciably 
different from one or zero is in the neighbourhood of Er over a 
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range of energies of the order of kT. The Fermi energy is itself a 
function of temperature but when kT << Ep(T) it can be 
shown that it changes only slightly from the value at T = o. 
At all temperatures f has the value $ when £; = Ey. For high 
electron energies where E; >> Ep, f behaves like the high 
energy tail of the classical Boltzmann distribution. 

The criterion for quantum statistics to be necessary is that 
the Pauli principle be important. This occurs when there are 
‘too many’ electrons for the available energy levels. If classical 
statistics are obeyed then almost all the electrons will have. 
energies up to ~kT' and at room temperature this is ~Z; eV. 
Silver is a ‘typical’ monovalent metal where one electron per 
atom is denoted to a sea of loosely bound conduction electrons. 
In one mole of silver, the volume of which is ~10 cc, the 
number of electron states with energies less than -; eV is 
~1020 (see 2.15). But one mole of silver contains 1023 conduc- 
tion electrons. There are far more electrons than states available 
and certainly quantum statistics must be used. This is the case 
for metals where to accommodate all electrons, states up to 
about 5 eV are filled. More formally, quantum statistics is 
necessary when kT << Ey and then the electron gas is said to 
be degenerate. A degeneracy temperature may be defined by 
kTo = Er(o) and the condition for degeneracy is T << To. 
For metals To ~ 104°K, 


§ 2.3 Variation of the Fermi energy with temperature 


The Fermi energy is not constant but decreases slowly as the 
temperature rises. The relation is obtained below. 

An important property of the Fermi—Dirac function is that 
its derivative differs from zero only in the neighbourhood of the 
Fermi energy. 

With this in mind consider the integral 


Te [ “S(E)-_G(B)AE (2.90) 
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where f (£) is the Fermi function and G(Z) is any function such 
that G——> 0 as E ——> o. Integrating by parts 


r=[s@owm | — Power (2.31) 


The properties of f(£) and G(£) make the first term on the 
r.h.s. vanish at the upper and lower limits respectively. If G(£) 


df 
dé 
O E; E 
____> 
kT 


Fic. 2.3 Derivative of the Fermi function 


does not possess a singularity at Ey it may be expanded by 
Taylor series and within its radius of convergence. 


G(E) = G(Er) + (E — Eg)G' (Er) 
+3(E — Er)?G" (Er) +... (2.32) 
Then 


 — f “[G(Er) + (E — Ep)G’(Ep) 
+ 3(B — Ep)*G" (Er) +.. J f'(E)dE (2.33) 
Assuming term by term integration 


= — IoG(Er) — 1,G’(Er) — L2G" (Er) —... (2.34) 
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= -—!I 


with ie : “fi (E)dE = [Fe) | 


ie i "(E — Ep) f'(E)dE = 0 


Le 4/ (E — Er)2f'(E)dE #0 (2.35) 
0 
To calculate Le it is useful to substitute x = (£ — Epr)/kT and 
then 


ly =— dx 


(KT)? { e neu 

2 Jeu (1 + e)* 
In metals kT << Ey and the lower limit may be replaced by 
— oo and so (see any standard table of definite integrals) 


Ly = —Z(kT)? (2.36) 
The integral J is 
2 
I = G(Er) + %(kT)?G" (Ep) +... (2.37) 


From this very useful result it is possible to show how the Fermi 
energy varies with temperature. To do this put G(£) equal to 
the total number of electron states (including spin) with energy 
less than E. 


E 
G(E) = 2 i N(E)dE so that G'(E) = 2N(E) 
0 
With this choice J represents the total number of electrons N. 
E. 2 
Nis i "N(E)dE + 7 (kT) 2N' (Er) +... (2.38) 
0 
N is also given by 
EF(o) 
N=2 | N(E)dE (2.39) 
0 


where Ey(o) is the Fermi level at absolute zero. From these last 
two equations 


EF are 
[0 Nae = —2(kT)N'(E®) 
Er(o) 
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As Er — Ep(o) is small then 
2 
[Er — Ep(0)]N(Ex(0)) = —%(kT)?2N’ (Er(0)) 


, 7 d 

ie. Ey = Ep(0) —-(k T)( + in N) (2.40) 
For free electrons N is given by (2.18) and Ep(o) is the maxi- 
mum electron energy (2.23). 


27k T\2 
Ep = En - “(z-) ] (2.41) 
The Fermi energy falls as the temperature increases but the 


change is only of second order in kT/Ey. At room temperature 
KT/Ep ~ 10-3 for metals. 


§ 2.4 Electron specific heat 


If electrons were classical particles then an energy of kT/2 
would be associated with each degree of freedom and the 
interna] energy of a gas of N non-interacting electrons would be 


U = 2 NkT 
The heat capacity is defined by 


and so Cy = 3Nk/2. This constant value is not observed and 
experiment shows that C, decreases linearly to zero as T ——> o. 
This latter effect is due to the fact that electrons are Fermi 
particles. 

The total kinetic energy (the energy zero is chosen here to be 
the base of the Fermi distribution) of the electrons in the 
volume Q is 


Us i “Ef (E)N(E)dE (2.42) 
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The factor two arises as each state may contain two electrons. 
Define 


E 
G(E) =2 [ EN(E)dE 
0 
and from (2.37) at sufficiently low temperatures 
Uxef EN(E)dE + (kT)? S(EN(E 
= 2 VEN(E)AE + = (kT)*55(EN(E)) or (2-43) 
The first term on the r.h.s. may be expanded by Taylor series. 
EP Er(o) 
2 | EN(E)dE~2f[  EN(E)dE 
0 


+ 2[Er — Er(0)|Er(0) N(Er(0)) 


The first term on the right is the kinetic energy Up of the 
electrons at absolute zero and the second is (2.40) 


2 
—7 (#7) *Ex(0) N’(Ex(0)) 
The total kinetic energy at temperature T is then 
2 
U = Uy + S(kT)®N(Ex(0)) (2.44) 


and the heat capacity of the electrons is 
Cy = §r*k*N(Er(0)) T (2.45) 


This expression is a good approximation for most metals at all 
temperatures for which they are in the solid state. This linear 
contribution of the electrons to the heat capacity is only ob- 
served at very low temperatures since at higher temperatures it 
is swamped by the contribution from the lattice vibrations 
which tend to zero as 73. For free electrons N(£) is given by 
(2.18) and this result together with (2.23) may be used to 
express Cy in terms of the total number of electrons. 


_3 =) 
Cpe 3.ye( * ZT (2.46) 
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This model can be used to represent a monovalent metal if N is 
made equal to the number of atoms in the volume. (For one 
mole of the metal Nk = R, the gas constant). Equation (2.46) 
is in quite good agreement with experiment. If there is some 
doubt about the number of free electrons per atom measure- 
ments of Cy, can be used to give information about the density of 
states at the Fermi level. 

Comparison of (2.46) and the classical value suggest that as 
KT << Em the number of electrons that may be considered 
free is much less than the total number of electrons. This result 
is a direct consequence of the Fermi—Dirac statistics. Only those 
electrons near the Fermi level are ‘free’ as they alone are able 
to be excited into higher energy states without an improbable 
energy jump. 


“2. Weak spin paramagnetism 


The magnetic properties of an atom are due to the orbital 
motion of the electrons and to their spins. The nuclear magnetic 
moment is negligible in comparison. The magnetic materials 
that occur in nature may be divided into four basic types: 


(a) ferromagnetic 
(5) antiferromagnetic 
(c) diamagnetic 


(d) paramagnetic \ include most metals 


These will now be discussed very briefly. 

Below a certain temperature (Curie temperature) ferro- 
magnetic substances are spontaneously magnetized and this is 
caused by an ordered coupling of parallel atomic magnetic 
moments. Examples are iron and nickel. In antiferromagnetic 
materials there is also an ordered coupling of atomic magnetic 
moments but arranged in antiparallel pairs so that there is no 
permanent macroscopic magnetic moment. There is also an 
intermediate class, the ferrimagnetic materials where the 
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atomic moments are antiparallel but not equal, so that there is 
a permanent magnetic moment. 

For both diamagnetic and paramagnetic substances the 
magnetic moment per unit volume M is proportional to the 
applied magnetic field H. 

M = x 


The constant of proportionality y is called the susceptibility per 
unit volume. 

_ If y is negative, the material is diamagnetic and the magnet- 
ization is in the opposite direction to the field. This is an 
example of Lenz’s law and may be explained by the induced 
Larmor procession of both the free and ionic electrons. Dia- 
magnetism may be expected to occur in saturated materials 
whose constituents do not possess permanent magnetic mo- 
ments. Induced magnetization of this type is generally very 
small (| x | ~ 10-5) and temperature independent. 

If y is positive, the induced magnetization is parallel to the 
field and the substance is said to be paramagnetic. A few sub- 
stances are strongly paramagnetic (x ~ 107%) and this is due to 
the orientation of the permanent magnetic moments of the 
constituents along the field. This can only happen if the dipoles 
can orient themselves freely in an applied field. Thermal 
motion acts to oppose this effect and it can be shown that the 
susceptibility obeys the Langevin formula and is inversely 
proportional to the absolute temperature 


xX = Nme2yo/3kT (2.47) 


N is the number of magnetic cores per unit volume and m,- the 
magnetic moment of the core. Crystals containing rare earth 
ions exhibit this type of magnetization. Their paramagnetism is 
due to the magnetic moments of the inner incomplete shell. 
Many metals have a weak temperature independent para- 
magnetic susceptibility of the order of 10-6. This is due to the 
spin of the conduction electrons. The Langevin formula (2.47) 
predicts that the susceptibility be temperature dependent and 
much larger than that observed at room temperature. The 
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reason why the Langevin theory fails is that the electrons can 
not orient themselves freely in an applied field but must obey 
Fermi-—Dirac statistics. Pauli (1927) gave the correct theory. 


Paramagnetism of free electrons 

The electron has a spin magnetic moment ms equal to one Bohr 
magneton. 

eh 
2m 


(2.48) 


In an external magnetic field H, the energy of a static electron 
can have one of two possible values. 


+msB B=pH (2.49) 
(uw is the permeability) 


8 = 


The minus and plus sign apply when the electron spin is parallel 
and antiparallel to the field respectively. If there was nothing to 
prevent it all the antiparallel electrons would orientate them- 
selves parallel to the field, each transferring an energy 2msB to 
the lattice vibrations. In a free-electron gas this is not possible 
for most electrons as the states with parallel spins are already 
filled. However, at non-zero temperatures, there is always a 
small number of electrons in excited states and so electrons near 
the Fermi surface may be able to change their spin direction. © 

If it is assumed that the spacial movement of the electrons is 
unaffected by the magnetic field then the two degenerate 
electron states with initial energy £ split into states with energies 


E += mB (2.50) 


The number of electron states per unit volume with spins 
parallel to the field in the energy range E to FE + dE is 


5 N(E + m,B)dE (2.51) 


Q is the volume in which the electrons are contained and N(£) 
denotes the density of states (neglecting spin (2.18)) in the 
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absence of the field. The number of states with spin anti- 
parallel to the field is 


a N(E — m,B)dE (2.52) 
The magnetic moment per unit volume is then 
M =o PME + mB) — N(E — msB)|f(E)dE (2.53) 


where /(£) is the Fermi—Dirac function (2.29). 
In practice the magnetic field is very weak and 


N(E + mB) — N(E — mB) = amsB (2.54) 


The susceptibility per unit volume is then 


x= "f (el CRE as mM ~ po (2.55) 


From (2.37), for sufficiently low temperatures 


a d2N 
AEN(Er)uo + wore RT)'| See], (256) 


— 


For all normal temperatures the first term is dominant and for 
free electrons N is given by (2.18). 


— 2 E,° (2°57) 


where n is the number of electrons per unit volume. If 
n ~ 1023 cm-3 and Em ~5 eV then y ~ 10-6 in quite good 
agreement with the alkalis. It should be noted that the dia- 
magnetic effects of the magnetic field on the electron motion 
have been ignored. Landau (1930) has shown that for free 
electrons this susceptibility is equal to —4 of the paramagnetic 
susceptibility obtained above. The total susceptibility of the 
free electron gas 1s 


ms2n 
X= E Lo (2.58) 
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Problems for Chapter 2 


1. The potential in a one-dimensional problem is defined by 
V=—W |x| <a V=o0 |x| >a, W>0 


Solve the time-independent Schrédinger equation using the 
correct boundary conditions 


(a) ¢——> 0 as | x | ——> ©. (Bound state.) 


(b) % and vt are continuous at | x | = a. 


Show that if , / am] EI . 4 >> | then the simplified boundary 


conditions, y =o at|x| = a may be used. Sensible values for 
a metal are | E| ~5 eV and a~1 cm. Confirm that in this 
case the simplified boundary conditions are valid. 
2. The potential in a three-dimensional problem is defined 
by 
V=-—-W r<a V=o0r>a Wr>o 


(r is a spherical polar co-ordinate). 
The eigenfunctions satisfy the conditions 
(a) & is finite at r = o. 
(6) ¢ and id are continuous at r = a. 
(c) ~——> 0 as r-—> o. (Bound state.) 


Show that the spherically symmetric bound solutions (E < 0) 
are of the form 


A [—2mE 
eats Sick r>a a= Re 

Bos 2m 
~ = — sin Br r<a B= [z(W + E) 


Verify that a necessary condition for a non-trivial solution is 


aa = —fa cot (Ba) 
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a and 8 must also satisfy 
2m 
(aa)? + (Ba)? = 2% Wa? 


Illustrate these conditions graphically. 
Explain why bound solutions exist.only if 


(ae seals 
Re 3 


| The spherically symmetric Laplacian is 5 Zr =) | 


g. An electron is confined in a cylinder of radius a and length 
I. The potential inside the cylinder is —W (W > 0) and is zero 
outside. If the co-ordinate axes are chosen so that the z-axis 
coincides with the cylinder axis and the origin is at the centre of 
the lower plane end then Schrédinger’s equation is, in 
cylindrical polars (p, ¢, 2) 

h?T1 of 0 1 02 02 
— Fal pap Pap) + preg + Gea) — 1 + EM =o 

p<a o<z<l 0<¢ < 27 


If | Z| is sufficiently large the boundary conditions reduce to 
y =o on the cylinder boundary. By the method of separation 
of variables show that the solutions finite and single-valued 
within the cylinder are of the form 


cos rp 
sin rh 


_ at 
yb = A sin Te: Fr(ap) 

where n is a positive integer, 7 is a positive integer or zero and 

Fr(aa) = o. Verify that the corresponding eigenvalues are 


h2 (i 


em\ [2 


+at)—W 


(The Bessel function 7,(x) is a solution of 
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4. The entropy of a system is a function of the total energy E, 
the number of constituents N and the volume Q. Assuming that 
S = k In tmax, show that for the free-electron gas 

= —Na — EB + Swi In (1 + e*eFF) 


From the thermodynamic relations 


as ae: (sy) - ie 
(Foy = T \ON/eq  T 


and the restrictive conditions (2.27) show that 


I 
= ET a= ep 
(a, 6 must be regarded as functions of Z, N, Q.) 

5. It has been stated in the text (2.2) that the probability ofa 
system being in a state with energy £ is proportional to 
exp (—£/kT). Illustrate how this is consistent with the Fermi— 
Dirac distribution (2.29). To do this consider two energy states 
of the entire system which differ only by the excitation of a 
single electron. Initially the ith one-electron state is occupied 
and the jth state empty. In the final configuration the jth state 
is filled and the zth state empty, all the other electrons re- 
maining unchanged. From the Fermi—Dirac law show that the 
ratio of the probability of the final state to that of the initial 


state is 
exp [(E: — E) [KT] 


(Note: the probability of the zth state not being filled is 
[1 —f(F:)].) 
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The Importance of Crystal Symmetry 


§ 3.1 The Bravais lattice 


The fundamental property intrinsic to crystalline solids is their 
degree of order, whereas gases and liquids are in a state of dis- 
order. The most important manifestation of this order is trans- 
lational symmetry. A crystal has a periodic structure. There are, 
of course, deviations in real crystals, both at the surface, where 
the periodicity condition obviously breaks down, and also due 
to dislocations, impurities, and other imperfections in the 
interior. In most cases these can be considered as perturbations 
of an ideal crystal. 

Solids possess symmetry properties apart from translational 
symmetry. These are the reflection planes and rotation axes 
which help to determine the external appearance of a crystal. 
From the crystallographic viewpoint, the rotation-reflection 
symmetry elements are most useful in the classification of 
crystals, but to the solid-state theorist the translational sym- 
metry property is of greater consequence. 

As crystals possess translational symmetry, it must be possible 
to define a minimum periodic displacement in any direction. 
Any two points separated by such a displacement must be 
equivalent in every respect. For a set of identical points placed 
at equal intervals along a straight line, the distance between 
any two consecutive points is the minimum repeating distance. 
In the case of a multi-dimensional periodic array of points, the 
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minimum periodic displacement is very long for most directions. 
However, for a two-dimensional array, it is possible to define 
two basis vectors, a1 and az, so that all the equivalent points are 
given by 

Rn = mai + neaz (3.1) 


where 71, nz are integers. Such an array is called a Bravais net. 
Similarly for a three-dimensional array of points, it 1s possible to 
choose three basis vectors, a1, @2, @3 so that all the equivalent 
points in the Bravais lattice are given by 


Rn = a1 + nea2 + nsas (3.2) 


where 71, no and ng are integers. 


e e es e 
e e e 
r o--—-~--—~ e 
af 
—* e 
! 
e e e e 


Fic. 3.1 Rhombic Bravais net 


The basis vectors for a lattice are not unique, but they must 
satisfy the requirement that all the equivalent points in the 
lattice are given by equation (3.2) with a suitable choice for the 
integers n;. In addition, all integral choices of m must represent 
a lattice point. A suitable set of basis vectors may be obtained 
as follows: 


ai is taken as the shortest possible periodic length in the 
lattice; 

az is the next shortest length not parallel to a1; and 

az is the next shortest length not co-planor with ai and az. 


The three basis vectors define a parallelepiped which is 
called a primitive cell. The primitive cell has a volume given by 
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the triple product a, . (az (A as) and is the smallest unit which 
will build up the lattice by periodic repetition. Different choices 
for the basis vectors yield different primitive cells but they all 
have the same volume. This can be understood if it is observed 
that one lattice point is associated with each primitive cell, and 
so for a given density of lattice points the primitive cell volume 
is fixed. 

A Bravais lattice has other symmetry properties apart from 
translational symmetry. A centre of inversion is a point such 
that the lattice is invariant under the operator r —> —r where r 


Fic. 3.2 Square Bravais net with line of symmetry 


is the vector displacement of any point in the lattice referred to 
the inversion centre. A plane of symmetry is a plane in the 
lattice such that the lattice points on one side of the plane are 
mirror images of those on the other. 

An axis of rotational symmetry is an axis such that the lattice 
takes up an identical position on rotation through a suitable 
angle about it. The angle of the rotation must be of the form 
2zj/n where n is an integer and such an axis is called an n-fold 
axis, A Bravais lattice may also possess a rotation-inversion axis 
if it is brought into self-coincidence by a rotation followed by an 
inversion. A two-fold inversion axis is equivalent to a mirror 
plane perpendicular to the axis (Fig. 3.3). 


43 


SOLID STATE PHYSICS 


An infinite periodic array of points is invariant under a 
translational of the form (3.2). This condition leads to certain 
theorems regarding the presence of rotation-reflection sym- 
metry elements. For example, it can be shown that a Bravais 
lattice has inversion symmetry about each lattice point and also 
about the midpoints of the lines joining two lattice points. 


1 Inversion axis 
x P 


AN é Mirror 


\ plane 
| 
| ‘ 
| 
| 
| 


Fic. 3.3. A two-fold inversion axis is equivalent to a 
perpendicular mirror plane 


However, actual crystal structures may lack this symmetry due 
to the distribution of atoms associated with each lattice point. 
An important result is that a Bravais lattice can only have an 
axis of two-, three-, four- or six-fold symmetry. There are no 
five-fold axes of symmetry in lattices although this restriction 
does not apply to physical systems lacking translational sym- 
metry such as molecules. 
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There are fourteen possible ways of positioning structureless 
points in space such that all points have identical surroundings. 
These are the fourteen Bravais lattices. Each of these fourteen 


Rhombohedral primitive celi 


Fic. 3.4. Face-centred cubic lattice 


Bravais lattices belongs to one of seven crystal systems which are 
characterized by the groups of symmetry elements associated 
with the lattice points. There are three Bravais lattices belong- 
ing to the cubic system, each of which has the same collection of 
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symmetry elements at the lattice points. Table (3.1) lists the 
different lattices and the corresponding crystal systems. 

The primitive cell has been defined as the smallest repeating 
volume in the lattice, and has only one lattice point associated 
with it. This parallelepiped often disguises the lattice symmetry 


TABLE 3.1 
Crystal system Bravais lattice 
Triclinic Simple 
Monoclinic Simple 
Base-centred 
Orthorhombic Simple 
Base-centred 
Face-centred 
Body-centred 
Tetragonal Simple 
Body-centred 
Hexagonal Simple 
Rhombohedral Simple 
Cubic Simple 


Face-centred 
Body-centred 


and it is sometimes useful to introduce a larger cell which 
illustrates the symmetry more clearly. Such a cell will be 
associated with more than one lattice point and is called a unit 
cell. Figure 3.4 shows a primitive cell for the face-centred cubic 
lattice and also the corresponding cubic unit cell four-times the 
size, containing four lattice points. 


§ 3.2 Symmetry operators 


When a body possesses symmetry, it must be possible to bring it 
into self-coincidence by some symmetry operations, such as a 
rotation or reflection. A physical rotation of a body through an 
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angle « in the clockwise direction about an axis, say the 
%3-axis, is equivalent to a rotation of the axes system through 
the same angle in the anti-clockwise direction, in the sense that 
the relation between the new and old co-ordinates of any point 
in the body is the same in both cases (Fig. 3.5). 
In each case | 
x1 = x1’ Cos a — x2’ Sin 
xq == %1' sin « + x2" Cos & (3.3) 


x3 == x3" 


*1, x2, x3 and x1’, x2’, x3’ represent the old and new co-ordinates 
respectively. 

In thinking about the symmetry properties of molecules and 
crystals, it is very useful to interpret the symmetry operations 
‘actively’ in the sense of actual rotations of the physical system, 
but mathematically it is far more convenient to use the 
‘passive’ viewpoint of axes rotation. In crystalline solids, the 
most general symmetry operation is a rotation, reflection or 
inversion about the origin, followed by a translation. Following 
the passive viewpoint, the new co-ordinates *1’, x2’, x3’ of a point, 
are given in terms of the old co-ordinates x1, x2, x3 by a 
linear relation of the form 


X= H11%1 + a12%2 + a1gv3 + fy 
Meit1 + aeere + aeax3 + be (3.4) 
x3 == ag1%1 -+ agex2 + agra -+ fs 


s 
I 


The matrix (a) represents the rotational part of the operator 
and & represents the 7th component of the translation. In matrix 
formulation 


x1’ 11 O12 a13\ / x1 ty 
xe’ | =| we. ae aes jl x2 | +] te (3.5) 
x3" 31 a2 ass / \xX3 tg 

or x’ =ax +t (3.6) 


The notations (3.4), (3.5) and (3.6) will be used as required. 
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For a pure rotation or reflection about the origin, the trans- 
lation components are zero. Then the distance of a point from 
the origin remains constant and 


r% == x42 4 xe? + xg2 = xy'2 + x02 + 3/2 (3.7) 
The matrix (a) is an orthogonal matrix and the determinant 
Jo] = 1 
A convenient notation for representing transformations of the 
type (3.5) is that introduced by Seitz. The symbol {a | t} denotes 


a rotation (inversion or reflection), about the origin, followed 
by a translation ft. {e | 0} denotes the identity operator. 


xy’ Io 0O\/*x1 
i.e. xe }=[([0 1: O]f xe (3.8) 
x3" 0 o 1/\x3 


{e | t} represents a pure translation and {« | 0} denotes a rotation 
about the origin. Two transformations {« | t}, {8 | t’} taken in 
succession, may be represented by (see (3.6)) 
x = a(Bx +t’) +t (3-9) 

i.e. {a | £}{8 | t’} = {a8 | at’ + 2} (3.10) 

The inverse of a transformation is that which returns the 
crystal to its origina] situation. This means that the product of 
an operator, with its inverse, must be equivalent to the identity 
operator. The inverse of {« | tf} written {a | t}-1 is {a-1 | —a-1¢} 
where a-! is the inverse matrix of «a The inverse of the in- 
version operator {I | 0} is itself. The same is true for the identity 
operator {e | o}. 

The effect of any operator on a function is to transform it into 


some other function. For example, the differential operator s 


changes x? into 2x. It is necessary to define the effect of a Seitz 
operator on any arbitrary function. The operator is taken to be 
a substitutional one and its action on a function f(r) is 


{a | Of(r) =f ({o | 4) (3.11) 


With this definition, the value of the function at each point in 


49 


SOLID STATE PHYSICS 


space is left unaltered under the rotation or translation of the 
co-ordinate axes. 


As an example illustrating equation (3.11), consider the two 
functions 


fz = cos 9, py = sin @ 


Old axes 
‘function p, 


New axes 
:function-p, 


FIG. 3.7 


When the co-ordinate axes are rotated through 7/2 radians 
about the z-axis in the anti-clockwise sense, the function pz 
becomes — fy if the value of the function at each point in space 
remains unchanged. 
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The symbol representing the axis rotation may be written as 
{42 | 0} and it has been shown above that 


{42 | O}pe = —py 


This result also follows from equation (3.11). The inverse of 
the operator {4, | 0} represents a clockwise rotation of the axes 
through 7/2 radians and 


{42 | O}f2 = cos [{42 | 0}-18] = cos (9 ++ =) = —sin 8 
Substitutional operators are linear 


{a | t}Lapa(r) + bysa(r)] = aya({o | thr) + bybo({a | 2}-1n) 
= a{a. | thha(r) + d{a | thp2(r), (3.12) 

where a, 5 are complex constants. 

The possession of symmetry by a body has important implica- 
tions for the solutions of Schrédinger’s equation. 

For any physical system studied within the one-electron 
approximation, the Hamiltonian for an electron may be 
written : 


f2 
H (x1, v2, x3) 7s * ie a V (x1, 2, x3) (3.13) 


where spin-dependent terms are omitted. If the system under- 
goes some reflection, rotation or translation, then the old co- 
ordinates x1, x2, x3 must be replaced by new co-ordinates 
x1’, xa, xg" as defined by an expression of the form (3.4). It is 
shown in standard textbooks on the calculus, that the Lap- 
lacian operator V? is invariant under such a linear transforma- 
tion. 
é 02 02 02 02 02 02 
Me Batt Gaat t Bast ~ Gaye b Gage + Gage (3-14) 
If the change in co-ordinates corresponds to a ‘symmetry 
operation ’, then the system is brought into self-coincidence by 
such an operation, and it is quite clear that the potential 
function is also invariant under the transformation. 
As an illustration, consider the potential which an electron 
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experiences when moving in the field of two protons. This 
occurs in the hydrogen molecular ion, 

One of the symmetry operations which brings the physical 
system into self-coincidence, is a rotation through 180° about 
any axis perpendicular to, and bisecting, the line joining the 
two protons. Suppose the arbitrary point P is rotated into P’ by 


P(x, X51 X3) 


PUK Ae) 
Fic. 3.8 


such an operation. It is obvious that the potential function has 
the same value at P and P’, 


i.e. V (x1, x2, %3) = V(x1', x2’, x3’) (3.15) 


The co-ordinates (x1, 2, x3) of P are related to those of 
P'(x1', x2’, x3’) by a transformation of the form (3.3). This 
transformation may be interpreted as a rotation of co-ordinate 
axes instead of as a physical rotation. Equation (3.15) expresses 
the result that the potential function is invariant when the 
transformation is applied. 

The whole Hamiltonian is, in fact, invariant under a sym- 
metry transformation. This basic result is of the utmost im- 
portance in the application of group theory to quantum 
mechanics, and is intimately related to the occurrence of 
degeneracies. 

Suppose for a given physical system a wave function cor- 
responds to an energy £, then 


H (x1, x2, X3)y(x1, X2, X3) = Evb(x1, x2, x3) (3.16) 


If the system is symmetric in some sense, then the Hamil- 
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tonian will be invariant under some transformation of co- 
ordinates, If the equation (3.16) undergoes this transformation, 
then 

H (x1', x2", x3") ap" (x1', x2", x3) = Exp’ (x1', x2", x3") (3.17) 


The new Hamiltonian # (x1’, x2’, x3’) has the same functional 
form as the original, and so the differential equation (3.17) is 
basically equation (3.16). In general, the transformed function 
ab’ (x1', x2’, x3’) is not the same as (x1, ¥2, x3) and they are two 
different wave functions denoting different states that solve the 
same equation, and belong to the energy E. 

Perhaps these ideas can be clarified by considering the par- 
ticular case of a free particle moving in one dimension. The 
energy eigenfunction e** belongs to the energy 42k2/2m. 


: h2 d2 h2k2 
1.€. om dx? etkz — (Sette (3.18) 


The Hamiltonian is invariant under the inversion through 
the origin, and so the transformed wave function e~t*2, obtained 
by changing x into —x, belongs to the same energy as e*?, This 
result is, of course, trivial, and is obvious from the form of the 
energy function. 


§ 3.3 Groups 


Group theory involves two basic ideas, groups themselves, and 
function spaces. 

A group is a set of elements {A}, which obey certain specified 
rules. In this book the elements will be taken to be transforma- 
tion operators of the type (3.4), but the concept is more general. 
The members of a group have the following properties. 


(1) It must be possible to combine the elements to form a 
new element of the group. 


€.g. A = BC (3.19) 
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If the elements are transformation operators then (3.19) 
means that if the operator C is applied first, and B second, the 
result is equivalent to applying operator A alone, where A is 
another member of the group. The order of the operators is 
important. Consider a point P on the x-axis, together with 
two lines of symmetry, m1 along the x-axis, and me bisecting 
the angle between the positive x and_y axes. 


Fic. 3.9 


If the operator mi, is applied first, followed by me, then the 
point P moves to P’. Alternatively, if mz is applied first and 
followed by mm, then the point P moves P’’. The two opera- 
tors m, and mz do not commute. 

Conversely, any two rotations about the same axis do 
commute; in this case, the order in which the operations are 
carried out does not matter. 

(2) One of the members of the group is a unit operator E 
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such that it commutes with all the other members of the 
group. 
i.e. EA = AE=A (3.20) 


If the elements are Seitz operators, the unit operator is {¢ | 0}. 

(3) Every element A in a group has an inverse A-! so that 
the combination of an operator with its inverse, results in the 
unit operator. 


AA = AIA =E (3.21) 


If A is the operator {« | £}, its inverse is {a7} | —a-1z}. 
(4) The combination of three operators must be uniquely 
defined. 


ive. A(BC) = (AB)C = ABC (3.22) 


The square group 

A simple example of a group of symmetry elements are those 
operations which leave a square invariant. It is assumed that the 
two plane sides of the square are not identical. 


Fig. 3.10 The square group 
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There are eight symmetry operations which leave the square 
invariant, they are: 


(a) The identity operator E, 

(6) The two perpendicular mirror planes mz, my. 

(c) The two perpendicular mirror planes ma, ma’. 

(d) The three rotations about the four-fold axis through o 
in the z-direction, 42, 2z, 423 corresponding to rotations of 
90°, 180° and 270° respectively. 


These eight operations may be shown to form a group, and 
Table 3.2 is the ‘multiplication table’ for the group: 


TABLE 3.2 


Multiplication table for the group of the square 4mm 


E mz: my ma Mma 42 22 423 Applied 
first 


42 42 Ma’ Ma My My 2z 423 E 

22 22 My Mz My Ma 423 E 42 
Applied 4,3 423 ma ma my mz E 42 2 
second 


In Table 3.2, the elements in the top row are applied first and 
then the operation in the first column. For example, if the 
‘operation mg is applied, and then followed by the rotation 4, 
then the result is equivalent to the reflection mz. If 4; is applied 
first and then mg, the result is my. This illustrates a point already 
mentioned, that the order of the operations is important, as they 
do not necessarily commute. 
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Point groups 

A ‘point group’ is a group of symmetry operations applied 
about a point which leave a body invariant. The group dis- 
cussed above is of this type. In infinite crystals there are thirty- 
two allowed point groups. This means the point group at any 
point in a crystal must be one of these thirty-two possibilities. 
This restriction to thirty-two follows the consideration that a 
crystal lattice must have translational symmetry. For example, 
all point groups containing a five-fold axis, or any axis greater 
than six-fold, are excluded. Seven of the groups define the 
seven crystal systems. The lattice points of all the empty 
Bravais lattices in a given crystal system have the same point 
group. The remaining twenty-five point groups arise by 
associating arrangements of atoms with each lattice point. 

If a basis of atoms is assigned to a lattice point in this way, 
the lattice point group will be reduced if the atomic arrange- 
ment has a lower symmetry than the Bravais lattice itself. 

In a crystal, each point has an associated point group. There 
is an arrangement of point groups in space. Such an arrange- 
ment is called a ‘space group’. There are two hundred and 
thirty allowed space groups in an infinite crystal. In addition to 
the point group symmetry elements, such as rotations and 
reflections, there is now the possibility of screw-axes or glide 
plane operations, which are composed of non-primitive trans- 
lations, together with a rotation or reflection respectively. 
The full set of translation operators always forms a group of 
its own. It is said to form a sub-group of the full space group. 


§ 3.4 Function spaces 


The transformation operators act on functions, so it is necessary 
to introduce the concept of function spaces, to realize the 
potentiality of group theory. A function space is defined by a set 
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of linearly independent functions. A set of functions {¢;} is said 
to be linearly independent if the relation 


aii + azpe +... t+ ann =O0 (3-23) | 


implies that the constants a; = 0 for all 7. This means none of 
these functions can be expressed in terms of the others. Any set 
of mutually orthogonal functions are linearly independent. 
Conversely, if a set of non-orthogonal functions are linearly 
independent and are also quadratically integrable, it is always 
possible to form a mutually orthogonal set from them by a 
suitable linear transformation. Suppose {¢4;} is the non-ortho- 
gonal set, then an orthogonal set {6;} can be defined as follows. 
Let 6; = ¢1 and then choose « with 62 = ¢2 + «61 so that 4 
is orthogonal to 6;. Then choose (1, B2 with 03 = $3 + 6141 
-+-B262 so that 63 is orthogonal to #1 and 42, etc. The set {64} 
may then be normalized. 

The function space defined by the orthogonal set of functions 
{6;} includes all functions of the type 


6101 + babe +... + bnOn (3.24) 


where 6,’s are any complex constants. The number of linearly 
independent functions used to define the space, gives the dimen- 
sionality. It is clear that either of the sets {¢;} and {6;} define the 
same space. This definition is in accordance with the usual idea 
of Cartesian space. The unit vectors é:, é2 and €3 are orthogonal 
vectors and are also linearly independent. Together, they define 
a three-dimensional space which includes all vectors of the type 


r = c18; + co€2 + cgé3 (3-25) 


The relationships between the function spaces and the trans- 
formation operators is very important. The function space {4} 
is said to be invariant under the set of operators {A} if all 
transformed functions of the type Ad; belong to the space. The 
Cartesian vector space, defined by the unit vectors {@} is in- 
variant under any point group. 

The space {6;} is said to be irreducible under the set of 
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operators {A} if with any function x, in the space, the set of 
transformed functions {Ax,} defines the same space where A 
runs over all the members of {A} in turn. The Cartesian space 
{&;} is irreducible under the set of three operators consisting of a 
rotation of go° about the x-axis and go° about the y-axis and 
go° about the z-axis. 

On the other hand, this Cartesian space is not irreducible 
under the set of operators consisting of the identity and rota- 
tions of +.120° about the (111) axis. To see this, consider the 
unit vector along the (111) axis. 

The space defined by the functions {i;} is said to be a sub- 
space of the function space {6;} if every function in the sub- 
space belongs to the space {6;} and if the converse is not true. 
The plane defined by the unit vectors &, & is a sub-space of 
the three-dimensional space given by &1, é2, é3. All the vectors 
in the sub-space are of the type 


r = 161 + c2€2 (3.26) 
and obviously belong to the full Cartesian space. Alternatively, 


any vector including é3, while belonging to the three-dimen- 
sional space, does not belong to the sub-space. 


Matrix representations 

If'a function space {6;} is invariant under the action of a group 
of operators {A} then the transformed functions may be ex- 
pressed in terms of the original set. 


i.e. AG: = YamiOm (3.27) 
The matrix (ams) represents the effect of the operator A on the 


basis functions {6;}. Each of the operators in the set has a 
matrix representation. If B is another operator in the set, then 


BO; = Y8miOm (3.28) 
The effect of the multiple operator AB can be seen to be 
ABO: = FD ajmB mid (3-29) 
: jom - 


if the operators are linear. 


59 


SOLID STATE PHYSICS 


The matrix corresponding to the product AB is the product 
of the separate matrices representing A and B. The matrices 
obey the same group multiplication table as the operators 
themselves, and are said to form a representation of the group. 
The order of the matrices gives the dimensionality of the 
representation. 

A different set of functions {6;’} can be used to define the same 
vector space as {0}, if the new basis functions can be expressed 
in the terms of {6;} by some linear transformation. 


6; = DSy% (3-30) 
and conversely ; 
0; = YSy18/ (3.31) 
The action of a transformation operator on 6;’, is given by 

Ad,’ = DAS: = > DSisemiOm 
= BD BSem-ameSiy9e' (3.32) 


3 


The new basis functions transform according to the repre- 


sentation 
SaS (3.33) 


This new representation is said to be equivalent to the repre- 
sentation «. A change of the basis function as described above, 
is equivalent to choosing a new set of unit vectors {é,’} in 
Cartesian space related to the old set {8} by some rotation. 

The matrices in equivalent representations corresponding to 
a given operator, need not be identical, but their traces, that is 
the sum of diagonal elements, are identical. 


Dd[SteS Je = ppp? Sim-tampSp, = 25> S7eStm—h omg 
Fi mij my it 
= > DSjmams 
m j 
= Do (3-34) 


A set of matrices representing a group of transformation 
operators with respect to a set of basis functions {8;}, is said to 
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form an irreducible representation if the function space defined 
by {64} is irreducible under the group of operators. A reducible 
representation can always be brought into irreducible form by 
an equivalence transformation of the type (3.33). 

The reducible matrices representing all the elements of the 
group are then reduced to block form, all with the same block 
structure. That is, each matrix is a direct sum of square 
matrices a1, a2, ... along the leading diagonal and the dimen- 
sions of corresponding sub-matrices are the same for every 
operator in the group. 


c.g. aj C0) oO 
fe) az ) 
a] Oo ag 
A set of sub-matrices forms one of the constituent irreducible 
representations of the group. 
The functions xe-", ye-", ze-* correspond to a reducible repre- 
sentation for the square point group. (Table 3.2 and Fig. 3.10.) 
The matrices representing the different transformation 
operators are given in Table 3.3. As an example, consider the 
effect of applying the operator mz, the mirror reflection through 
a plane perpendicular to the x-axis. 


Ma(xe", ye", Ze) == (xe, ye", Ze) [—-1 0 O 
Oo IF QO 
o oi it! 

(3-35) 


On closer inspection it will be observed that none of the eight 
group elements mixes the functions xe~*, ye-* with ze". In fact, 
the two sets (xe-', ye-") and (ze~") span separate spaces which 
are invariant under the square point group. The two separate 
irreducible representations are listed in Table 3.3 under A 
and £, 

The letters A, B are commonly used to label one-dimensional 
representations, E is used for a two-dimensional representation 
and T is used for a three-dimensional representation. 
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In this simple instance, there was no need to carry out an 
equivalence transformation to bring the matrices into the 
characteristic ‘block’ form. 

All groups have one identity representation, such as A, in 
which all the ‘matrices’ are unity. This representation is not 
faithful in that the representation possesses more relations than 
the group itself. (Compare Table 3.2 with Table 3.3.) 

There are other irreducible representations of the square 
point group apart from A and E£. Altogether there are five 
irreducible representations of the square point group, four one- 
dimensional, and one two-dimensional representations. The 
trace of a matrix corresponding to an operator is called the 
‘character’ of the operator in the representation, and is inde- 
pendent of the choice of basis functions (3.34). Table 3.4 is a 
character table for the square point group and includes all the 
irreducible representations. 

Another basis for the irreducible representation E is 
(x + iy)e-", (x — iy)e-". Clearly these functions are a simple 
linear combination of the previous basis and produce an 
equivalent representation of the square groups. The transforma- 
tion matrix is 


s=(j a and sai! a. 


1—i I 1 


The new matrices are listed in Table 3.3 and are related to the 
old by (3.33). 
e.g. for mz 


Gee dete: ) 


The group characters are the same for both sets of basis func- 
tions and are given in the last row of Table 3.4. 

The reduction of a space into sub-spaces defined by sets of 
functions forming irreducible representations for the group 
considered, can always be done, and is of the utmost importance 
in the application of group theory to quantum mechanics. The 
concept of symmetry class is useful in this context. 
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The elements of a group may be divided into distinct sets 
called ‘classes’. Given any element A of a group, then the set of 


TABLE 3.4 
Character table for square point group 


4mm E me my ma may 42 22 428 


Ay : I I I I I I I 
Ag I o—-I —I —I -—I I I I 
By I I I ~—I -—I —I I —I 
Be I —I —I!I I I —I I —I 
E 2 oO re) oO re) o —2 .0 


elements B-1AB, as B runs over all elements in the group, is 
called a class. 

The operation B-14B, is the net result obtained by first 
rotating the system to some equivalent position by B, next 
carrying out the operation A, and then reversing the initial 
operation by B-1, The members of a class must be operators of 
the same sort, such as rotations through the same angle, but 
performed about different axes which are related to each other 
by members of the group. 

The characters of all the elements of a class in a given re- 
presentation are obviously identical (3.34). In a group 
character table it is necessary only to specify the classes. In 
Table 3.4 there are five classes 


E 
Maz, My 
Md, Mq’ 


42, 423 
22 


The number of classes in a group is equal to the number of 
irreducible representations. For the square point group, there 
are five classes and five irreducible representations. If n; is the 
dimensionality of an irreducible representation, and N is the 
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number of elements in the group, then the following relation 
can be derived: 

Dn? =N (3.36) 
t 


The sum is over all the possible irreducible representations. 
If the number of classes in a group is known, and if the number 
of elements in a group is also known, then using (3.36) the 
dimensionality of the irreducible representations can be found. 

For the square point group, the number of classes, and hence 
the number of irreducible representations, is five. 

There are eight elements in the point group, and from (3.36), 

22-12-72 + 724 72 =8 (3-37) 

There are four one-dimensional and one two-dimensionalirre- 
ducible representations for the group as shown in the Table 3.4. 

The relation (3.36) has important consequences for Abelian 
groups. An Abelian group is one in which all the elements 
commute with one another. The elements which compose a 
four-fold rotation axes form such a group. Each element in an 
Abelian group forms a class in itself. This means the number of 
irreducible representations is equal to the number N of ele- 
ments in the group. From the equation (3.36), it is obvious that 
all these irreducible representations must be one-dimensional. 
This result is referred to again in Chapter 5 with regard to the 
Bloch theorem. 

Sufficient background knowledge will now have been gained 
by the reader concerning the theory of groups and their repre- 
sentations. The next step is to relate the group theory to 
quantum mechanics. 


§ 3.5 Symmetry and quantum mechanics 


Section 3.2 has shown that there is a close link between the 
occurrence of degeneracies and symmetry operators. If y1(r) is a 
solution of some Hamiltonian #, then 


Hn(r) = Eyi(r) (3-38) 
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If the Hamiltonian is invariant under the transformation 


operator A then 
AH =H# 


Operating with A on the left side of equation (3.38), # and A 


commute. 
Addy =#Ads 


The transformed wave function 


Ady = de 


is also a solution of the Hamiltonian, corresponding to the same 
energy. 

If #1 and ye are linearly independent, they represent de- 
generate eigenfunctions. 

Corresponding to each energy level there is a set of degenerate 
eigenfunctions, which span a space which is invariant under the 
action of any group of transformation operators which com- 
mute with the Hamiltonian. The eigenfunctions produce some 
representation of the group. The basic theorem relating group 
theory to quantum mechanics postulates that the set of all 
degenerate wave functions corresponding to a given energy 
level, forms an irreducible representation of the group of all the 
operators which commute with the Hamiltonian. 

In general, this statement is true, but sometimes it happens 
that non-equivalent irreducible representations correspond to 
the same energy. This occurrence is termed an ‘accidental 
degeneracy’. 

For a multi-electron wave function, account must be taken 
of the permutation operators which commute with the multi- 
electron Hamiltonian. This additional complication does not 
arise in the one-electron approximation. 

It is a well-known fact that when a physical system is per- 
turbed, degeneracies are often destroyed. This is easily under- 
stood if it is observed that in general, the perturbed Hamil- 
tonian possesses a lower symmetry than that which is unper- 
turbed. The erstwhile degenerate wave functions will now form 
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a reducible representation for the new smaller symmetry group, 
and will split into two or more irreducible representations. 
Hence a degenerate level may be expected to divide into two or 
more levels of diminished degeneracy. 


§ 3.6 Some applications of group theory 


Axial rotation group 


The group composed of all the rotations about some axis is 
called the axial rotation group. If all angles of rotation are 
allowed, the group has an infinite number of elements. The 
group is Abelian, since a rotation ¢1, followed by ¢2, is exactly 
equivalent to ¢2 followed by ¢1. Each rotation forms a class 
in itself, and so all the irreducible representations are one- 
dimensional. The representation must satisfy 


T(pi) (ge) = T(pi + ¢2) 
The solutions of this equation are 
Pom (g) = ems 
The index m is restricted to the values 
m=o0, +1, +2,... 
by the requirement that 
T@) (27) = P)(0) = 1 


The basis functions for these irreducible representations are 
of the form 


Ym(r, 8,6) =f(r, Oem  m =o, +1, +2,... 


The eigenfunctions for a Hamiltonian whose symmetry 
group is the group of all rotations about an axis, are of this 
type. This expresses the result that angular-momentum is con- 
served along an axis of symmetry. 
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Three protons positioned at the vertices of an equilateral triangle 


The Hamiltonian for a single electron moving in the field of 
three protons arranged in the form of an equilateral triangle is 


ree sy e2 e e2 I 
FN ge lr—nlt[r—nl *[r—vrs] 4meo 


where ri, re and rs denote the positions of the protons. 


IN 


The origin is at the centre of the triangle 


Fic. 3.11 The group of the equilateral triangle 


This system is not of any great importance in itself, but the 
symmetry is related to that of an ion surrounded by three 
molecules of water in the hydrated crystal of a salt. (The ozone 
molecule has almost this symmetry.) 

There are twelve symmetry operations which leave the 
system of three protons invariant and they form the group 6m. 
They are: 


(a) the identity operation E; 

(6) the three mirror planes of the type my; 

(c) the three rotations of the type 2, through 180° about 
the vectors rn, re and rg respectively ; 

(d) the mirror plane mz perpendicular to the z-axis and 
including all three protons; 


68 


THE IMPORTANCE OF CRYSTAL SYMMETRY 


(e) the two rotations 32, 3.8 corresponding to a positive 
rotation of 120° and 240° about the z-axis respectively; 

(f) the two operations 6, involving a rotation through 60° 
and 300° about the z-axis followed by an inversion through 
the origin, 


TABLE 3.5 
Character table for the point group 6m2 


Number of elements 
in Class I 3 3 I 2 2 


Particular element 


in class E My 22 Mz 32 6z 

Representation 
At’ 4 I I I I I 
Ad’ I —I —I I 1 I 
A,” I —J I —I I —I 
Ag” I I —I —I I ~—I 
E 2 re) re) 2 —I ~—I!1 
Ee” 2 re) o —2 —I I 


The point group 6m2 is composed of six classes, so there are 
six irreducible representations of the group. This observation, 
together with the equation 3.36 determines the dimensions of 
the representations. 


12 +4 12 4+. 12 + 1% 4 92 4+ 92 = 12 


There are four one-dimensional and two two-dimensional irre- 
ducible representations of the group. 

The eigenfunctions of the Hamiltonian belonging to the same 
energy must form a basis for one of the representations. It is to 
be expected that the ground state has the maximum symmetry, 
and so the state is non-degenerate and the eigenfunction is 
a basis of the representation Ai’. The function z? has this 
symmetry. 
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The hydrogen atom in an electric field with symmetry 6m2 


In this section is considered the effect of an electric field with 
symmetry 6m2 on the degenerate (2) level of the hydrogen 
atom. 

Before the electric field is applied, the single hydrogen 
electron moves in the spherically symmetrical field due to the 
‘fixed’ proton. 

The electron Hamiltonian is invariant under all rotations 
about any axis through the proton taken at the origin; 

i.e. the Hamiltonian isinvariant under the full rotation group. 

The (2p) states 


(2pz) =f (r)x, (2hy) =S (ry, (2f2) =f (r)z 


with f(r) a function of 7 only, transform according to an irre- 
ducible representation of this group and belong to the same 
energy. The (2p) level is triply degenerate if spin degeneracy is 
ignored, 

It is quite clear by inspection of Table 3.5 that when the 
electric field is applied and the electron Hamiltonian has the 
reduced symmetry 6mg then at most a double degenerate level 
remains. The three (2) wave functions transform according to 
a reducible representation of the group 6m2 which i is given in 
Table 3.6. 


TABLE 3.6 


Pry Py, pz 3 +1 aq cha? Zor exe 


This representation is composed of the two irreducible Eepre: 
sentations Ae’ and E’. 

The previously triply degenerate (2p) level has split up into a 
single level and a doubly degenerate level under the action of 
the electric field. 

A similar type of argument is used in ‘crystal field ery 
problems. 
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Problems for Chapter 3 


1. Let a1, az and ag be a particular set of basis vectors for 
some three-dimensional Bravais lattice. Define three vectors 
ui, uz and us by 


Wy = la, + maz + mas 
2 = lear + meae + neas 
ug = Ilsa, +- m3a2 + n3a3 


where /;, mi, . . . ms, ng are integers. The condition that these 
vectors form a new basis is . 


@ = 11in + sive + tins 


with 71, 51, ¢; integers and similarly for @z and as. Show that this 
condition is satisfied if the determinant 


bh le ds = +I 
m1 mao m3 
my Ne 7g 


Verify that the primitive cell volume remains constant under 
the change of basis. 

2. (a) Show that there is a two-fold axis going through every 
lattice point of a two-dimensional Bravais net. 

(b) Show that if a Bravais net has an axis of n-fold symmetry 
about some arbitrary point it has at least an axis of n-fold 
symmetry about any lattice point. (Hint — Consider the funda- 
mental rotation about the arbitrary point P which takes the 
lattice points A and B to B and C respectively. Then apply the 
translation which returns A to its original position and note the 
final position of B.) 

3. Prove that a Bravais net can only have an axis of two-, 
four- or six-fold symmetry, about a lattice point. Show also that 
a three-fold axis may exist at an arbitrary point. (Hint — Use 
the results of Problem 2 and consider the fundamental 
rotation of the shortest possible basis vector a; about a lattice 
point.) 
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4. Prove that the elements of the rotation transformation 
matrix (see equations (3.5) and (3.7)) 


11 G12 013 

M21 M22 23 

O31 O32 33 
satisfy the relations 


3 * 
> aj? =1 for alli 
i=1 


3 
2 mite =o jx#m 
= 


Prove also that the determinant of the transformation matrix Is 
ET. 

5. Describe the four symmetry operations under which a 
rectangle is invariant. It may be assumed that the two faces of 
the rectangle are not identical. 

Draw up a group multiplication table. (This is the group 
2mm.) 

6. Describe the six symmetry operations under which an 
equilateral triangle is invariant. It may be assumed that the 
two faces of the triangle are not identical. Draw up a group 
multiplication table. (This is the group 3m and is a sub-group 
of the group 6m2 studied in the text.) 

7. The permutations of n objects form a group. This is called 
the ‘symmetric group’. There are six permutations of the 
objects 1, 2, 3. A permutation that replaces 1 by 2, 2 by 1, and 
leaves 3 invariant may be denoted by the symbol 


& 2 ) 

2 1 3 

The product of two such permutations is another of the six 
possible permutations. 


<6 (EMT 9-029 
2 3 1/\2 1 8 32 1 
Draw up a group multiplication table for these permutations. 
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Show that this table can be made to have the same form as that 
of the group 3m of question six by pairing off each of the ele- 
ments of the permutation group with a suitable element of the 
group 3m. The two groups are said to be ‘isomorphic’. (Since 
all electrons are identical, a many-electron Hamiltonian is 
invariant under the group of permutations of the electron 
co-ordinates.) 

8. The Hermitian adjoint At of a real operator A may be 
defined by 


[Add dr = [tary ar 


where ¢, % are arbitrary functions and the integration is over 
all space. If A is a linear operator which possesses an inverse 
A-1, it is said to be a unitary operator if At = A-}, 

A change of co-ordinate axes can not affect the values of the 


integra] over all space J ¢*y dr. That is, the integral is in- 


variant under a symmetry transformati®n of the type {a | Rn}. 
Explain how this implies that {« | Rn} is a unitary operator. 

g. (a) Show that the function space (x?, y2, x2 + 2) has the 
dimension two and not three. 

(6) Similarly show that the function space (xy, 2, 2K, 
x(y — z)) has dimension three and not four. 

10. Describe the 48 symmetry operations which leave a cube 
invariant. These form the full cubic point group mgm. Divide 
the elements into the ten classes. Find the dimensions of all the 
irreducible representations. _ 

11. List the 24 symmetry operations which leave a tetra- 
hedron invariant. This is the point group 43m. Divide the 
elements into the five classes and give the dimensions of all the 
irreducible representations. 

(The carbon atom site in the methane molecule CHa and in 
the diamond lattice has the symmetry of this point group.) 

12. In Chapter 3 it was shown that there are eight elements 
in the square point group 4mm. The group character table is 
given below. 
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No. of elements 1 I 2 2 2 
in class 


Typical element 


in class E 22 42 Mz ma 

Irreducible 

representation 
Aj I I I I I 
Ag I I I —I —I 
By, I I —I I —!I 
Be I I -—-I —I I 
E 2 —2 a) re) oO 


The character of the group element A in the «-irreducible re- 
presentation is denoted by (A). Verify the ‘orthogonality 
relations’. 


ZxO(A )x(A)=0 a #8 
SrA x(A) = 


The summation is over all the N elements in the group. 
13. The identity E and the two-fold rotation 2, form the group 
2. The character table is 


E 22 
A 1 1 
B I a 


A projection operator may be defined as 
O™ = Yy(A)A 
A 


where the sum is over all the operators in the group and 
x (A) is the character of A in the a-irreducible representation. 
The effect of this operator on an arbitrary function f(x) is to 
‘project’ out that part which transforms according to the 
a-representation. 
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Apply the projection operator method to the function f(x) 
and obtain functions which transform according to the repre- 
sentations A and B of the group 2. 

14. The functions (xy, xz, yz) define a three-dimensional 
function space. By the projection operator technique pick out 
those functions in this space which transform according to the 
irreducible representations Bz and E of the square point group 
4mm. (See Problem 12.) 

15. The atomic orbitals (2s), (262), (2py) and (2) form the 
basis for a reducible representation of the tetrahedral point 
group 43m. (See Problem 11.) Find five matrices of this repre- 
sentation each corresponding to a group element from one of 
the five classes. Show that this representation is composed of a 
single one-dimensional and one three-dimensional irreducible 
representations. (Hint — Choose the co-ordinate axes so that a 
3-fold axis lies along the (111) direction). (The (25), (2p2), 
(2hy), (262) carbon atomic functions are available for bond 
formation in the CH, molecule.) 

16. Consider an electron moving in the field of four protons 
positioned at the vertices of a tetrahedron. An eigenfunction of 
the Hamiltonian must be associated with one of the five irre- 
ducible representations of the tetrahedral point group 43m. 
(See Problems 11 and 15.) A uniform electric field is applied 
along one of the three-fold tetrahedron axes from the centre to a 
vertex. The reduced symmetry group is the point group 3m 
composed of six elements. The group character table is given 
below. 


No. of elements 


in class 2 3 
Element in class E A B 
Ay I I I 
Az I I —I 
As 2 —I Oo 


Identify the elements A and B. 
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The three atomic functions fz, py, Pz form a basis for an irre- 
ducible representation of the tetrahedral group. (Problem 15.) 
Show that this representation is composed of two irreducible 
representations of the new symmetry group 37. 

17 (a) For a real Hamiltonian the process of taking the com- 
plex conjugate commutes with #. 


i.e. (4)* = H*$* = #4* 


This is an additional operation which may produce degeneracies 
between eigenfunctions that belong to different representations 
that are irreducible under the symmetry group of the Hamil- 
tonian. This operation of complex conjugation may be denoted 
by XK. Show that K is not a linear operator but has a similar 
property, and also that 


| K[ pap de | = |[ pry dr | 


K is said to be anti-linear and anti-unitary. What is the inverse 
of K? 3 

(In the general case of a many particle Hamiltonian in- 
volving spin, the corresponding operator is the time-reversal 
operator T.) 

' (b) Suppose a uniform electric field is applied to a hydrogen 
atom. The eigenfunctions of the new Hamiltonian transform 
according to the representations of the axial rotation group ©. 
In Chapter 3 these irreducible representations were shown to be 
one-dimensional. Explain why the new energy levels need not 
all be non-degenerate. (This is not a true example of extra 
degeneracy arising from K as it is also a consequence of the 
group com, composed of all reflection planes through the axis 
of symmetry.) 


Recommended further reading 


Heine, V., Group Theory in Quantum Mechanics, Pergamon Press, 
London (1960). 
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Johnston, D. F., ‘Group Theory in Solid State Physics’. Reports 
on Progress in Physics, XXIII, 66, 1960. 

Tinkham, M., Group Theory and Quantum Mechanics, McGraw- 
Hill, New York (1964). 
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The Theory of Lattice Vibrations 


§ 4.1 Longitudinal waves in continuous media 


The atoms in a real crystal are never actually at rest. In this 
chapter the elastic vibrations of solids are discussed. For wave- 
lengths that are large compared with the interatomic spacing, 
the crystal may be regarded as a continuous distribution of 
matter. 

A longitudinal wave is one in which the particle displacement 
€ is in the direction of wave propagation. The equation of 
motion of such waves in a continuous medium (Coulson 
“Waves’) is for propagation along the x-axis 


02g x 02 2 A 

Ox2 2 OR where € p (4.1) 

A is a constant of elasticity and p is the density. The general 
solution of this one-dimensional wave equation 1s 


S(% — ct) + g(x + et) (4.2) 


where f, g are arbitrary twice differentiable functions of x and t. 
The constant ¢ is the velocity of all the waves. Particular solu- 
tions are the harmonic waves exp i(wt + kx) with ¢ = w/k. 
In addition to the longitudinal waves there are two types of 
transverse waves which are propagated with a velocity which is 
less than c. Of course the distinction into longitudinal and trans- 
verse waves can only strictly be made in an isotropic medium. 
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§ 4.2 Classical theory of lattice waves 


When the wavelength of the elastic waves is of the order of the 

interatomic spacing it is clear that the microscopic structure of | 
the crystal must be important. In a crystal the waves are trans- 

mitted by the vibrations of the atoms. The binding force be- 

tween atoms can arise in several ways. For ionic crystals the 

force of attraction is the Coulomb field, for metals it is due to the 

free electrons and for valence crystals it is the covalent bond. 

When the atoms approach each other closely, short-range re- 

pulsive forces become dominant and prevent the collapse of the 

lattice. 


Wave motion along a line of identical atoms 


To a first approximation the force between any two atoms is 
proportional to their relative displacement from their equili- 
brium separation. Several important effects such as thermal 
expansion and heat resistivity are due to deviations from this. 
Nevertheless it is instructive to study the transmission of longi- 
tudinal waves down a line of atoms within this approximation. 
An additional simplification is to consider nearest neighbour 
interactions only. 

Consider a line of identical atoms, each of mass m, which are 
a distance a apart in their equilibrium positions. x, denotes the 
displacement of the nth atom from its equilibrium position. The 
resultant force (positive when acting to right) on the nth atom is 


Fy = &(%n+1 —~ Xn) — &(%n — Xn-1) (4.3) 
where « is an elastic constant. The equation of motion is 


os d xn 
Min = &(%n41 + Xn-1 — 24m) in = “diz (4-4) 


There exist travelling wave solutions in which the atoms 
vibrate with the same frequency and for these 


Xn = Agiut+kna) (4.5) 


k is called the wave ‘vector’. A is a complex constant, but it is 
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understood that the real part of the r.h.s. is intended. Substitu- 
tion of (4.5) into (4.4) gives 


4 
Ww == +(#) sin ka 


= (4.6) 
Equation (4.6) is the dispersion relation between the angular 
frequency w and the wave vector k. As w is a periodic function 
in k with period 27/a then for a given frequency there are an 
infinite number of k-values such that 


m 


kim = k! + m= m integer (4.7) 


with | k’ | < z/a. From equation (4.5) it is quite clear that each 
of the k-values related by (4.7) represent the same atomic 
displacements and are therefore equivalent. This is explained 
when it is remembered that the wave description is here em- 
ployed only to describe the motion of a discrete set of points 
rather than a continuous line of matter. 


Fic. 4.3 Two equivalent lattice waves 


Figure 4.3 illustrates the two equivalent cosine waves with 
k’ = w/2a and k; = 57/2a. Both waves give the same atomic 
displacements. It is common practice to restrict k to lie within 
the range —z/a <k < a/a and this region is called the Bril- 
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louin zone. (Brillouin zones will be discussed in detail in 
Chapter 5.) Also, for each w, there are two values of k inside the 
Brillouin zone, which are equal in magnitude but opposite in 
sign. These correspond to waves moving in opposite directions. 
This degeneracy occurs also in two- and three-dimensional 
lattices. 


At the ‘maximum’ frequency wm, the solution is a standing 
wave. 


Xn = Aetot cos na (4.8) 


There are in fact solutions corresponding to frequencies greater 
than wm but then k is complex and the waves are exponentially 


2 
attenuated. For w > wm then (sin =) > 1. The solutions are 
k =< + if 6 real 
2 2 
with (cosh =) = (4.9) 


B is called the attenuation coefficient and increases rapidly with 
dw. 


Fic. 4.4 Attenuation coefficient 
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The phase velocity of the waves is a function of k and 


__w _ 7, sin (ka/2) 
iy oe Voge (4.10) 
a\t 
where Vp = (=) a. 
m 
Y 
M, 
| 
| 
| 
| 
| 
| 
| 
@) 
us k 
qg 


Fic. 4.5 Phase velocity 


In the limit of long wavelengths (small k) the atomic de- 
scription must be identical with the continuum treatment. 
From (4.10) the limiting value of the velocity as k —-» 0 is Vo. 
This should be compared with equation (4.1). 


c? = Np 


The modulus A = «a and the average linear density p = m/a. 


c= — — Py 


as required. 

In a crystal a reasonable value for the k giving the maximum 
frequency is k = w/a ~ 108 cm-—!, The wave velocity ~10° cm/ 
sec and so the corresponding maximum frequency is ~10}8 rad/ 
sec and lies in the infra-red region. 
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Finite linear lattice 
Consider a line of N + 1 atoms, the first and last atoms held 
fixed. For a frequency w, the general solution of (4.4) is 

Xn = etut( Agtkna + Be~tkna) (4.1 1) 


A, B are arbitrary constants. The appropriate boundary condi- 
tions are 


Xn =O n= 0 (a) 
tn =O n=N (8) 
The first condition is satisfied if 4 = —B. The second condition 
requires 
sin (kNa) =o 
lar 
k= Wa integer (4.12) 


The solutions are standing waves and 
; : L 
Xn = A’etot sin (n 7) (4.13) 


It is clear from (4.13) that the negative values of | give the same 
solutions as the positive values and in fact there are N — 1 
allowed standing waves corresponding to the N — 1 particles 
which are free to vibrate. 


Lif 2,2 iy (Nh 4) (4.14) 


Periodic boundary conditions 


It is often more convenient to obtain travelling wave solutions. 
The infinite crystal is then assumed to be composed of identical 
microcrystals each containing N lattice points. The boundary 
conditions are | ; 


xn = Xn+n for all n (4.15) 
A travelling wave of the type (4.5) satisfies this condition if 


In 


k=277 


i integer (4.16) 
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The allowed values of k are spaced twice as far apart as those 
of (4.12). If the number of lattice points in the microcrystal is 
odd 
N— 
b= 41, 42,..,, 4a) (4.17 (2) 
and if WN is even 


f= +1, +2,..., +2 (4.17 (2) 


In this latter case, the states with / = 45 are equivalent (4.7). 


The total number of vibrations is N — 1 as before as now the 
positive and negative values of / correspond to independent 
states. Note that now in addition there is a non-zero solution 
with J = o (4.16) and that this corresponds to a translation of 
the crystal without change in potential energy and with zero 
frequency. , 


Frequency distribution 


A frequency distribution g(w) may be defined so that g(w)dw is 
the number of allowed frequencies between w and w + dw. 
The number of states in the element dk is Nadk/m (4.12). (With 
periodic boundary conditions the number of states between 


|k| and | & + dk | is (22) dk = we as before.) 


‘Qn 
Na dk 


For a line of identical atoms (4.6) 


k= 2 sin (=) 


Wm, 
ee 
dw ~~ a(wm? — w?)t 
2N I 
gw) =—. eo (4.18) 
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This distribution function has a singularity when .w = wm and 
the frequency modes are crowded near the cut-off frequency. 


g (w) 


@) Wm w 


‘Fic. 4.6 Frequency distribution for the monatomic linear 
lattice 


Interactions beyond nearest neighbours 


The exclusion of interactions beyond the nearest neighbours is 
only a convenient approximation and a more accurate de- 
scription is obtained if interactions over larger distances are 
considered. When the Mth neighbours are included the obvious 
generalization of (4.4) is 


M 
min = 2» t(tnts + Xn-s — 2Xn) (4.19) 


The corresponding dispersion relation is 
= . 
w?m = 2> as(1 — cos (ska)) (4.20) 
s=l 
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Once more w? is periodic in & with period 27/a and inside the 
Brillouin zone w? is an even function and there is two-fold de- 
generacy (except at the origin). As cos (ska) may be expanded 
as a polynomial of degree s in cos ka the frequency may be 
expressed as a polynomial of degree M in cos ka. For a given 
frequency there are M solutions for cos ka and consequently M4 
values for k in the range 0 to w/a. It is possible that some of these 
values are complex and correspond to damped solutions. 


pc ll O aK 
a 


Fic. 4.7. Dispersion curve with many neighbour interactions 


The maximum frequency does not necessarily occur at the 
points -+-7/a although the dispersion curve must have a 
stationary value at these points. 


Transverse waves | 
If the atoms can vibrate in a plane perpendicular to the line | 
then transverse modes are possible. For small displacements it | 
can readily be shown that the equations of motion are very 
similar to those for longitudinal motion (4.4). The lateral dis- 
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placement y, of the nth atom satisfies (nearest neighbour 
interaction only) 


min = a'( ynti + Yn-1 — 2yn) 
There is of course a new force constant a’, All the analysis of the 


previous sections will apply and for a linear lattice with 
cylindrical symmetry there will be two dispersion curves. 


Fic. 4.8 Transverse and longitudinal waves 


§ 4.3 A linear lattice with two types of atoms 


Consider a line of atoms equally spaced at a distance a apart. 
Alternate atoms have masses M and m (M > m). The primitive 
cell has a length 2@ and contains two atoms. 


2n-| 2n 2nt+| 2nt+e2 
° e ° e re) @ 
M m M mM M m 
+> 
a 


Fic. 4.9 Linear lattice with two types of atom 
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The masses M are located at the even sites 2n and the masses m 
at the odd sites an + 1. When only nearest neighbour inter- 
actions are considered the equations of motion for the two 
types of atom are (longitudinal modes) 


Mitan == &(Xan+1 + xan-1 — 2%n) 
Miton+1 = a(Xen+2 + Xen — 2Xen+1) (4.21) 
There exist travelling wave solutions of the type 
Xan = Azetlwttk2na) xoq 4) == Agetwttk2nt+ia) (4,99) 
The equations of motion are then 


Ai(w?M — 2a) + Az2a cos ka = 0 
Ai2« cos ka + A2(w2m — 2a) = 0 (4.23) 


These two simultaneous equations have a non-trivial solution 
only if the determinant of coefficients is zero. This is a quadratic 
expression in w? with 


Oar : 2k 
OR gga a 7 a =) als al (a7 oh =) = 4 ae aoe] (4. 24) 


The dispersion relation (4.24) is periodic in k with period r/a 
and the Brillouin zone is defined by —7/2a < k < m/2a. These 
two sets of solutions are known as the optical and acoustical 
branches. For long waves (small k), the two roots are 


wo? = 2a a + =) optical 


2% 


eee : 
Wm: ka acoustical (4.25) 


w2 = 
For k = +72/2a, the two roots are 
optical 


ok 
m 
wi? = a acoustical (4.26) 
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Acoustical 
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Fic. 4.10 Optical and acoustical branches 
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Fic. 4.11 Atomic motion for small k 
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From (4.23) and (4.25) it can be shown that in the limit of 
small k 


Ai oe | acoustical branch 
Ag 

Ay a m : 

yr 3 optical branch 


If the particles are ions with opposite charges then this latter 
type of motion may be excited by an electromagnetic wave. 
This maximum frequency represents the fundamental infra-red 
absorption frequency. 

There is a frequency gap in the dispersion diagram (Fig. 4.10) 
in which k is complex and the waves are damped. 


§ 4.4 Normal co-ordinates 


Consider a crystal (or molecule) composed of N atoms. 3N 
position co-ordinates q; are necessary to specify the positions of 
all the particles. If the interatomic forces are conservative a 
potential function V may be defined which depends upon the 
position co-ordinates. This potential has a minimum in the 
stable equilibrium configuration specified by 9g; = a and all 
the first order derivatives are zero. If the system undergoes a 
small displacement from equilibrium so that g; = «; + ¢; then 
by Taylor’s theorem to second order in ¢; 


Vs al (Ga: toe ft lw zee) V |. (4.27) 


Here the zero of potential is the minimum. In this Harmonic 
approximation the potential energy is a quadratic function of 
the displacements. This corresponds to the Hooke law approxi- 
mation in Sections 4.2 and 4.3 (see Problem 1). 


3N 
V= + 2 fubebs 


: 02V 
with Su = Sr = (scan). = constant (4.28) 
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The position co-ordinates may be chosen relative to a Car- 
tesian set of axes and the kinetic energy is 


3N dé 
Ls a2, moi? c= di (4.29) 


The normal co-ordinates Q; are defined as the linear com- 
binations of the co-ordinates £; which bring the kinetic and 
potential energies simultaneously into the diagonal forms 


T 2s +> Oe V = FS otQe (4-30) 


These co-ordinates can always be found and the process is 
briefly explained below. The kinetic energy (4.29) is easily 
dealt with by introducing ‘mass-weighted’ co-ordinates 
™m = Viti and 


3N 
T=43%2°  V= LaF (4-31) 


where Fi = fylV mi. mj 


Any further orthogonal transformation will preserve the kinetic 
energy as a simple sum of squares. For example, define a new 
set of co-ordinates by an orthogonal transformation Ay. 


3N 
m= 2 Au) (4.32) 
* 


In matrix formalism.» = A@ where 7, 6 are column matrices. 
As A is orthogonal A-! = A’ (A’ is the transpose matrix) and 
SO 


9 = Ay = A’y (4-33) 
The kinetic energy term in (4.31) is 
T = 49'9 = 3(A6)’(A6) = 46’A-146 = 466 
and the sum of squares is preserved by the orthogonal trans- 
formation. When the transformation is applied to the potential 


energy 
V = 4n/Fy = }(A0)'F(A0) = 40'(A-1FA)O 
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It can be shown that the orthogonal transformation can always 
be chosen so that (A-1FA) is a diagonal matrix with 


(A-1FA)ig = wi?Sz (4-34) 


(A complete discussion is given in Synge and Griffith’s Principles 
of Mechanics.) 
The co-ordinates @; are then the normal co-ordinates of (4.30). 


The momentum P; conjugate to Q; is defined by P; = 5, = 0: 


and the Hamiltonian is 


3N aN 
H= t2P et $d oO? (4.35) 
oH 


From Hamilton’s equation P; = 


— 30; 
GO, +02Q:=0 i=1,2,...3N (4.36) 


For stable equilibrium V is positive definite and w;?2 > 0 for all 
z. The solutions of (4.36) are 


Qs = Q:° sin (wit + «) (4.37) 


As the particle displacements are real the normal co-ordinates 
are real and vary harmonically in time with period 27/w;. There 
are 3N such co-ordinates and frequencies corresponding to the 
3N degrees of freedom of the system. 

A general solution in terms of the co-ordinates 7; may be 
expressed as a linear combination of normal co-ordinates (4.32). 
The amplitudes and the phases of the normal co-ordinates are 
6N constants chosen to fit the initial conditions of the system. 
If these are such that the 3N — 1 normal co-ordinates Qu, 
Qi, ~~ +5 Qr-1, Qr41, +++; Qsw are all permanently zero and 
Q, # o then the system is vibrating in a normal mode. The 
motion of each particle is periodic with the same frequency wy 
and with the same phase. This is the type of solution obtained 
for the linear lattice in section 4.2 (equation (4.5)). . 

The allowed frequencies may be obtained from (4.34) or 
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from the equations below which are derived from (4.31) in the 
same way that (4.36) is obtained from (4.30). 


ie + Shim =0 8 =t=1,..53N (4.38) 
3 


When the system is in the rth mode = AirQy and (4.38) 
becomes 


2 [Fu — wSy]Ajyp =O ti=1,...,3N (4.39) 


For a non-trivial solution the determinant of coefficients is 
zero. 
| Fy — wry] = 0 (4.40) 


This is a polynomial of degree 3.N in w? and gives the 3.N values 
for w which are characteristic of the particle system. If two or 
more of the roots are coincident the system is said to be de- 
generate. For a given wy the equations (4.39) determine the 
relative values of the coefficients Aj, j = 1, 2,..., 3N. When 
this set is normalized it forms a column of the orthogonal trans- 
formation matrix. 

It has been stated that V is positive definite and w;? > 0 for 
all frequencies. This is not strictly true and in fact six frequencies 
are zero, (Five for a linear molecule.) These correspond to modes 
which represent a translation or rotation of the crystal. Clearly 
if each atom is given the same translation there is no change in 
potential energy and hence no restoring force. The system is in 
neutral equilibrium and V is semi-definite with respect to the 
six zero-frequency co-ordinates. 

If the position of the centroid and orientation of the system 
are specified then there are 3N — 6 ‘internal’ co-ordinates all 
with real, non-zero frequencies. For the longitudinal vibrations 
of a linear lattice (with periodic boundary conditions) there is 
only one zero frequency. However, if the ends of the crystal are 
fixed this translational mode is not allowed (4.14). 
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§ 4.5 Group theory and normal modes 


A normal mode can be described by the 3N co-ordinates 
m° = AirQr° giving the maximum displacements of the atoms 
from their equilibrium positions. The molecular symmetry may 
demand that several different modes possess the same frequency. 
For example consider an equilateral triangular molecule. 
Figure 4.12 (a) illustrates the atomic displacements for one of 
the normal modes. The distortion in 4.12 (b) is related by a 
mirror plane and so must represent a mode with the same 
frequency. It is in fact the same mode. 


Mm, | 


| 
3 


(a) (b) 
Fic. 4.12 
For a non-degenerate frequency w; the normal co-ordinate must 


transform according to a real one-dimensional representation of 
the symmetry group of the molecule. 


i.e. RQ = £9: 

where R is any member of the group. For an n-fold degenerate 
frequency, the n normal co-ordinates Qi,...,; Qn (say) must 
transform according to a real n-fold representation of the group. 


RQ: =, CQ: (4.41) 
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This representation will be irreducible apart from accidental 
degeneracies. 

As an illustration consider the triangular ‘molecule’ of Fig. 
4.12. If the two faces are not identical the symmetry group is 
3m containing six elements. (It is suggested that the reader 
extend this treatment to the case of the ozone molecule with 
symmetry 6m2. See Section 3.6.) Table 4.1 is the character 
table. 

TABLE 4.1 


Character table for group 3m 


3m E 32 Mea 
Aj I 1 I 
Ag I I —I 
E 2 —I re) 
Xtot 9 oO I 


Nine co-ordinates are necessary to specify the positions of all 
the atoms and so the representation of the group given by 


L, 


x 


Translational modes Rotational mode 
E A> 


Fic. 4.13 Some zero-frequency modes 


the normal co-ordinates is nine-dimensional. This ‘total’ repre- 
sentation is reducible. The characters of a representation are 
invariant under a change of basis functions (3.34) and so the 
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atomic displacements may conveniently be used to find the 
total representation. These characters are listed under ytot in 
Table 4.1. Note that for a given symmetry operation only those 
atoms which remain fixed can contribute to the character. 
By inspection the total representation is 


2A; + Az + 3E 


The translations account for representations A; -+ E and the 
rotations for Az ++ E. These are the six zero frequency modes. 
The remaining three vibrational modes form a basis for the 


A E 
Totally symmetric 
mode 


Fic. 4.14 Vibrational modes 


representation A; + £. One of the vibrational modes is non- 
degenerate while the other two form a degenerate pair (Fig. 


4.14). 


§ 4.6 Normal modes and the linear lattice 


The vibrations of the linear lattice have been discussed in some 
detail in Sections 4.2 and 4.3, but it is instructive to relate the 
results obtained with the theory of normal modes, 

Consider an infinite line of identical atoms distance a apart. 
Using periodic boundary conditions, the lattice is assumed to be 
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composed of microcrystals containing N atoms. The N normal 
frequencies are given by (4.6) (or more generally by (4.20)) with 


= 


RS PS a N odd 


Na 
1 = o gives the zero-frequency translational mode. Each of the 
normal co-ordinates must be real so that the atomic displace- 
ments are real when the system is vibrating in a normal mode. 
However, it is customary to introduce complex co-ordinates Ry 
defined by 


Nx1 99 
R= TN2 SS etknay, or Ry =— VN Dette (4.42) 
Conversely, from the result 
N—1 
2 on ; 
2, Gen Nan: (4-43) 
aS 
N-1 
I a 27 
*n = TH 2 Oras (4-44) 


(Compare (4.32) and (4.33).) 
On substitution of (4.44) into equation (4.4) (or (4.19)) it is 
easily verified that 

Ry + wR, = 0 (4-45) 


Compare this with equation (4.36). The general solution for R; 
is 
Ri = Aj-e-tt + Aytelot w > 0 (4.46) 

A necessary condition to ensure the reality of the atomic dis- 
placements is R; = R*_; so that the complex co-ordinates are 
not completely independent. The correct number of inde- 
pendent normal co-ordinates is obtained by taking the real and 
imaginary parts. The kinetic and potential energies of the linear 
lattice are 

T =— > tn? = 2S (ee — ene)? (4.47) 


2 n=0 2ne a ( 
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V is correct for first neighbour interactions only. It is left as an 
exercise for the reader to confirm that 


oP RRs (4.48) 


These equations remain true when further neighbour inter- 
actions are included. (Compare (4.30). The present co- 
ordinates are not mass-weighted.) 

The linear lattice is invariant under the group of N transla- 
tions {e | ma}. This group is Abelian and the N irreducible rep- 
resentations are one-dimensional. The complex co-ordinates 
are each a basis for one of these irreducible representations. As 
the co-ordinates are complex there is now no need to demand 
real representations. The operator {e | ma} is defined so that 


{e | ma}xn = Xn—m 
and by operating on (4.42) 


{e | mayRi = ¢tn'™R, 


i.e. {e | ma}Ry = ¢tkmaR, 


The irreducible representations may be labelled by the wave 
vector so that in the kth representation {e | ma} is represented 
by «ma, The representations of the translation group will be 
considered in more detail with reference to the Bloch theorem 
in Chapter 5. 


§ 4.7 Quantum theory of lattice vibrations 


Up to this point the treatment of lattice vibrations has been 
classical. Certainly the wave vector takes discrete values so that 
the frequency band is quasi-continuous but this is not a quan- 
tum effect. However, it is a relatively simple task to reformulate 
the problem in quantum terms. The quantum energy levels are 
necessary to give the theory of specific heat. 
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It has been shown that for a linear crystal the total energy is 
N-1 
Paks ie 
E=- a (Rit Rr + wi2Ri*R1) (4-49) 


fa 


2 
N real normal co-ordinates are defined by 


Qu= (Ri + Ri) Qu = a(R — R) (4.50) 


and the total vibrational energy of the lattice is 


> > (Q2i + w12Q?41) (4.51) 


Note that the sum is over the positive half of the Brillouin zone 
excluding the origin. The single / = o mode gives the transla- 
tional energy. It is not difficult to show that the new inde- 
pendent co-ordinates represent standing wave modes. From 
(4.44) and (4.50) the atomic displacements due to the single 
co-ordinate Qi are 


Xn = ft Cos (Fin Qu sin (wit + €) (4.52) 


The momentum conjugate to Qu is Pu = 5Or = mOiu and 


the classical Hamiltonian is 
Nai 
Peay 202 
HH = > aE — tf un Qu (4.53) 


The Schrédinger Hamiltonian corresponding to (4.53) is 


#= 5 >|- x a + Fon? Qu | (4.54) 


This Schrédinger equation is separable in the normal co- 
ordinates and the many particle wave function is the product of 
(N — 1) single particle wave functions satisfying 


[am Tout Feta la = Be (458) 
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That is, in the harmonic approximation the crystal may be 
treated as N — 1 independent harmonic oscillators with fre- 
quencies given by the classical dispersion relation (4.6) (or 
(4.20)). It is well known that the solutions of (4.55) are the 
Hermite functions and the single particle energies are 


Eu = (nu + 4)hor (4.56) 


where mj is a positive integer or zero. The total energy of the 

lattice vibrations is obtained by summing (4.56) over the N — 1 

vibrational modes and the lowest allowed energy is } fiwi. 
I>0 


This is called the zero-point energy and is an expression of the 
uncertainty principle. 

It is possible to transform to real normal co-ordinates which 
represent running waves, but to do this it is necessary to use a 
contact transformation which mixes R; and R;. The new co- 
ordinates are defined by 


Q =a] (RHR) + F(R —R)] (457) 
Conversely 


R=YQ+ e+ 2-ey]  (s8) 


Clearly Q: satisfies equation (4.36) (see 4.45) and is given by 
(4.37). From (4.44) and (4.58) the displacements due to the 
single co-ordinate Q; are 


Xn = a sin (— Fyn + wit + a) (4-59) 


In terms of these normal co-ordinates the classical Hamiltonian 
is 


N-1 
ys Fo + ZarQit| with P; = m0; (4.60) 
= Lam ' 2 
i 


The quantum mechanical Hamiltonian is again separable in 
the normal co-ordinates and the many-particle wave function 
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is the product of N — 1 Hermite functions. (The / = 0 mode 
is omitted). Once more the modes are divided into degene- 
rate pairs and the allowed lattice energies are the same as 
before. 

There is an interesting alternative quantum description of 
the lattice vibrations. Regarding the real running-wave normal 


co-ordinates as Schrédinger operators, new operators may be 
defined by 


_ fatig, . ; [2 
a=,/ ah Q +i Shae (4.61) 


wm A I 
att = /——Q, —1i | ———P, 
’ ah Q: 2hwym : 


a;* is the Hermitian conjugate of the operator a; (See Landau 
and Liftshitz Quantum Mechanics). Equations (4.61) may be in- 
verted to give 


= [ast +a;) Pi=i MOTT ays — az) (4.62) 


The quantum mechanical Hamiltonian for a single mode is then 


H, = fet aytay + ayayt) (4.63) 


With the infinite set of Hermite function solutions of the 
travelling-wave oscillator as basis the operators may be repre- 
sented by matrices. Obviously the Hamiltonian is represented 
by the diagonal matrix. 


[Hi] = tr a (4.64) 


oO OO él 
O wie O 
vis O © 


The eigenfunction with energy (n + 4)fw1 is represented by a 
column matrix 
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momentum’ requires further discussion. This is referred to 
again in Chapter 6. 

The higher non-harmonic terms in the Taylor expansion 
(4.27) produce interactions between the phonons, and the 
phonon has a finite lifetime (van Hove, 1959; Peierls, 1929). 


§ 4.8 Three-dimensional lattices 


Apart from the general discussion of normal modes (Section 
4.4) the above theory is concerned with linear lattices. No 
attempt is made in this text to give the complete three-dimen- 
sional theory but some important results are indicated below. 
In analogy to 4.5 there are solutions in which the atoms vibrate 
with the same frequency. 


Xin = Ainet ott k-Ra) (4-73) 


xin is the displacement in the zth direction of the atom at 
Rn = mai + neaz + nsa3. The crystal is assumed to be com- 
posed of identical microcrystals with edges Nai, Naz, Nas. 
These periodic boundary conditions require 


N(k. a) = har I, integer 


For simplicity the lattice basis vectors will be assumed ortho- 
gonal and the components of the wave vector are 


2 
ki = lina, (4.74) 
There are N distinct k-values and these may be chosen so that 
_ Wer <h< ‘- = Nodd = (4.75) 


These values of & lie within a rectangular parallelepiped with 
corners at 


k = 4[-+b: + be + Bs] (4.76) 


where the vector 5; is parallel to a; and b; = 2m/a;. This 
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is the product of N — 1 Hermite functions. (The / = 0 mode 
is omitted). Once more the modes are divided into degene- 
rate pairs and the allowed lattice energies are the same as 
before. 

There is an interesting alternative quantum description of 
the lattice vibrations. Regarding the real running-wave normal 


co-ordinates as Schrédinger operators, new operators may be 
defined by 


_ fitig, [3 
a | Oh Qi +i ahcoumn | 482) 


wim ‘i I 
agha f ONG e fp 
: 2h Q: 2hwim : 


ay+ is the Hermitian conjugate of the operator a; (See Landau 
and Liftshitz Quantum Mechanics). Equations (4.61) may be in- 
verted to give 


a= A/ salar" + ar) Pr= ‘| POH ays — a1) (4.62) 


The quantum mechanical Hamiltonian for a single mode is then 


A, = Pel aytay + aja;*) (4.63) 


With the infinite set of Hermite function solutions of the 
travelling-wave oscillator as basis the operators may be repre- 
sented by matrices. Obviously the Hamiltonian is represented 
by the diagonal matrix. 


{i} =fwift o o.... (4.64) 
o #0 
o o 2 
ee 


The eigenfunction with energy (n; + 4)fiw; is represented by a 
column matrix 
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n> = (4.65) 


--0 0 
ewe | 
mH O 


8) | ni -—~ I 
ni 


| m-+i 


The matrices representing Q; and P; are Hermitian and (see 
Sneddon’s book Special Functions of Mathematical Physics and 
Chemistry for properties of Hermite functions) 


h OVI O O 0O.... 
wi 2wm|r/I 04/2 0 Oo 


02 0 V3 0 


a ~ sone .66 
P= i ve as ais ; i 
o 2 0 —V2 


OH 


Then 
ar=[o +/1 ro) 
0 0ov2 Oo 
0 0 0V3 
airs oO oO re) 
VI oO oO 
0 1/2 oO 


0 0 V3 
(4.67) 


Observe that a; and a;*+ are not Hermitian. By matrix multi- 
plication 
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arttaa=[o0 0 0 O.. 
o 1 0 0 
0 0 2 0 
(4.68) 
and ajaj+ — ay+ta; = 1. The Hamiltonian (4.63) becomes 
A = hey (naz*ay -+- 4) (4.69) 


in agreement with (4.64). The operators a; and a;+ have the 
property 


alu>=Vm|m—1> at|m> = Vm +1lm+1> 
(4-70) 
and so may be regarded as quantum annihilation and creation 
operators, respectively. 
The many-particle state is completely described by the list of 
quantum numbers m and the state function may be written 
m...nmv—-1> with 


arta:|m...mv-1>=m|m...nw-1> (4.71) 
my =0,1,2,... 


The total Hamiltonian is obtained by summing (4.63) over / 
(exclude / = 0) and the eigenvalues of the vibrational energy 
are 


E= Dhon(m + $) (4.72) 


The state is defined by the number of quanta of excitation in 
each running wave mode. These quanta may be regarded as 
particles and are called phonons. (Compare the theory of 
photons.) There can be any number of quanta in each mode 
and the phonons obey Bose-Einstein statistics. 

In essence, the crystal lattice is replaced by a non-interacting 
gas of phonons with energies fiw;. It can be shown that in 
collisions with electrons the quantity 4k; manifests itself as a 
type of linear momentum. As f is not uniquely determined for a 
lattice mode (4.7) the strict interpretation of the ‘phonon 
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momentum’ requires further discussion. This is referred to 
again in Chapter 6. 

The higher non-harmonic terms in the Taylor expansion 
(4.27) produce interactions between the phonons, and the 
phonon has a finite lifetime (van Hove, 1959; Peierls, 1929). 


§ 4.8 Three-dimensional lattices 


Apart from the general discussion of normal modes (Section 
4.4) the above theory is concerned with linear lattices. No 
attempt is made in this text to give the complete three-dimen- 
sional theory but some important results are indicated below. 
In analogy to 4.5 there are solutions in which the atoms vibrate 
with the same frequency. 


Xin = Ajnet ott Rn) (4.73) 


xin is the displacement in the ith direction of the atom at 
Rn = mai + neag + nsa3. The crystal is assumed to be com- 
posed of identical microcrystals with edges Nai, Naz, Nas. 
These periodic boundary conditions require 


Mk. a1) = han i; integer 
For simplicity the lattice basis vectors will be assumed ortho- 
gonal and the components of the wave vector are 
2 
hy = lige (4.74) 


There are N? distinct k-values and these may be chosen so that 


~Wed cy <2 Nodd (4.75) 


These values of & lie within a rectangular parallelepiped with 
corners at 


k = 3[4b1 + be + bs] (4.76) 


where the vector 5 is parallel to a; and &; = 2n/a;. This 
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volume is called the Brillouin zone. (See Chapter 5 for a more 
complete discussion.) When there is one atom per primitive 
cell substitution of (4.73) in the three-dimensional analogue of 
4.4, results in a cubic expression for w2 and each root is a con- 
tinuous function of k, periodic in the Brillouin zone. For a given 
value of the wave vector there are three frequencies all of which 
go to zero as k ——> o. These are the three acoustical branches. 
Note that there are 3N? normal modes as expected. At ‘sym- 
metry’ points in the Brillouin zone there may be degeneracies 
between the different branches and it may be possible to 
distinguish a longitudinal mode and two transverse modes. 
However, for a general wave vector there is no degeneracy 
and there is no definite distinction between longitudinal and 
transverse modes. For each branch of the spectrum w?(k) = 
w?(—k), 

If there is a basis of 7 atoms per cell there will be 37 modes for 
each k. The additional modes are called optical modes and do 
not go to zero as k = o. 

The quantization of the vibrational modes follows as in the 
previous section and the crystal may again be regarded as a 
system of independent quantum oscillators. 

Information about the lattice spectrum has been found 
experimentally by the scattering of X-rays and thermal 
neutrons. The neutron method is most appropriate as the energy 
of the thermal neutrons is of the same order as the lattice waves. 


§ 4.9 The lattice specific heat 


Classical theory 


It is well known that classical theory predicts that the heat 
capacity at constant volume is independent of temperature. In 
the harmonic approximation, the vibrations of a crystal com- 
posed of N3 identical atoms can be described in terms of 3N3 
normal co-ordinates each with a characteristic frequency. 
(Strictly six of these frequencies are zero.) From classical 


108 


THE THEORY OF LATTICE VIBRATIONS 


statistical mechanics it is simple to show that the mean energy 
of an oscillator in thermal equilibrium is kT where T is the 
absolute temperature. The total internal energy of the crystal is 


U = 3N%T (4.77) 


(Here & is the Boltzmann constant. Later it will be used also for 
the wave vector). For one mole, N? = 6-025 . 1023 (Avaga- 
dro’s number) and 

U = 3RT (4.78) 
where & is the gas constant. The lattice contribution to the 
molar heat capacity at constant volume is 


oT 

So classical theory predicts that the heat capacity per mole is 
the same for all substances at all temperatures. At room tem- 
perature and above many materials obey this law, but the 
agreement fails at lower temperatures. Experimentally it is 
found that the heat capacity approaches zero as T——> o. Near 
absolute zero the heat capacity of metals decreases as T. This is 
due to the conduction electrons and was mentioned in Chapter 
2. For non-metals the heat capacity approaches zero as T, This 
effect is characteristic of the lattice vibrations. 


Cy = (Sr) = 3R ~ 6 cal/deg mole (4.79) 
v 


Quantum theory 


In Section 4.7 it was explained how the normal mode oscillators 
are quantized with energy levels 


Ey = (m + $) hor 
There is no restriction upon the number of phonons in a given 
mode and they obey Bose-Einstein statistics. On average the 


number of phonons in the mode with angular frequency w at 
the temperature T is 


fi(w) = 1/(cto? — 1) (4.80) 
The mean energy of the normal mode is 
E(w) = [fi(w) + $]iw (4.81) 
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The total internal energy is obtained by summing over all 
modes. As the frequencies lie close together it is convenient to 
replace the sum by an integral. A frequency distribution is 
defined so that g(w) dw is the number of allowed frequencies 
between w and w + dw. This function must satisfy 


[e(o) dw = 3N8 (4.82) 


g(w) is a sum of contributions from the different branches of the 
vibration spectrum. The total internal energy is 


U = [e(w)E(w) do (4-83) 


The difficulty arises in the derivation of g(w). 

Einstein considered the simplest possible model (1907, 1911). 
He assumed that all the modes have the same frequency wo. 
The internal energy is then 


U = 3N8E (wo) 
= 3N3 al 1s 3 [hes (4-84) 
The heat capacity at constant volume is 
oU Or\2 eox/T . 
Co = aT => 3N°K( 7) [ealT — TP | (4.85) 


where the Einstein temperature @z is defined by fiwo = kOz. Ox 
for a material is chosen to make equation (4.85) give a reason- 
able fit to the experimental results and typically is of the order 
of a few hundred degrees absolute. Einstein’s formula is in fair 
agreement with experiment for JT’ >> 6z but at low tempera- 
tures when J’ << 6g the agreement is poor. Instead of giving 
a T8 variation 


C= BN8K( SE) een” T<<0_ (4.86) 


To obtain the heat capacity at very low temperature it is only 
necessary to know g(w) for small frequencies as the mean 
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energy E(w) becomes small when w is large. The motions of 
adjacent atoms are certainly not independent for these im- 
portant long wavelength modes. The Debye model (1912) 
successfully explains the T® variation at low temperatures. In 
the Debye theory the acoustical spectrum of a solid is treated as 
if the material is homogeneous and all the lattice waves travel 
with the same phase velocity; the ordinary velocity of sound. 
This is clearly reasonable for very small frequencies and these 
are the important modes at low temperatures. For a diatomic 
solid, the Einstein model may be used for the optical branches. 

A cut-off frequency is defined so that the total number of 
acoustical frequencies is 3N%. The number of modes in a given 
branch of the spectrum with wave vector less than k is (from 
4-74) : 

M (k) = Soke [SP (4.87) 

where Q is the crystal volume. So the number of modes with 
wave vector between k and k + dk is 


k2 
N(k) dk = —,dkO (4.88) 


and the number of modes with frequency between w and 
w -+- dw is 


g(w) dao = N(k)SAde (4.89) 
7 depends upon the form of the frequency spectrum. Debye 
took 
k= z i.e, ot = a (4.90) 
where C is a constant velocity and so 
ew) = “ay (4.91) 


For a three-dimensional crystal there are three acoustical 
branches. If it is assumed that the medium is isotropic (no 
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crystal is) so that the spectrum can be described by one longi- 
tudinal and two degenerate transverse branches with phase 
velocities C; and C; respectively then 


w2QO/ 1 2 
g (w) = las a ca) (4.92) 
The cut-off frequency wp is defined by (4.82) so that 
I 2\Q oN3 
(as so ta)eqt ae (4-93) 


g (w) 


Fic. 4.16 Debye frequency distribution 


Note that there is a common cut-off frequency for all three 
branches and so the longitudinal and transverse branches have 
different cut-off wavelengths. 

From Brillouin zone theory it would seem more appropriate 
for the branches to have the same cut-off wavelength and hence 
different frequencies. Brillouin discusses this in his book Wave 
Propagation in Periodic Structures. The internal energy is from 


(4.83), (4.92) and (4.93). 
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_ 9N8 
Te eppe | hw Excess a | dw (4.94) 
and the Debye heat capacity is 
T 3 8p/T exx’ 


where the Debye temperature is defined by iwp = k 0p. The 
integral (4.95) can not be expressed in terms of known functions 
but has been tabulated (Landolt—Boérnstein tables). 
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Fic. 4.18 Debye heat capacity 


8p may be calculated by a suitable average of sound velocities 
(Blackman, 1941) and is not really as simple as (4.93) would 
imply. It is not difficult to see that the Dulong—Petit law (4.79) 
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is reproduced at high temperatures. Of course the correct 7? 
variation is given at very low temperatures when T' << 8p. 
At these low temperatures only the long wavelength modes are 
excited (4.80) and for these the Debye distribution is exact. If 
the upper limit on the integral is made infinite 


Cy = 3Ne( 5.) 42 
D/ 5 


T<< @p (4.96) 
At intermediate temperatures the fit is not so good. A standard 
test for the accuracy of the Debye specific heat function (4.95) 
is to calculate @p at different temperatures from the observed 
heat capacity. If the Debye theory was exact 6p would be a 
constant. Often the variation in the Debye temperature is as 
little as 20%. 


§ 4.10 Anharmonic effects 


There are important phenomena for which it is necessary to 
consider the cubic terms in the potential function expansion 
(4.27). Examples are: 


(a) thermal expansion 
(b) thermal conductivity and 
(¢) specific heat at high temperatures. 


In the harmonic approximation, the mean displacement of 
the atoms from their positions of equilibrium is zero. When the 
cubic terms are included it is found that there is an average 
displacement proportional to the absolute temperature. (See 
Kittel’s Introduction to Solid State Physics.) 

Also, in the harmonic approximation the lattice waves travel 
freely without attenuation. This does not happen in reality 
owing to anharmonic terms and also due to crystal imper- 
fections. The phonons interact with each other. 

At high temperatures (T > 0p) the specific heat deviates 
slightly from the Dulong—Petit value in that it includes a linear 
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term. This may be accounted for by considering the cubic and 
quartic terms in the potential expansion. 

The reader is referred to Peierls’ book Quantum Theory of 
Solids for further discussion. 


Problems for Chapter 4 


1. Write down the one-dimensional analogue of equation 
(4.28). The force on the nth particle is obtained by taking the 
negative derivative of the potential energy with respect to the 
displacement xn. Hence obtain 


Min = —> famXm 
m 


By considering the forces between particles explain why 


Snm = fmn = A(|n — m|) 
and > fam = 0 


Assuming nearest neighbour interactions only, derive equation 
(4.4) with 2« = A(o). Why is « positive? 

2. Find an expression for the group velocity of the lattice 
waves in a line of identical atoms. Show that this velocity is 
zero at the edges of the Brillouin zone. (In a lattice with no 
absorption the group velocity is equal to the ‘energy velocity’.) 

3. Consider a line of four equally spaced particles of equal 
mass m, the end ones being in fixed positions. Write down the 
equations of motion for longitudinal vibrations of the two free 
particles assuming that the force between any two is propor- 
tional to their relative displacement and that nearest neighbours 
only need be considered. Show that the angular frequencies of 
the two normal modes are in the ratio. 
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Write down expressions for the kinetic and potential energies 
and verify that the mass weighted normal co-ordinates are 


Qi = [Mees —*2) Qo= NE + x2) 


where x, and x are the displacements of the two free particles. 
4. The energy of a classical one-dimensional harmonic 
oscillator with angular frequency w is 


where m and p are the mass and momentum of the particle. 
When in thermal equilibrium the average energy of the 
oscillator is 


E = [[® exp (—EkT) dp ax/f fexp (—E/kT) dp dx 
The integration is over phase space. Show that E = k7. (You 
may use [(xe= dx = +/n/2 and | “ett de = Vr. 


5. Consider one branch of the frequency spectrum. The locus 
of all points in the Brillouin zone with a given frequency is 
called a constant frequency surface. Explain why the spacing 
between two such surfaces with frequencies w and w + dw is 
dk = dw/| gradjw |. From equation (4.74) the density of 
allowed k-values is Q/(27)® where Q is the crystal volume. 
Show that the frequency distribution is 


Q I 
Ew) = ceayel | Tezadea co 
where the integration is over the surface of constant frequency. 


One approximation for the frequency in the vicinity of the 
bottom of an acoustical branch is 


a= y (kz? + ky? + k,2)% 
Show that g(w)aw? and determine the constant of propor- 
tionality. (Compare equation (4.91) in the Debye model.) 
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6. Consider a nondegenerate branch of the lattice-spectrum. 
As w? is an eigenvalue of a matrix whose elements are analytic 
then w? is also analytic, i.e. w? is continuous and has continuous 
derivatives. Write down the relation between grad,w? and 
grad,;w. Show that the group velocity is necessarily zero at a 
maximum or minimum of w? except at the bottom of the 
acoustical branches (Fig. 4.10). Employ the periodicity of w2 
in k-space to deduce that the normal component of gradxw is 
zero at the boundary of the Brillouin zone. (This last result is 
not true if two branches are degenerate at the zone boundary.) 

7. Consider a monatomic solid. A model similar to Debye’s 
for calculating the lattice specific heat is as follows. It is again 
assumed that the phase velocity in a given branch is a constant 
and that the transverse branches are degenerate. Suppose a 
unique cut-off wavelength is defined by replacing the Brillouin 
zone by a sphere of equal volume. Then different cut-off fre- 
quencies are defined for the longitudinal and transverse 
branches so that each branch has N? modes. (In the Debye 
model the imposed condition is that the sum of all modes is 
3N8.\ Define characteristic temperatures 9; and 9; for both the 
longitudinal and transverse branches and show that the internal 
energy is 


a U = (Nk) r[D(+) Es 2D() | 
where D(y) -3"l5 = =P i} dx 


8. Consider a pair of atoms, You may assume that when the 
atoms are displaced a distance x from their equilibrium separa- 
tion the potential energy is, 


V (x) = ax? — bx8 — cx' 


Explain the significance of the terms in this expansion. The 
average displacement at a temperature T is 


—2a 


CO [- o} 
z= | xe-VIkT Ax i | eVkT Ax 
a) 
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Show that for small displacements # = 3k7b/4a8. (You may 
use the integrals given in Problem 4.) Note that the anharmonic 
terms are necessary for a non-zero thermal expansion. 

g. The Schrédinger equation for a one-dimensional harmonic 
oscillator of mass m and frequency w is 


(BE Sey— a 


By a suitable transformation of the form x = ky (k constant) 
obtain 


d 
a + (A — »%)b = 0 A constant 


Quadratically integrable solutions exist only if A = 1 + 2n 
where 7 is a positive integer and the solutions are multiples of 
the Hermite functions Yn(_y). 

These are orthogonal and satisfy © 


i “Wall's dy = 2%! 4/n 


Given the recurrence relations 


on n—1 = Tn + Pa’ 
ay, = anFr-1 + Pat 


show that 
mAnt+i 
[¥m¥n dy = ie I)!\/rm m=n+1 
-  1an-l(n)la/a m==n— I 
and 
re) 1) mA~An+i 
i Yn sn’ dy = fers + 1)l/r m=n+1 
i Qt—-l(n)l4/a m=n—I 


Multiply 'n(_y) by a suitable constant to normalize the wave- 
function ¢(x) to unity and confirm the matrix representations 


(4.66). 
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Electron Energy Bands in Solids 


§ 5.1 Introduction 


The lattice periodicity of a crystal has important consequences 
for the dynamics of the electrons. At infinite separation, NV 
atoms considered as a system, have the same energy levels as a 
single atom except that the degeneracy of each level is multi- 
plied by N. As the atoms are brought together to form a crystal 
the degeneracy of the levels is removed by interactions. Each 
level of the isolated atom spreads out to form a band of allowed 
electron energies. Between the bands there may occur forbidden 
energy gaps. The resulting band structure to a large extent 
determines the electrical properties of the crystal. 

From an electrical viewpoint a material may be classified as a 
metal, semiconductor or an insulator by the value of its re- 
sistivity. At room temperature metals have a resistivity of about 
10-6 ohm cm, semiconductors have values in the range 
10-? ~ 108 ohm cm, and insulators have resistivities greater 
than about 101° chm cm. The main factors which decide this 
classification for a material are the valencies of the constituent 
atoms and the magnitude of the energy gaps between the 
bands. 
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§ 5.2 The reciprocal lattice 
Each Bravais lattice has an associated reciprocal lattice, defined 
by the three basis vectors bi, be and bs with 

by . ay = 2784 (5-1) 


a1, 42, a3 are the basis vectors of the direct (Bravais) lattice. 
Vectors 5; that satisfy these relations are 


2 
hi = 7, (a A as), be = V, (a A a1) (5-2) 
bs = a A @2) 


where Va = a1. (az | as) is the volume of the primitive cell 
of the direct lattice. 

The reciprocal lattice is composed of the end-points of all the 
vectors 


Km == mzyb1 + mabe + msbs (5-3) 
m1, M2, msg integers 
The primitive cell has a volume 
Vo = bi. (b2 A bs) (5-4) 
and it can be shown that 
Vo = (27)*/Va (5-5) 
A function f(r) which is periodic with the period of the direct 
lattice may be expanded in a Fourier series in the reciprocal 


lattice vectors Km. To see this observe that the reciprocal lattice 
vectors have the important property 


Km. Rn = 271 (5.6) 


with / an integer or zero. The converse of (5.6) is also true. That 
is if a vector K is such that the scalar product K. Rn is an in- 
tegral multiple of 27 for all Rx then K is a reciprocal lattice 
vector. If f(r) is given by 


F(t) =X Amettn-t (547) 
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then 
S(r + Ra) = DY AmethKm «0 + Rn) 
km 
and as 
eiKm «Rn — 5 
then 


S(r + Ra) =f(r) 
The reciprocal lattice of a reciprocal lattice is the original 
direct lattice. This can be verified by considering the vector 


b 


I 


a,(also b,) out of the paper 


Fic. 5.1 Direct and reciprocal lattice vectors 
2m(bz A bs)/V» and showing that it is equal to a1. From (5.2) 
Fi lba A bs) => (as A as) A (ar A a) 
= 7; {la - 2 /\ @3)a1 — (a1. a2 A ai)as3} 
= ay 
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The reciprocal lattice is invariant under the same point group 
as the direct lattice. For if « is a point group symmetry operation 
which leaves the direct lattice invariant then «Rn is a lattice 
translation if Ry is. From (5.6) 


Km. (aRn) = 2al l interger 
i.e. (21Km) . Rn = 27l 


The latter equation shows that («-!Km). Rn is an integral 
multiple of 27 for all Rx and so («-1Km) is a reciprocal lattice 
vector. 

This final result demands that a direct lattice and its reci- 
procal lattice belong to the same crystal system. In fact ten of 
the fourteen Bravais lattices are their own reciprocals. Examples 
which are not are the face-centred cubic and body-centred 
cubic lattices which are the reciprocals of each other. 


§ 5.3 The Bloch theorem and energy bands 


One-dimensional lattice 


This central theorem shows how the allowed energy eigenstates 
of a perfect crystal may be classified by a real wave vector k 
and will first be proved for the one-dimensional case. Later a 
group theoretical proof will be given for three-dimensional 
lattices. 

The Schrédinger equation for a line of identical atoms is 


— ESE + Vey = (5,8) 
with V(x — a) = V(x) 


Suppose ¥i(x) and ye(x) are two real linearly independent 
solutions each satisfying prescribed boundary conditions and 
belonging to the energy EF. Any other solution with energy E 
may be expressed as a linear combination of them and in par- 
ticular yi(x — a) and e2(x — a) are solutions with 
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ya(x — a) = aarpi(x) + aroye(x) (5.9) 
yro(x — a) = aerpi(x) + aeaibe(x) 
Key real 


If the ¥1 Schrédinger equation is multiplied by y2 and the ge 
equation by y then subtracted and the result integrated by 
parts then | 


W(x) = hh = — te SB = constant (5.10) 


W(x) is the Wronskian of the solutions. From (5.9) and (5.10) 
W (x _— a) = (11022 ~ a12%21) W (x) (5.11) 
and so (x11%22 — X120%21 ) = iI 
It is possible to find solutions ¢ which are linear combina- 
tions of ys, and #2 so that they are multiplied by a constant on 
translation by a. 
1,e. d = Bid + Baye B1, B2 may be complex 
with p(x — a) = Ad(x) (5.12) 
A is a constant. From (5.12) and (5.9) 
(11 — A)Bi + a2182 =o 
41281 + (22 — A)Be =o 
The condition for a non-trivial solution is 


“11 — A 21 =0 
o12 a22 — A 
i.e. Az — A(o11 + a2) +1=0 


The product of the two roots of this quadratic is unity and as 
(a11 + ag2) is real their sum is real. A real function of the 
energy is defined by 


w(E) = $(a11 + a2) = aa ct As 


where A; and 22 are the two roots. It is convenient to write 


Ai oo etka Ag = etka 
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where f is not necessarily real. Then 


p(E) = cos ka 
There are three distinct cases to consider. 
(i) | w(E)| >1 


X is not real but the roots of the quadratic are real. 
Ar = +e84 Ag = ea 
8 real 


Both these solutions are divergent for an infinite crystal and do 
not represent allowed eigenstates. 


(ii) | u(E) | <1 


k is real and the roots are complex and lie on the unit circle in 
the complex plane. They represent allowed eigenstates and the 
wave functions are denoted by & and are called Bloch waves. 


u(x — a) = e-thay(x) (5-13) 


This is the Bloch theorem. Note that the states for +k are de- 
generate and that & is defined only to within a multiple of 
2m/a. It is conventional to restrict k to lie in the range 
—ala <k <-a/a and there will be a series of wave functions 
for each & transforming as in (5.13) but with different energies. 


(ii) | (Z) | =1 


In this case the quadratic has a double root equal to +1. This 
is an allowed solution. 
Kramers (1935) studied the form of »(£) for a non-degenerate 
band (Fig. 5.2). 
In particular 
Lim p(£) = + 0 
E—>—© 


p(E) ——> cos (a\/E) as E —> 


In between these limits u(Z) oscillates between values greater 
than one and Jess than one and so there are alternate allowed 
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and forbidden bands of energy. This essential property is main- 
tained in three-dimensional crystals. 


Three-dimenstonal lattice 


A rather different approach will now be used to prove Bloch’s 
theorem for a three-dimensional lattice. The one-electron 
Schrédinger Hamiltonian is 


f2 


KH = — — 
2m 


Mer VA) (5-14) 
The crystal potential V(r) is unchanged if the point of observa- 
tion is displaced by a lattice vector Rn. 


V(r — Rn) = V(r) (5-15) 
Rn = ma. + meag + nsas 1, Nz, Ng integers 


The substitutional operator representing this translation is 
{e | R2}. 
fe| RV(r) = V(r — Rn) (5.16) 


For an infinite crystal the number of lattice translations is 
infinite. A real finite crystal is of course not truly invariant 
under {e| Rn} but if surface effects are negligible it may be 
taken to be so. To simplify the discussion it is convenient to 
keep the number of distinct translations finite by the introduc- 
tion of periodic boundary conditions. 

Consider a microcrystal of size Nai by Naz by Nas. Two 
points of view may be taken. Either consider an infinite number 
of identical microcrystals stacked together or imagine opposite 
sides of the crystal to be joined in some sense. The latter view- 
point is easily conceived in one-dimension in the form of a ring 
of atoms or in two dimensions as a torus surface but there is no 
such physical picture in three dimensions. The three transla- 
tions Na:, Naz, Naz are considered to produce no change at all. 


{e | o} = {e | Nai} = {e | Naz} = fe | Naz} (5.17) 


The periodic boundary conditions may be used whenever in- 
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terior properties of crystals are studied and the surface to volume 
ratio is very small. 

The translation operators now form a finite group of N? ele- 
ments. This group is of course a sub-group of the full crystal 
space group. The translations commute with one another, so 
forming an Abelian group, and consequently all the irreducible 
representations are one-dimensional and the number is equal 
to the order of the group N%, 

Clearly the Hamiltonian is invariant under the translation 
group and so the energy eigenfunctions must transform accord- 
ing to one of the irreducible representations of the group. 
Another way of stating this is to say that the Hamiltonian and 
the translation operators commute with each other and so wave 
functions may be chosen to be simultaneously eigenfunctions 
for # and all the operators {e | Ry}. 

Let y be the basis for one of the irreducible representations 


{e | ai}(r) aT Avb(r) i= I, 2,3 (5.18) 


where 2; are constants which may be complex. The periodic 
boundary conditions demand 


AYN = dgN = AgN = 1 (5-19) 
Each A; must be one of the N roots of unity. 
Ae = eo tanry/N fe = 1,2,...N (5.20) 


For a general translation {¢ | maz + m2az + ngas} A is given 
by 


e-t a (71%, + Teg + £323) (5.2 1) 


This may be written e-#-®» where the real ‘wave vector’ k 
takes the discrete values 


k= = (rab + reb2 + 13b3) r= 1,2,..5N 
(5-22) 


bi, be, bg are the reciprocal lattice basis vectors. These k-values 
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lie very close together and for many purposes are taken to form 
a continuum. 

In this way, the irreducible representations and hence the 
energy eigenfunctions may be characterized by a wave vector 
and for these functions 


fe | Ru}pe(r) = e-*- Fuby,(r) (5-23) 
This is Bloch’s famous theorem. 
A function satisfying this condition is 
Par) = efh-ruy(r) (5.24) 
if u,(r) has the lattice periodicity. 

If the wave vector of an eigenstate is defined by (5.24) then 
it is obvious (from (5.6)), that & is not unique since & and 
k + Km are equivalent. 

etRn(k + Km) = etRn-k for all Ra (5.25) 


To label k uniquely it is necessary to restrict k to a reciprocal 


Xx x x 
x X xX 
x x x X 


Fic. 5.3 Brillouin zones of square reciprocal lattice. The 
dotted square shows an alternative primitive cell 


lattice primitive cell. This value of k is called the reduced wave 
vector. The symmetric primitive cell, called the Brillouin zone, 
is obtained by bisecting the lines joining k = o to the nearest 
reciprocal lattice points. The analogous cell in the direct lattice 
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is the Wigner—Seitz cell. Each point exterior to the Brillouin 
zone is equivalent to some interior point in the sense of (5.25) 
and similarly, a point on the surface of the zone will be equiva- 
lent to at least one other point on the surface. 

There are N%k-vectors of the type (5.22) inside the Brillouin 
zone and this gives correctly the total number of distinct irre- 
ducible representations of the translation group. 


E(k) 
Band 2 { 


Band |! 


——— ae ee oe oe ln 


0 Kk 
pig 


Fic. 5.4. Energy band structure 


The energies of the electron states may be regarded as a 
function of the ‘ quasi-continuous’ k-vector. In this reduced 
zone scheme, this function is multi-valued and is normally 
written E,(k), n denoting the branch. It may be shown that 
within the Brillouin zone, the energy is a continuous function of 
k with a continuous derivative. A single continuous branch is 
called an energy band. Discontinuities in the energy may occur 
only on the surface of the zone. 

The restriction of & to lie within the Brillouin zone is only a 
convention. Alternatively A may be allowed to take all values 


131 


SOLID STATE PHYSICS 


and then the energy £(k) will be a periodic (still multivalued) 
function of k. 


En(k) = En(ke + Kn) (5.26) 


If the periodic boundary condition is relaxed then there are 
other possible irreducible representations. The operator {e | Rn} 
may be represented by e#-%» where » is a complex constant 
vector depending upon the basis function. 


{e | Ra}py(r) = ev-Rnp,(r) 


The commutation relations of the translation group are 
satisfied, but these solutions do not correspond to physically 
acceptable wave functions as they diverge. The eigenfunctions 
in the energy gap between bands are of this form. 

The major problem in band theory is to find as good an 
approximation as possible to Z(k). An important quantity that 
can be derived is the density of states N(E) defined so that 
N(E) dE is the number of states in a band in the range E to 
E + dE. Consider a surface in k-space on which the energy has 
the constant value E and construct a right cylinder with base 
area dS on this surface and with length dk. The number of 
states in this cylinder is (5.22) Q dk dS/(27)3 where Q is the 
crystal volume. As dk = dE/| grad,£| then the number of 
states per unit energy is 


Q 
NE) = aa ff era (5.27) 


where the integration is over the surface of constant energy. 
Each of these states may be occupied by two electrons. Clearly 
N(£) may have a singularity when the integration includes a 
point at which | grad,# | = 0. These are called critical points 
and their occurrence has been studied by Morse (1938) and 
van Hove (1953 — lattice spectrum). There are four types of 
analytic critical points about which the energy can be ex- 
pressed as a Taylor series in & with the first derivative terms 
omitted. These are a maximum, minimum and two kinds of 
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a b c d 


lf 
(a)— Minimum W(EI@E* 


(b) 


(o) fSadate points— discontinuity in slope 
Vy, 
(d)— Maximum W(£) a (£,-£)* 


Fic. 5.5 Singularities in the density of state function 


saddle point. It can be shown that a non-degenerate band must 
include singularities of all four types. 


§ 5.4 Symmetry properties of the energy function 


The operations which carry a crystal into itself form the space 
group. The Seitz notation for a general rotation-translation 
operator is {a | t} and has previously been discussed (3.8—-3.11). 
t is not necessarily restricted to the lattice translations. For 
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crystals possessing screw axes or glide planes, t may be a 
fractional lattice displacement. 

The point group of a crystal is the collection of rotations (in- 
versions or reflections) «, which appear in the elements of the 
space group. There may exist rotations «’ which always occur in 
the form {«’ | t’} where tf’ is not a primitive translation. Even 
though {«’ | o} is not a member of the space group the rotation 
a’ is a member of the point group. It is quite possible that no 
atom site in a crystal has the full point group symmetry. For 
example, the point group of the diamond crystal is the full 
cubic group mgm but the carbon atoms are in an environment 
with the reduced tetrahedral symmetry 43m. Diamond, of 
course, has a screw axis. A space group which does not involve 
non-primitive translations is called symmorphic and in this case 
the point group is a sub-group of the space group. 

In Section 5.3 it has been shown how the lattice periodicity 
means that the energy eigenfunctions may be specified by a 
wave-vector k. If the crystal point group contained only the 
identity (as in a triclinic lattice without even inversion sym- 
metry) and the complex conjugation operator (Section 5.5) is 
ignored, then all the eigenstates would have different energies. 
In general though, the elements of the point group produce 
degeneracies between different k states. 

The operator {« | 0} transforms a Bloch wave with vector k 
into the Bloch wave with vector «k. From (5.24) 


{x | o}fa(r) = eth. (oy (atv) (5-28) 


The scalar product k . (ar) is invariant if the same rotation is 
applied in both k-space and real space. 


(ak) .r =k. (a-Ip) 
{oe | OF (7) = fh) -rugy’(r) = pox(r) (5.29) 


The function U,,’ maintains the lattice periodicity. 
It is left to the reader to show that 


{a | thf, (r) = eek) -ru,,’(r) (5-30) 
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even in the case where f¢ is not a primitive translation. The 
rotation-translation operators commute with the Hamiltonian 
and so all the wave functions given by {« | t}}, must be de- 
generate with #,. The energy spectrum in k-space must have 
the full point group symmetry. 


En(k) = En(ak) (5-31) 


The wave vectors generated from k by the crystal point group 
form a set which is called the star of the k-vector. 

For symmetry points in k-space «k = k for some point group 
operators, These operators form the ‘small group’ of that par- 
ticular k. There is the possibility of degeneracies between bands 
at such symmetry points and then Bloch waves from the 
different bands must belong to an irreducible representation of 
the small group of k. 

A simple two-dimensional example will help to clarify the 
above. A square direct lattice has a square reciprocal lattice. 
If a is the minimum repeating length in the real lattice then 
2m/a is the minimum length in the reciprocal lattice. The Bril- 
louin zone is a square of side 27/a. 

The square point group is 4mm with eight elements. (See Fig. 
3.10.) Each of these operations will take a general wave vector, 
say ki into a new position. The star will include the eight points 
ki to kg (Fig. 5.6). There are six types of special points in the 
zone. The point I at the centre of the zone with k =o has a 
small group identical with the square point group 4mm. The 
star of I is the single point k = o. Degenerate Bloch functions 
at this point must transform according to one of the irreducible 
representations of the square group (Table 3.4). Consequently 
only singly or doubly degenerate states are possible. 

Two elements of the square group, E, my leave the point A 
invariant and so form its small group. (The group m.) The re- 
maining six elements produce the star of A consisting of four 
points on the positive and negative kz, ky axes. 

At first glance it may appear that the point Z is left invariant 
only by the unit element £. This is not in fact true as Z and Z’ 


135 


SOLID STATE PHYSICS 


ky 
eM 
a 
@ 
y Ks x Ko 
/ 
Z e e7 
Xk’ x x 
I ky Kk, 
xX 
stash 2 h 
7 x 
X kg 
x k. 


Fic. 5.6 Square lattice Brillouin zone 


are really equivalent as they differ by a reciprocal lattice vector. 
The small group of Z consists of the two elements E, mz. 


It is left as a simple exercise to the reader to show that the 
small group of 


M i 4mm with eight elements, 
X is 2mm with four elements, 
> is m with two elements. 
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Fic. 5.7 


137 


SOLID STATE PHYSICS 


It is easy to show that an energy band has a relative maxi- 
mum or minimum as it crosses the zone boundary. The point 
ky’ is equivalent to ki and so 


E(ki’) = E(ki) = E(ka) 


It has already been mentioned that the energy function is con- 
tinuous with a continuous derivative and so by taking k close 
enough to the zone boundary the result is obvious. This means 
that, for the square lattice, the normal component of grad,£ (k) 
is zero on the zone boundary. 

Often this result is true in three dimensions. If there is 
sufficient symmetry so that the two opposite sides of a Brillouin 
zone face are symmetrically equivalent then in general the 
normal component of grad,#(k) is zero on the boundary. 
Possible exceptions arise when two bands stick together at a 
zone boundary. This may happen when a crystal contains 
screw axes or glide planes. 


§ 5.5 The complex conjugation operator 
If the crystal point group contains the inversion then clearly 
E(k) = E(—k) (5-32) 


In fact, if spin is ignored it is an easy matter to show that this 
result is always true. The Schrédinger equation is a real dif- 
ferential equation and so if #, is an eigenfunction with energy E 
then ,* belongs to the same energy. 


P(r) == etk-ru,(r) 
py*(r) = ek uy! (1) 


‘ 


with u-,/ = u,* 


The complex conjugate function y,* has the same energy as #, 
but transforms according to #-,. The energy surfaces in k-space 
have inversion symmetry even if this symmetry is not included 
in the crystal point group. 
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In this spin-free case, each k-state may be filled by two elec- 
trons and so in fact each band is two-fold degenerate. When 
spin-orbit coupling is introduced, the situation is more complex 
and will not be considered further. 


§ 5.6 Nearly free electron model 


In this model the crystal potential V(r) is assumed to be very 
weak compared with the electron kinetic energy. This does not 
correspond to the motion of valence electrons in any crystal as 
the variations of potential energy are generally comparable to 
the kinetic energy. Even so, the method is worth describing as it 
is one of the limiting cases (the other is tight-binding) ; and is 
the basis of more sophisticated approaches. (See O.P.W., etc.) 
The approximation will apply to the motion of a beam of fast 
electrons ejected into a crystal. 

For simplicity, the one-dimensional case is discussed. It is 
convenient to choose the zero of energy so that the mean value 
of the potential function is zero. 


je. "V(x) dx = 

i.e i (x) dx = 0 (5-33) 
where a is the lattice period. With this choice of zero, V repre- 
sents the variations of the potential energy. The Schrédinger 
equation to be solved is 


eg i, ee (5.34) 


As Vis assumed to be small compared with the kinetic energy, 
perturbation theory is employed. 

The unperturbed wave functions corresponding to V = 0 are 
the plane waves 


r(x) = Na othe (5-35) 
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These wave functions are normalized over a microcrystal con- 
taining N cells. The unperturbed electron energies are 


E%(k) = h®k2/am (5.36) 


Of course in the free electron model k takes all values. Suppose 
account is taken of the lattice ‘structure’ only by restricting k 
to lie in the Brillouin zone. 


Then 
be (x) se ~/ 


Weta eit EDK) = x(k bn az) 
[A| < . n integer 


The energy bands can be drawn straight away. This is the 
‘empty-lattice’ test (Shockley, 1937). Observe that the bands 
do not join smoothly at the zone boundary. 

When the small periodic potential is introduced the different 
plane waves are no longer independent. However, if V is weak 
it may be expected that each perturbed wave function is com- 
posed primarily of a single plane wave ¢,x and so the real wave 
functions are labelled by k. | 


he = be + Y Avi by (5.37) 


From perturbation theory, the constants A,’(k) are small and 
‘correct to first order are given by 


Av) = Bape (5-38) 
with <k | V|k> = [$u*Vox dx (5.39) 


The perturbed energy, correct to second order is 


E(k) = BMB) +5 Beant (5.40) 


(The first order term is zero by virtue of (5.33).) 
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The real periodic potential may be expressed as a Fourier 

series. 

V(x) = ze Vne-tn = x n integer (5-41) 

nF 

V,* a Ven 
Clearly, the integral 

<k |V|k> =Vaifk—k =21= 

= o otherwise (5-42) 


The state k is mixed only with the states k — nazja by the 
periodic perturbation. The new wave function is now seen to be 


I 


dn(s) = ete] a Aalbers | (5.49) 
with Ap = 1 


Ay, = 


Vn 


ae AM a 
E(k) — 9 k = n=) 


The quantity in brackets is just the periodic part of the Bloch 
wave. The perturbed energy is 


E(k) = Bok) + ¥ —LM@P yy 
nO ROCK) — °( k _ =") 


This final expansion for the energy is satisfactory if the 
Fourier components V, tend to zero rapidly as n increases and 
if there are no degeneracies among the plane wave states of the 


type 
E(k) = E(k — n 22] (5-45) 
The condition necessary for this type of degeneracy to occur 
is (from (5.36)), 
|k| = E nee 
a 


(5-46) 
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At or near a Brillouin zone boundary, where k = nr/a the per- 
turbation expansion is no longer valid. The perturbation 
method fails because it is no longer permitted to assume that 
Ay, is small compared with unity (5.43). The incident and 
diffracted rays are now on an equal footing. Near a zone 
boundary it is necessary to write 


pe = Aode + Andi-n (5-47) 


This expression is substituted into the Schrédinger equation 
(5.34) and the resulting equation multiplied on the left by $;%* or 
$* rn 2%. On integration, the following simultaneous equations 


are obtained. 
Ao[E(k) — E(k)] — AnVa* =0 (5.48) 


—AVn + Aa| E(k) S Bok cs naz) | 95 


For a non-trivial solution, the determinant of coefficients must 
vanish. The resulting quadratic has the two solutions 


E(k) = 4] Bk) + E(k — 02) 
«(ee (to) Herne T] Gow 


At the zone boundary, the two plane wave states have the same 
energy. | 
7 =) = F0 =) 
E ( a\ mE ( a) + | Mn 


e(2)—(2)-1m 


Near & = o the difference in unperturbed energies is so large 
that the perturbation effect is very small. On the other hand, 
the above analysis shows that there is a discontinuity at the zone 
boundary and the energy gap is 


Figure 5.9 (a) indicates clearly how the energies are derived 
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by perturbation of the free-electron parabola. The two points - 
marked C represent the same state in the two schemes. 
For points on either side of the energy discontinuity 


An 
Ao 
Suppose Vj is negative. This corresponds to a periodic potential 
which is negative in the region of each atom (x = na) and so 
attracts electrons. At the zone boundary zr/a the wave function 
belonging to the lower energy E- is 


j= (« a + e132) |y/(2Na) 
= I (ia )os 


+ = il (aa)in Za 


The lower eiegenfunction has its maximum amplitude at the 
atom sites where the potential is most attractive. The higher 
state has its nodes at these points. 

If an external electron beam is incident on the crystal with 
energy lying in one of the gaps it will be totally reflected. 

The generalization of this treatment to three-dimensional 
crystals is straightforward. Energy gaps will occur when 


E°(k) = E(k — K) 
with K a reciprocal lattice vector. 
i.e. |A| =|kK—K| (5-51) 


Geometrically this means that k lies on the perpendicular 
bisector of the reciprocal lattice vector K. (This is simply one 
statement of the Bragg condition.) 

So energy gaps occur at the Brillouin zone boundaries and 
are essentially 2 | Vx |. It is to be expected that at every zone 
boundary (except for certain symmetry points), there will be a 


=1 (5-50) 


Similarly 
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gap between energy bands. In three-dimensions this does not 
necessarily imply that there is a forbidden band of energies 
because the gaps at different parts of the zone occur at different 
energies so that the bands often overlap. 


Fic. 5.10 


The plane wave approach may be expected to be useful only 
if the Fourier components Vx converge rapidly. This is certainly 
not the case. The potential V(r) is due to the periodic array of 
ions and will have singularities at the ion sites. It has Fourier 
components of very short wavelength so that Vx may well be 
considerable for large values of K. 

There are several improved methods such as the ‘ortho- 
gonalized plane wave’ technique. (O.P.W. — Herring, 1940) 
but these are beyond the scope of this monograph. 


§ 5.7 Tight-binding model 


In the other extreme of tight-binding, it is assumed that the 
lattice constant is so large that the wave functions of adjacent 
atoms do not overlap to any appreciable extent. Near the ions 
cores it is to be expected that the crystal wave function re- 
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sembles an atomic wave function. In 1928 Bloch pointed out 
that a crystal wave function with the correct symmetry could 
be constructed from the atomic orbitals. 


Yar) = A(k) > eB B. $(r — Ro) (5-52) 


where ¢(r — Rn) is an atomic orbital centred on the Rath site. 
The constant A is chosen to normalize 4, over the microcrystal 
containing NV? cells. 

The expectation value of the energy for this wave function is 


= fur|—Zve fe Vie) be dr 


f bib, Ar 
E(k) = A2N8y ek . Rie(Ry) (5-54) 


E(k) (5.53) 


with «(R) = [orer — Ry{— Eve + vl 60 ae (5.55) 


Equation (5.54) expresses the periodicity of the energy function 

in k-space and the Fourier coefficient «(R;) is the matrix ele- 

ment of the Hamiltonian between the atomic orbitals on sites 

separated by R:. These Fourier coefficients involve three centre 
V=0 


> “ 
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Fic. 5.11 Atomic and crystal potentials 
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integrals as the crystal potential is normally expressed as a 
superposition of atomic potentials. 

A simplified treatment is given below. The atomic functions 
satisfy the Hamiltonian for a free atom. 


{— Eve + valr)|$ = Bad CEG 


or 
a 
{— V2 4 vin} = [En + (V(r) — Va(r))]¢ 


am 


V, is the potential for a free atom and £, is an eigenvalue. 
If the overlap between different atomic functions is neglected 
A? = 1/N3 and 


e(0) = Ea +a (5-57) 


with c= for - Valo dr 


(Ri) = [o*(r — RLV — Val(r) dr (5.58) 


In fact, some overlap between ¢(r — R:) and ¢(r)is assumed 
so that e(Rz) is not zero but e(R:) will certainly be very small 
except for nearest neighbours. If the atomic wave functions 
have spherical symmetry, the integral (5.58) will be the same 
for all the nearest neighbours. 


i.e. e(Ri) = 8 for nearest neighbours (5-59) 


In this nearest-neighbour interaction approximation the 
energy Is 
E(k) = Ea + a+ Bd e-tk.R (5.60) 
Rn 


As (V — V4) is negative, the integrals «, 8 are both negative 
and may be of the order of a few electron volts. 

For a simple cubic lattice, there are six nearest neighbours at 
the sites (a, 0, 0), efc. 


E(k) = Ea + « + 28[cos ake + cos aky + cos akz] (5.61) 
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(b) Energy plot along section AB 


Fic. 5.12 Simple cubic lattice 
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The Brillouin zone is a cube and Fig. 5.12 illustrates the energy 
contours in the kz, ky plane. 

The band has a width 12 | @ |. 

Each atomic level in the free atom will split into a band of 
energies as the atoms are brought together to form a crystal. 
The N? atoms originally contain N? degenerate levels. The final 
band contains N? different levels. ‘The integral 8 and hence the 
width of the band depends upon the overlap of the atomic wave 
functions, The band will tend to be very narrow if the atomic 
orbitals are closely localized. In practice, the band widths of 
inner core electrons are vanishingly small. 

If two atomic levels broaden so much that they almost inter- 
sect, then it is clear that the method described above may not 
be sufficient. If Ey, — Eg, the separation between two atomic 
levels, is of the same order or less than the band width 126 it is 


Bands overlap 


Ey 
N°- fold degenerate 


0 Interatomic distance 
i] 1 
Demian rs 
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necessary to assume that the wave function is a linear combina- 
tion of atomic orbitals (L.C.A.O.). 


p(n) = ZA;(k) Ze "a bi(r — Rn) (5-62) 


The sum over 7 includes all the atomic wave functions con- 
sidered necessary. For example, a triple degenerate p-state in an 
atom is split into three energy bands degenerate at the centre of 
the Brillouin zone at k = o. The j-sum in (5.62) must include 
all three p-functions. If an s-state occurs close to a triply de- 
generate p-state, as in carbon, then the L.C.A.O. method in- 
volves a sum over all four atomic orbitals. In general, the more 
atomic functions in the sum, the more accurate the final energy 
band will be. The expression (5.62) is used as a trial wave 
function and the energy is minimized. This will result in a set of 
linear equations, the number being equal to the number of 
atomic orbitals and the coefficients A; are the solutions. 

The L.C.A.O. method has certain disadvantages. The bound 
states of an atom do not form a complete set of functions as the 
ionized plane wave states are not included. The overlap inte- 
grals between atomic orbitals on different sites are not negligible 
and this complicates the algebra. Slater and Koster (1954) have 
used the L.C.A.O. method in an interpolation scheme, fitting 
the interaction integrals (5.55) to values either observed em- 
pirically or calculated from the results of other methods. 


§ 5.8 The Wannier function 


The atomic orbitals used in the L.C.A.O. method are not 
orthogonal and calculations involving them are often tedious. 
Wannier (1937) introduced a set of orthogonal localized func- 
tions which may be used to describe a band. The Wannier 
function is defined by a sum over the Brillouin zone of the 
Bloch functions for the band. 


a(r — R,) =— 


/ N38 2 eth Ra by (Pr) (5-63) 
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a(r — R,) denotes the Wannier function on the Rpth site. It is a 
trivial matter to demonstrate that the N? Wannier functions 
centred on different sites are orthogonal. Similarly, the 
Wannier functions belonging to different bands 7 and j are also 
orthogonal. 


fact(r — Rn)aj(r— Rm) dr 
= ya Let Reet Re [a (eau (0) ar 


I P 
= 882m 


as detk «(Ra —Rm) — N38nm 
k 


If equation (5.63) is multiplied by e**’.®» and summed over 
R,, then from the result 


it is simple to verify 
I . 
b(t) = Saya Bee Rear — Rn) (5.64) 


n 


The Wannier functions may be considered as the Fourier com- 
ponents of all the Bloch functions belonging to a band. The 
Wannier functions centred on different sites are all identical. 
From (5.63) 

{e | ~Raja(r — Rn) = 


I 
V NS 


The usefulness of the Wannier function arises because it is 
localized 


2 e~ th. Rngth Ruy, (r) = a(r) 


I 
VN 


Except where r — Ry is small, the exponential term is rapidly 
varying in k-space and the summation is small. The Wannier 


a(r — Rn) = 


> etk(r _- Rn) uy, (r) 
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function on the nth site has an appreciable value only in the 
vicinity of the site. Blount (1962) showed that for bands which 
do not have degeneracies the Wannier function may be chosen 
to fall off faster than any finite power of r. Where there is a 
point degeneracy the Wannier function may decay as r-3 and 
for a line degeneracy it may decay as r-?. 

The Bloch waves satisfy the Schrédinger equation 


Hy, = E(k) by, 
The energy is periodic in k-space and can be expressed as a 


Fourier expansion 


E(k) = Sete Rne(Rn) (5.65) 


If (5.64) and (5.65) are used to substitute for 4, and E(k) in the 
Schrédinger equation then 


> oh Rnta(r — Rn) = > > e(Rm)et*- (Ri -Rma(r — Rz) 
Ra Ri Rm 


This equation may be multiplied on the left by e-t#.® and 
summed over &. It follows that 


Ha(r — Rj) = 2<(Ri)a(r — Ri —R;) (5.66) 


The Wannier functions satisfy a non-diagonal form of the 
Hartree~Fock equations. 

There is an intrinsic ambiguity in the Wannier function 
definition (5.63). The Bloch waves as solutions of Schrédinger’s 
equation are only defined to within an arbitrary phase factor 
e*, However, it can be shown that the phase can always be 
chosen so that 

p* = b-~ 


In this case 
a*(r) = aya B® = Says Se = alr) 


The Wannier function may always be chosen to be real and for 
a non-degenerate band this function decays exponentially. 
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Wannier functions have proved themselves useful in the study 
of impurities in crystals, but in general it is difficult to find the 
Wannier function as it is necessary to know all the Bloch waves 
in a band. However, if the Bloch functions are given by the 
tight-binding approximation the Wannier function is easily 
found. Consider a one-dimensional lattice of like atoms with one 
electron per atom. The Bloch wave is 


u(x) = Ak) > e*mag(x — na) 


where ¢(x — na) is the atomic-orbital on the nth atom and A(f) 
is a constant normalizing the Bloch function over N cells. If the 
atomic orbitals are orthogonal A = 1/4/N and the approximate 
Wannier function is simply the atomic orbital. More generally 


A(k) = [N(QetenaS(o, n))]~# 
with S(m, n) = [orc ~— ma)¢(x — na) dx 
The Wannier function at the origin is 
a(x) = yp bve(*) 
; -1 
ie. ae we] Sermes(o, m) | Yetknag(x — na) 


If the overlap integrals are small and satisfy 


> S(o, m) <1 
m#z90 
then 
(Selma (o, m))-?==1 — 4 PO ai) (o,m) +... 
m mz 


To this approximation the Wannier function is 


id ia A) cama Pea) 
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§ 5.9 The cellular method 


This method was first used by Wigner and Seitz (1933, 1934). 
The cellular method is a completely different approach to the 
band problem than either the plane-wave method or the tight- 
binding approximation, In the latter methods, the wave func- 
tions are chosen from the start to be Bloch functions satisfying 
the conditions imposed by the lattice periodicity, and then the 
ionic potential is considered as a perturbation. The cellular 
method, on the other hand, attempts to find as exact a solution 
as possible near the ion cores and then to approximate to the 
boundary conditions. A more complete discussion will be given 
in Chapter 7. 


§ 5.10 Motion in steady fields 


In general each Bloch wave is complex and describes a state 
with definite mean velocity. Formally the velocity is (in the 
absence of a vector potential) 


h 
V (k) = |e grad de — Yu grad y*] dr 
The integral on the r.h.s. may be evaluated to give 
; . 
V(k) = 5 grad (k) (5.67) 


E(k) is the band energy function. This equation refers to a 
single Bloch state but a simplified wave-packet argument gives 
the same result. The complete time-dependent wave function 
for a Bloch state is 


Pelt t) = u(r) [e-r—F] (5.68) 
A wave packet with peak at ko may be constructed from Bloch 
functions belonging to a given band 


wr, t) =f a(k)pa(r, #) dk (5-69) 


a 
BZ, 


155 


SOLID STATE PHYSICS 


The integral is over the Brillouin zone. If the packet is suf- 
ficiently concentrated so that E(k) may be replaced by its 
Taylor expansion to terms linear in k then 


E(k) = E(ko) + grad,E(ko) . (k — ko) 
and 


ab(r, t) — et [kar 2h) .f a(k) uyet* —*o) [+ — Farad, eh | dk 
BZ. 


The probability density 4* depends only on the integral and 
so will essentially remain constant along a trajectory for which 


t 
r= % grad,E 


The velocity of the wave packet is 


Vik) = 5 grad.E (5.70) 
When a static electric field is applied to the crystal, the 
electrons will be accelerated and gain kinetic energy. Of course, 
total energy remains constant and stationary states do exist. 
The study of the electronic motion in a steady uniform field is 
surprisingly quite complicated. The Schrédinger equation now 
includes an extra potential 


U=eé.¥ (5.71) 


& is the electric field strength and —e the electron charge. 
Ordinary perturbation theory cannot be used to solve this 
problem as the extra potential becomes arbitrarily large for 
sufficiently large distances. 

Free electrons would be accelerated off to infinity by an 
electric field. In one-dimension, the operator which has eigen- 


values fk is just the linear momentum operator f = see (see 
(5-35))- (This is not the case in crystals.) If the wave function % 
represents a wave packet with mean momentum fk, then the 


rate of change of this momentum is 
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tk = 5 [ys (a — pat) dr (5.72) 
had 
= 7 aan dx2 + e&x 


The field @ is in the positive x-direction. Equation (5.72) be- 
comes 
hk = —e& (5.73) 


This final result is an expression of Newton’s second law — the 
applied force is proportional to the rate of change of linear 
momentum. 

Equation (5.73) is also true for electrons in crystals, but this 
does require some amplification. In crystals, k defines a Bloch 
wave, that is a function which is an eigenfunction of both the 
Hamiltonian and translation operators. But the new term e& . r 
does not have the lattice periodicity and & is no longer a con- 
stant of the motion. The effect of the field is to induce transi- 
tions from one energy state to another. For weak fields 
(<106Vm-1) band transitions are improbable and so the main 
effect is to alter k. If k is the mean wave vector of a wave packet 
constructed from a single band, the change in band energy of 
the ‘electron’ in unit time is 


grad, E .k 


The work done on the electron by the field in this time is 
—e&V(k) = —erad.E 6 


These expressions must be equal and 


hk = —eé (5-74) 
is a solution. 

For the mean k-vector of a packet to have significance it is 
necessary that the packet be spread over several crystal primi- 
tive cells in real space. Consequently the argument is only valid 
in general if the field is sufficiently slow-varying. A more formal 
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proof would take account of inter-band transitions (Zener 
effect). 

This result is not quite so obvious as the corresponding free 
electron case. Here Ak does not represent the real momentum 
but the ‘crystal momentum’, and will be considered further in 
the next section. Equation (5.74) is a modified Newton’s law 
for Bloch wave packets — the external force is equal to the rate 
of change of crystal momentum. 

In the presence of a combined electric-magnetic field the 
appropriate result is 


hk = —e[@+V A B] (5-75) 


B is the magnetic flux vector. 

For a magnetic field alone it would seem reasonable to write 

hk = —eV \ B (5-76) 

This cannot be proved in the same way, as the force is always 
perpendicular to the velocity and the work done by the field is 
zero, Even so, the formula is correct up to a point, but it does 
neglect the fact that in a magnetic field the electron describes 
closed orbits and a set of discrete energy levels is formed. 


§ 5.11 The crystal momentum 


The solution of Schrédinger’s equation for a free-particle is 
Aétk.r (A is a normalizing constant). This wave function is also 


e e h « 
an eigenfunction of the momentum operator "4 grad with 


momentum eigenvalue fk. 

By analogy the quantity 4k where k defines a Bloch wave, is 
called the crystal momentum. This is not the true momentum 
of the electron. Linear momentum is not conserved in a 
periodic field. 

The expectation value of the real momentum of a Bloch state 
is 


B= [bet? grad yy de (5.77) 
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The periodic part of the Bloch wave may be expanded in a 
Fourier series (plane wave expansion). 


tnt) = Sate W(k + Km) (5-78) 
From (5.77) and (5.78) 
B=Z|Wk+ Km) |?H(k + Km) (5-79) 


The real momentum may take the values i(k + Km) with the 
respective probabilities | W(A + Km) |?. 


Equation (5.75) illustrates one aspect of how the crystal 
momentum replaces linear momentum in crystals. Another is 
obtained by considering a collision between two Bloch electrons. 
For classical particles, linear momentum is conserved in a colli- 
sion. An analogous result holds true for electrons in crystals. 

A system of two electrons in a periodic field is described by 
a two-electron wave function. The two-electron Schrédinger 
equation must be invariant under a lattice translation and so 


{e | Ra}W(ri, re) = et-Raf(r, rz) (5.80) 


The wave function may be characterized by the vector k and 
hk is the crystal momentum for the two electrons. When the - 
two electrons are sufficiently far apart for their interaction to be 
negligible, the wave function is the product of the one-electron 
Bloch waves with wave vectors k; and ke respectively. 


Pri, re) = $x, (11) bx, (2) 


Clearly 
{e | Radda Pe, = 61m thy dy, (5.81) 
Comparing (5.80) and (5.81) 
k=h-+hke+ K (5.82) 


The crystal momentum of the system is equal to the sum of the 
momentum of the electrons, apart from an arbitrary reciprocal 
lattice vector. As & is a constant of the motion it follows that if 
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the two electrons collide, the sum of their crystal momenta 
before is equal to the sum of the momenta after, apart from a 
reciprocal lattice vector. 


i.e. kitke=hki’+ ke’ + K (5.83) 


§ 5.12 The electric current 


Each Bloch wave represents a moving charge with a mean 
velocity V(k) (5.67). The total current carried by a full band is 
clearly proportional to the integral 


— [ _S (RV (k) ak (5-84) 


where {(k) is the probability that the state k is occupied. In 
equilibrium this is the Fermi—Dirac function. For a full band 
S(k) =1 for all k. (Each k-state contains two electrons.) As 
energy bands always have inversion symmetry V(k) = 
—V(—k) and the current integral (5.84) is zero. 

To produce a current, an electric field must excite some 
electrons into states representing a net current. If there is an 
energy gap between a filled valence band and an empty ‘con- 
duction’ band, a finite excitation energy is required to lift the 
electron into the next band. Normally this cannot be supplied 
by a constant electric field and the solid will be an insulator. 
If the energy gap is small then at temperatures above absolute 
zero there will be a small but not zero number of electrons 
excited in to the upper band and a corresponding number of 
‘holes’ in the lower band. These electrons will be able to carry 
a current and the conductivity will increase with temperature. 
Such substances are semiconductors. If the band is not full then 
although the true ground state will represent a state with zero 
current, electrons are easily excited into the infinitesimally 
higher energy states. The conductivity is high and may be 
expected to decrease with temperature as the scattering pro- 
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Fic. 5.14 A band has inversion symmetry 


cesses become more effective. These substances are conductors. 

A Brillouin zone contains exactly as many allowed k-states 
as there are primitive cells in the micro-crystal. Each state may 
contain up to two electrons of opposite spin and if the crystal 
contains N3 cells then each energy band may contain up to 
2N3 electrons. 

Typical metals are the alkalis Li, Na, K, Rb, Cs and the 
noble metals Cu, Ag, Au. Each primitive cell contains one 
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atom and hence contributes one free-electron. There are N3 
electrons for 2N3 places and so all these metals have a half- 
filled energy band. If the solid has an odd number of available 
electrons per cell it is always a metal. Insulators may arise when 


TABLE 5.1 
Metal Semiconductor _ Insulator 
Substance Copper Silicon Diamond 
Resistivity 
Ohm cm 3°6 . 1076 ~2. 108 1022 ~ 1018 


there is an even number of electrons per cell (e.g. diamond has 
eight). This is not necessarily so as bands may overlap in energy 
and parts of a lower band may be unfilled and some electrons 


Fic. 5.15 Movement of Fermi surface in a field 
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may be in the upper band. In fact the divalent elements are 
metals because of this reason. 

The movement of a Fermi surface in a conductor is worthy of 
mention. In a zero field, the electrons are in an equilibrium 
distribution filling k-space up to a Fermi energy or if the 
temperature is above absolute zero, the occupation will de- 
crease rapidly near the Fermi energy. Immediately the electric 
field is applied the electron states begin to move in k-space with 
a velocity given by (5.74). The electron distribution will move 
bodily in the direction opposed to the field and in a time 6¢ the 
centre will move a distance —e& 5t/h. Some of the electrons in- 
crease their band energies E(k), others show a decrease. 

In a perfect crystal this process would continue indefinitely. 
A given state k would move towards the edge of the zone and 
then re-appear at the opposite side (C, C’ represent the same 
state, Fig. 5.14), and the state would continuously traverse the 
band in a direction opposed to the field. In real space the 
electron would oscillate both with and against the field (the 
velocity is zero at A, B, Fig. 5.14), undergoing a Bragg reflec- 
tion at the zone edge. The time for one complete cycle is the 
time T for k to change from —X/2 to K/2. If K is taken to be 
2a7/a where a is the separation between atom planes per- 
pendicular to the field, then 


T = h/e€a (5.85) 


For normal fields & ~ 10°Vm-1, T ~ 107° sec. 

If this oscillation took place in conductors then no net 
current would be carried. 

As conductors do allow current flow, then clearly there must 
be some other factor not yet taken into account. The process 
ignored up to this point has been the electron scattering due to 
the deviations of the crystal from true periodicity. In a real 
' crystal, electrons in the ‘high energy’ region A (Fig. 5.15), 
undergo transitions to B. These transitions, which normally do 
not involve appreciable loss in electron energy, are brought 
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about by the thermal vibrations of the lattice and by im- 
purities. A state of dynamic equilibrium is set up and if 7 is the 
average lifetime in the region A above the Ferrni-energy before 
scattering takes place, the centre of the electron distribution 
will move a distance e&7/f in k-space. The mean shift in the dis- 
tribution represents a drift velocity superimposed upon the 
random electron motion, i.e. a current. 

In metals + ~ 10-12 secs (at room temperature) and so the 
oscillations described previously do not take place. On the other 
hand, scattering does not take place in insulators as there are no 
final empty states, and so the electrons may be considered to 
carry out the oscillations. | 

This periodic wave-packet motion is quantized and the 
energy levels can be shown to be e&a apart. This system of dis- 
crete levels forms a ‘Stark ladder’. For normal fields the 
splitting is extremely small ~10~teV but it has been observed 
in p-n junctions where the local field may be ~10!°Vm-1. For 
further information the reader is referred to the literature 
(Wannier, 1962). 


§5.13 The Zener effect 


The most important effect of a weak uniform electric field is to 
accelerate the electrons and their motions in k-space are given 
by equation (5.74). In this simplified theory the electron oscil- 
lates back and forth in a bounded state. (In a perfect crystal.) 
But bounded states do not really exist in an electric field. Con- 
sider the potential energy diagram (Fig. 5.16), which shows the 
effect of an electric field on a one-dimensional potential well. 
When there is no applied field, the state with energy E is 
bound, but when the field @ is applied as shown there is a finite 
possibility for the electron to tunnel to the left. For weak fields, 
the likelihood of tunnelling is very small and so the lifetime of 
the bound state is correspondingly large. 

Tunnelling is the basis of the Zener effect. The band structure 
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in the presence of an electric field is illustrated in Fig. 5.17. 
The potential energy at the base of the conduction band has a 
slope e&. The relative positions of the other band edges are not 
altered by the field. As & varies, an electron in the valence band 
oscillates between A and B and the total energy remains con- 


é 


Conduction Valence 


Fic. 5.17 Band structure in an applied field 


stant. Now, however, there is a small but definite possibility that 
the electron will tunnel the forbidden gap to C in the conduction 
band. No change of total energy is involved. 

It can be shown that for the probability for such a transition 
to be appreciable, the field must be of the order of 109Vm-1, 
This is one theory of breakdown in insulators, but it is probable 
that other causes are more important. 

A true Zener effect has been observed at p—n junctions in 
semiconductors. 
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§ 5.14 The effective mass-holes 


In many respects the motion of an electron wave packet in a 
crystal is related to that of a free-electron. In particular, the 
velocity remains a constant when there is no external applied 
force. Near the top or bottom of a band it is often possible to 
extend this correspondence by introducing the concept of 
‘effective mass’. This is particularly useful when studying the 
properties of the charge carriers in semiconductors, 

In the vicinity of an energy minimum or maximum, if there 
are no degeneracies, the function E(k) is analytic and may be 
expanded in a power series with no linear terms. Let ko be a 
stationary point. 


E(k) = Elko) + 42 Sau be — hoi) (ks — hoy) +. 


By a suitable rotation of the co-ordinates the cross-terms may be 
removed. 


E(k’) = E(ko) + — 


— 


h2(ky? — ke’? — kg’? 
2\m* ' mo¥* * m3* 


\+. .. (5.86) 


The orthogonal co-ordinates k,;’ are linear combinations of the 
co-ordinates (k; — koi) and 


o2k 
m* = ml(se7), (5.87) 


Compare equation (5.61) for the simple cubic structure. For 
small k 


E(k) = (Ea + « + 68) — Ba2(kyz2 + ky? + kz?) 
and 
m* — —h?/2Ba? 


If there is a minimum at ko the constants m;* are all positive. 
Equation (5.70) is still valid in the new co-ordinates and 


i(k’) = (5.88) 
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The acceleration of the state in a field is from (5.74) 
m,* Vi(k’) = —ed; i= I, 2,3 (5.89) 


This final equation is analogous to Newton’s law. The motion of 
the state is classical except that three different masses m:* are 
used for the three directions. This effective mass description is 
well suited for discussing the properties of those excited elec- 
trons in a semiconductor which are near the energy minimum 
in the conduction band. 

At the top of a band, equation (5.86) will again be correct, 
but in this case the three masses m;* must be negative. It is often 
convenient to redefine the masses to be positive and then 


F h2(ky'? ke’? kg’? 
E(k’) = Elko) — (5 + Fa + Fa) 6-90) 


ok 
mt = — te/(sr), 


The equations analogous to (5.89) are 
mit Vilk’) = e& (5-91) © 


A state near the top of a band behaves as a particle with either 
positive masses m;* and positive charge or negative masses and 
negative charge. 

The former interpretation is usual for semiconductors. At the 
top of the valence band there are unoccupied electron states or 
holes. Instead of studying the nearly full band in terms of the 
many occupied states it is simpler to give a description in terms 
of the holes. A little thought should convince the reader that a 
hole will travel through space at the same rate as its neighbours 
and so the velocity of a hole is the same as that of an electron in 
that state. A full band carries no current at all and so a band 
with one electron missing must carry the negative current which 
that electron would carry. Clearly if the hole description is to be 
used then a hole must produce the same current as a positive 
electron in the same state. So, in semiconductor theory, a 
nearly full valence band is conveniently described by a few 
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positive charge carriers with positive masses and the current 
carried by the band is proportional to the integral 


f ef _ Lt —F(R)IV (A) dk 


E (k) ; 
Conduction 
band 
Positive 
curvature 


Negative 
curvature Valence 


band 


O k 
Fic. 5.18 


§ 5.15 Cyclotron resonance 


The experimental study of the motion of crystal electrons in a 
magnetic field has proved to be a most fruitful source of in- 
formation about the band structure. 

When a steady uniform magnetic field is applied to the 
crystal an electron state moves in k-space so that the vector k 
is perpendicular to both the applied field and the velocity. 


th = —eV AB (5.92) 
A wave packet in k-space does not gain kinetic energy from a 
magnetic field and so the ‘orbits’ are the intersections of the 
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constant energy surface in k-space with planes perpendicular 
to the field. The component of k& parallel to B remains constant 
at its initial value. In real space the electrons move in helices. 
(In fact the orbit will only be a true helix if the energy surfaces 
are axially symmetric about the field B.) 

The frequency of rotation of the electrons depends upon the 
magnitude of the field and the structure of the energy bands. 
Consider two neighbouring energy contours in k-space in a 
plane perpendicular to the field. 


Fic. 5.19 Magnetic orbits 


The distance between these two orbits is dE/| (grad,£) , | or 
dE/AV,. Vis the component of the velocity perpendicular to 
the field. If dk is an increment of length measured along an 
orbit then the area of the annulus between the orbits is 
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dé ¢ dk 
dé = —- } v, (5-93) 
oe Oe 
dE AI Vy, 
The periodic time of the orbit is 
T = 2 rap a ~ (5-94) 
From (5.93) and (5.94) the cyclotron frequency is seen to be 
27eB [| (dA 
o= | (ae) (5-95) 
If the constant energy surfaces are spheres 
h2k2 
B= ont 


(For free electrons substitute m for m*.) The orbits are now 


Fic. 5-20 Magnetic orbits for spherical band 


circles and if k is the component of the wave number in the 
plane normal to the field 


A = tk? = w(2m*) E/h? + constant 


“ = mam* /h2 
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The cyclotron frequency is the same for all orbits. 
w = eB/m* (5-96) 


In general, however, w is different for each orbit. A ‘cyclotron 
mass’ m,* for the orbit is defined by the equation 


h? (dA 
me* = =(aa) (5-97) 
so that 
eB 
o= Mme* (5-98) 


The effects of scattering have been ignored in this discussion. 
If B ~o-1Wbm- a free electron describes a circle in ~10710 
sec. At room temperature, many collisions per cycle may be 
expected. However, with pure materials at low temperatures 
the collision times may be increased so that several circuits may 
be made with a suitably large field. Experimentally, the fre- 
quency is found by microwave absorption. In semiconductors 
where the charge carriers occupy a small region of a band, the 
energy surfaces are usually ellipsoidal and it can be shown that 
all the orbits have the same frequency for a given field direction. 
The observed resonance line gives information about the energy 
surface. 

The application to metals is more complicated and the 
reader is referred to the literature for further information 
(Pippard, 1960). 

However, it should be pointed out that the above ‘wave- 
packet’ treatment is semi-classical and is not the full story. 
Elementary quantum ideas suggest that as the motion is 
periodic perpendicular to B then quantization may be expected 
to occur. An exact treatment for crystals presents some diffi- 
culties but the free-electron case has been examined rigorously 
and it has been shown that only certain orbits in k-space are 
allowed and that the energy difference between levels is fw. 
These are the ‘Landau’ levels (Landau, 1930). 
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§ 5.16 The effective mass approximation 


Real crystals contain many imperfections, such as dislocations, 
vacancies, impurities, etc. Some crystal properties are little 
effected by these, others are profoundly altered. The cohesive 
energy and specific heat are among the former, whereas. the 
conductivity of metals at low temperatures and the transport 
properties of semiconductors are examples of the latter. In 
general these imperfections may be regarded as localized per- 
turbations of the periodic potential. 

In Section 5.14, it was shown how close to an energy maxi- 
mum or minimum a wave packet moves as a classical particle 
but with different effective masses along the three principal 
axes. Intuitively, in the vicinity of an impurity, it would seem 
to be a reasonable guess to solve the Schrédinger equation: 


_h 02 I 
cla a2 + ® mt mi ot oe ms3* es $+ Vib = Eb (5.99) 


If V is the localized potential of an impurity atom, then E may 
be taken to represent the energy values near the impurity. The 
crystal periodic potential is ‘contained’ in the effective mass. 
This equation is correct when the right interpretation is given, 
but it certainly requires some justification. 

If #o is the perfect crystal Hamiltonian the equation to be 
solved is 


(Ho + Vb = Eb (5.100) 


The perturbed wave function % may be expressed as a sum over 
all the Bloch functions in all bands. If the perturbation is 
slowly varying it should be a good approximation to restrict the 
sum to the Bloch waves from a single band. As V is a localized 
potential it would seem sensible to express % in terms of 
localized functions. The set of Wannier functions for a band 
span the same space as the Bloch waves and so in this approxi- 
mation 


= pE (Radar — R,) (5.101) 
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If equation (5.101) is substituted into (5.100) and the result 
multiplied on the left by a*(r — R,,) and integrated, then 


> F (Rn) a*(y — Rn)[#0 + Vja(r — Rn) dt = EV(R,,) 
" (5.102) 
Using (5.66) these equations may be written 
ple(Rm — Rn) + V(Rn; Rn)|¥( Rn) = E¥(Ry) (5-103) 
V(Rm Rn) = fore — Rn)V(a — Rn) dr 


The difference equations (5.103) may be simplified by usefully 
employing the localization of the Wannier functions, When V is 
slowly varying over several cells the matrix element V(R,,, Rn) 
is nearly zero when Rm # Rn. All off-diagonal terms are 
ignored and the remaining term is, to a good approximation, 
simply V(Rm). 

The set of difference equations (5.103) describes the im- 
purity problem, but is very difficult to solve. Wannier showed 
how this could often be reduced to a differential equation. The 
function ‘'(.Rn) is defined by (5.101) only at lattice points. Let 
Y(r) be a continuous function which takes the values required 
by (5.101) at the lattice points. By Taylor’s expansion 


Y(r — Rj) = V(r) — Ry. grad ¥(r) ... (5.104) 
= exp (~R;.grad)V(r) 
The first term in (5.103) may be rearranged: 
2e(Rm — Ra) ¥(Rn) = p(y) F(R — Rj) (5.105) 


An operator may be defined from E(k) by pp ene k by 
grad/z in (5.65). 


1(E3) = Se(R) exp [—Ry. grad] (5,106) 
5 
and so the first term of (5.103) is 


£(EFS)¥ (Rn) (5.107) 
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The difference equations may be replaced by the differential 
equation 


[=(=*) as rn) eo —E¥(r) (5.108) 


The wave function ¥ is an ‘envelope’ of Wannier functions. 
The differential equation is correct to the same approximation 
as the difference equation, namely that the potential is slowly- 
varying and that the wave function can be expressed in terms 
of a single band. 

When the band energy can be a as 


h? 
CO ae ee 


the equation to be solved is 


_f 1 @ 1 0 
ear an mi tat t aes aml’ 4+ VY = (E— B)¥ 
(5-109) 


The equation is identical with (5.99), the energy of the im- 
purity state being measured from Ep. 

If an atom, such as arsenic (valence 5), replaces one of the 
atoms in the semiconductor germanium (valence 4), there will 
be one valence electron left over. There will be an excess posi- 
tive charge at the impurity and one ‘free’ electron to go into the 
conduction band. This type of impurity is a donor. 

The perturbing potential is Coulombic 


V = —e?/4reeor 


¢ is the dielectric constant of the medium. If this potential is 
considered slowly varying and if a scalar effective mass m* may — 
be used to describe the bottom of the conduction band, then 
equation (5.109) is of the hydrogen atom type. 


Y= (FE — Eo)¥ (5,110) 
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The energy levels are given by the Rydberg formula. 


—e4m* 


sie a at Bh2e2%eq2n?2 


There are discrete bound states at this depth below the bottom 
of the conduction band. For germanium a sensible value of m* 
is ~o-1m and « ~ 15. The corresponding ionization energy of 
the impurity is ~o-005 eV. (This should be compared with the 
hydrogen atomic ionization energy of 13-6 eV with m* = mand 
e = 1.) These results are in fair agreement with the donor 
levels in semiconductors. 

If there are trivalent impurities such as boron or aluminium 
in a quadrivalent semiconductor they will act as acceptors as 
they take up electrons from the valence band leaving holes. 
The impurity levels lie immediately above the valence band. 


Donor level 


Gap 


Acceptor level 
— Valence 
PITTED ETT TLE TEL OT ELT TT band 
Fic. 5.21 


In the absence of impurities the number of electrons in the con- 
duction band is equal to the number of holes in the valence 
band. If donor atoms are present in any number the electrons 
will be predominant and the material is said to be n-type. On 
the other hand, if acceptor atoms are present, then the positive 
carriers will control the conductivity and the material is p-type. 

The extrinsic conductivity of an impure specimen is in general 
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considerably higher than the ideal intrinsic conductivity of the 
pure material as the acceptor and donor states are ionized much 
more easily. 

The above treatment of donors and acceptors has been 
greatly simplified to bring out the principles of the problem 
rather than the complications. Considerably better results are 
obtained if a non-scalar mass is used (equation (5.109)), but 
then the Schrédinger equation cannot be solved analytically. 
(For germanium the absolute minimum of the conduction bands 
is not at the origin but occurs at the centres of the hexagonal 
faces of the zone. It is then necessary to define two effective 
masses. ) 

An F-centre in an ionic crystal is a similar problem. This may 
be regarded as an electron bound to a negative ion vacancy. 

An electron and a hole may form a bound pair and so generate 
a set of hydrogen type levels. Frenkel called such a pair an 
exciton. It can be shown that they would produce a set of levels 
just below the conduction band in the forbidden gap. These 
levels have been observed in highly purified germanium. 


Problems for Chapter 5 


1. The reciprocal lattice primitive cell has a volume 
Vo = b1. (be A bs). Show that Vp, = (27)3/V_ where Va is the 
volume of the direct lattice primitive cell. 

2. The reciprocal basis vectors are given by equations (5.2). 
Invert these relations to give 


ar = F(b2 A bs), a2 = F(bs A bi), as = F(bi \ bs) 


(This result follows in fact from the symmetry of the definition 


(5.1).) 


3. The two vectors (ja _ 703) and (ja —_ #2) may be 
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used to define the set of planes with Miller indices (Akl). Find 
the scalar product of the reciprocal lattice vector K = hh, 
+ kbe + lbs with each of these vectors. Hence show that K is 
normal to the (hkl) planes of the real lattice. 

4. Given that the function f(r) has the periodicity of the 
lattice verify that 


[Foe dr =o 


unless k is a reciprocal lattice vector. (Note equation (5.7).) 

5. The potential V(r) of an electron and the probability 
density are related by Poisson’s equation V2V = ep/eo. Both V 
and p are periodic functions and may be expressed as Fourier 
series. 


V(r) = 2M Ky) ¢—tKn.7, p(r) — dp( Kn)e~tke- 
* Kn 
Show that 


V(Kn) = — ap ( Kn) K, #0 


6. A cubic unit cell containing four lattice points may. be 
chosen for a face-centred cubic lattice. If this cube has a side of 
length 2a a rhombohedral primitive cell may be defined by the 
hasis vectors 


a = a(0, I, 1), dg = a(I, 0, 1), ag = a(1, 1, 0) 
(Fig. 3.4) 

Using equations (5.2) find a set of basis vectors for the reciprocal 
lattice. Verify that this lattice is body-centred cubic. (The 
basis vectors for a b.c.c. lattice may be chosen so that 
a, =a(—I, I, 1), a2 =a(1, —1I, 1), aj = a(1, 1, —1). The 
cube of side 2a is twice the size of the primitive cell.) 

7. A simple cubic lattice with parameter a has a simple cubic 
reciprocal lattice with minimum distance 27/a. The Brillouin 
zone is a cube. Show that: 


(a) the small group at k = 0 is the full cubic group m3m; 
(6) the small group at a general point A on the k;-axis is 
the tetragonal group 4mm containing eight elements; 
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(c) the small group at a general point A on a cube diagonal 
is the group 3m containing six elements. 


Find the star of the k-vector for each of these cases. Show that 
the point Y where the 4,-axis meets the zone face and the point 
R where the diagonal meets the zone corner have the small 
groups 4/mmm (with 16 elements) and m3m (48 elements) re- 
spectively. 

8. A normalized hydrogen type (15) wave function may be 


written 
a3\4 
— — —ar 
(15) = (=) e 


The normalized plane wave e*k-r/(Q) can be made orthogonal 
to the (15) function by the Schmidt process. Q is the volume of 
the crystal. Show that the resulting orthogonalized plane wave 
when k = 0 is 


to.r.w. = vo (1 — 8e-a") 


Note that this new function has a single radical node like a (25) 
wave function. Investigate the behaviour of this function when 
r>> 1/a, 

(Herring (1940) introduced the O.P.W. method. In crystals 
the plane waves are orthogonalized to all the core Bloch waves 
with the same k-value.) 3 

g. Consider the hypothetical problem of a one-dimensional 
lattice of like atoms. If there are two electron states close to- 
gether in the free atom and if these states have s- and p-type 
symmetry respectively, then on crystal formation each level 
will form a band. In the tight-binding approximation the wave 
functions are 

Ye = Artix + Aawber 


ux and 2, are Bloch waves composed of s and p atomic func- 
tions respectively. It may be assumed that all the atomic 
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functions are mutually orthogonal. Show that the two energy 
bands are given by 


H, H. Ay — Hee)? Hy, |?]# 
E(k) = 11 + Hee) + [( 7 22)? + 4 | Mie |?] 


where Ay = J Pint yk dr 


By using the substitution tan? 6 = 4| Ais |2/(Hi — Hee)? 
show that the lowest band energy function is 


Ex(k) = $[Aii(1 + sec 0) + AHee(1 — sec 6)] 
with Bloch waves 


V1 = cos oak + isin wk 


Find the corresponding solutions for the upper band. 
10. In Problem 9, investigate the two cases: 


(a) Hi # Hee for all k and 
(6) Air = Ae for some particular k = k’ 


Show that in case (a) the band edge functions have the same 
symmetry as k =o whereas in case (b) the wave function 
symmetry changes from s-type (or p-type) at k = 0 to p-type 
(s-type) at the band edges. Note that the bands themselves do 
not cross. _ 

11. Show that the periodic part of a Bloch wave satisfies the 
equation 


. h2 | heh? 
[~ —itk, grad |u,. = | Ea(h) _= en 
n denotes the band. For sufficiently small k the term 
. he h 
—i me . grad un,x = = . Pun, k 
may be treated as a perturbation and the energy found in terms 


of the solutions at k = o. Obtain the expression for the energy 
(correct to second-order), 
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a A? | <nik.p|n'>(|2 
HONE) og PMO) Ee Eto) ey 


‘én 


<n|k.p|n'> = [uno* (ke. p)uno de 


Verify that the effective mass tensor is given by 


m\ _ 2 ,<ai| pila’ ><2'| pla> 
(a), = 0 tae En{o) — En’(o) 


(The k . p approximation can be extended to include a general 
point in the Brillouin zone. When energy bands are degenerate 
the expression is more complicated.) 

12. Show that the momentum eigenfunction (5-78) and the 
Wannier function are Fourier transforms of one another. 


Wk’) = “ [ater de 
space 

13. Prove that for a symmorphic space group the Wannier 
function may always be chosen to be a basis for a one-dimen- 
sional irreducible representation of the crystal point group if 
the band is non-degenerate. 

(Hint: Equation (5.29) may be assumed.) 

14. Consider the linear lattice with spacing a with N cells in 
the repeating micro-crystal. If the Bloch functions are taken to 
be normalized plane waves show that the Wannier function at 
the origin is 

a= V4 sin ax/a 
7 x 

(Hint: Convert the sum over the Bloch waves into an integral 
over the Brillouin zone.) 

15. If #(r) is an eigenfunction with energy E for an electron 
in a steady uniform electric field & show that %(r — Rn) is also 
an eigenfunction with energy (E ++ e& . Rn). (This illustrates 
the ‘Stark Ladder system’.) 

16. Ifa magnetic field Bis applied along the z-axis the motion 
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of an electron (charge —e) is periodic in the plane z = con- 
stant. By the Bohr-Sommerfeld method the phase integral 


bp dq = (n+ ¢)h 


nis an interger, ¢ a phase factor. The momentum p conjugate 
to the position vector r is fk — eA . (B=curl A). Show that 
the permitted orbits have projected areas (n + ¢)h/eB in real 
space corresponding to 


(n + $)27eB/h in k-space 


(Hint: integrate equation (5.76) to find & in terms of r. This 
is a simplified explanation of Landau levels. See also Problem 
19.) 

17. Show that near the bottom of an s-band in a simple cubic 
structure the density of states is 


Q 
N(E) = qe [Bl [E ~ E(o)}* 

where £(o0) is the energy at k = o. (Note equation (5.61).) 

18. In general the Fermi level Er in a semiconductor 
lies in the gap between the valence and conduction bands. 
' At moderate temperatures it may be assumed that 
(E — Er) >> kT for energies in the conduction band. Show 
that the electrons obey the distribution 

Se ~ ¢@(Ey—E)/kT 
Verify that if (Er — E) >> kT the holes in the valence band 
obey 
Sa= 1 — fe & eB-Epikt 


If scalar effective masses me, mn can be used to describe the 
conduction and valence bands respectively prove that the pro- 
duct of the concentrations of electrons and holes is a constant at 
a given temperature independent of the position of the Fermi 
level. (As long as the stated assumptions are satisfied.) 


ney = 4( 2k T/h?)3(memp)2e-Z9/kT 


E, is the energy gap between bands. 
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(Hint: Take the zero of energy as the top of the valence band 
and write 


roe] 


ne = 5] Ne(E)fe AE, m = 2 [ No(E) fa dE 
‘Eg —00 


where N.(£), Ny(E) are the density of states in the conduction 
and valence bands respectively. The factor 2 arises by con- 
sidering spin.) 

19. For a spherical energy band, the effective mass equation 
(5-109) for an electron in a magnetic field may be taken to be 


h 2 
saw? grad + cA) ¥(r) = E¥(r) 


(The energy zero is chosen so that Ey = 0.) 

Ifthe magnetic field is in the z-direction it is possible to choose 
the vector potential A = (—By, 0, 0). Separate the variable » 
by using the substitution 


¥ir) — et(ksx+ kz) h( y) 


Prove that the resultant equation in J is that of a simple 
harmonic oscillator of frequency w = ¢B/m* with equilibrium 
point at yo = fkz/eB. Hence derive the following expression for 
the energy 
E= ey 
2m 


+ (n+ 4)fw 1 positive integer 


This result shows that the energy is made up of the kinetic 
energy parallel to the field and the quantized energy of the 
motion in the plane perpendicular to the field. These discrete 
levels are the Landau levels. 
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Electron Transport Phenomena 


§ 6.1 Introduction 


This chapter deals with conduction processes in metals. The 
much smaller conductivities of insulators depends on the ionic 
motion and is a completely different problem. Semiconductors 
are not dealt with specifically although it is not difficult to 
extend the treatment given below to cover most of them. The 
essential difference is that in metals the electrons are highly 
degenerate whereas classical statistics frequently apply to the 
carriers (holes and electrons) in semiconductors. Only the 
effects of steady fields and of thermal gradients will be discussed 
and neither photoconductivity nor superconductivity will be 
mentioned. 

In the previous chapter it was shown that in the absence of 
applied fields an electron in a stationary state in a perfect 
crystal has an average velocity | 


V(k) = 5 grad, E(k) (5.67) 


This implies that once an electric current was set up no field 
would be required to maintain it. The electric resistance would 
be zero. Clearly this is not correct for real crystals and so the 
deviations of the lattice from true periodicity must be im- 
portant in a discussion of transport effects. 

When a constant electric field is applied to a perfect crystal 
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an electron is accelerated so that its wave vector passes periodi- 
cally through the Brillouin zone in a direction opposed to the 
field. In Section 5.12 it was stated that in fact this does not 
happen in a real metal as electron scattering results in a steady 
state in which the electron distribution is displaced by a fixed 
amount in k-space. The two main causes of the scattering are: 


(i) the thermal lattice vibrations; 
(ii) the impurities and lattice imperfections. 


(The electron—electron interaction will have some effect, but 
this is generally small. However, this interaction is important in 
the theory of superconductivity.) 

A detailed study of the mechanics of the electron-lattice inter- 
action is beyond the scope of this book, but a brief discussion 
will give some insight. (Jones’ review article, 1956.) The 
thermal lattice vibrations give rise to a perturbation which in- 
duces electron transitions between different k-states. A vibra- 
tional state of the lattice is described by a set of quantum 
numbers giving the number of phonons in each state with wave 
vector q (Section 4.5). The scattering of an electron is inter- 
preted as the emission (or absorption) of a phonon by the 
electron. The electron loses (or gains) the quantum of energy 


E(ky) = E(ki) F hw —emission (6.1) 
. -+-absorption 


ky and k; are the final and initial electron states respectively. 
This energy change is very small compared with the electron 
kinetic energy and may be neglected for most purposes. The 
electron-lattice collisions are essentially elastic and scattering 
takes place from one point of an electron constant-energy sur- 
face in k-space to another. 

Equation 6.1 expresses the conservation of energy. It can also 
be shown that when a phonon is created (or destroyed) the 
‘crystal momentum’ satisfies the conservation relation, 


kp=kiFqtk (6.2) 
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where K is a reciprocal lattice vector. All the wave vectors are 
confined to the first Brillouin zone. (Reduced zone scheme. 
Note (5.82).) A transition for which K = 0 is called a normal 
process (N-process) whereas one with K #0 is termed an 


N-process 


/\ 


| Zone bounda ry 


U-process 


Fic. 6.1 Electron transition on the Fermi surface. 
Phonon destruction 


Umklapp process (U-process). If ki =F q lies inside the Brillouin 
zone the process is normal, but if the vector sum lies outside the 
zone so that a reciprocal lattice vector K must be added to 
make k; lie within then the process is ‘Umklapp’. This classi- 
fication of transitions into normal and Umklapp is somewhat 
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arbitrary and depends upon the choice of Brillouin zone. In the 
case of a spherical Fermi-surface (Fig. 6.1) it is convenient to 
choose a zone with the origin of the sphere coincident with that 
of the k-axes. 

The U-processes produce difficulties in the theory, but if they 
are neglected many possible transitions are ignored. In the 
example of the spherical energy band illustrated in Fig. 6.1 if 
U-processes are ignored then in monovalent metals any scatter- 
ing through angles greater than ~80° is not taken into account. 
(Problem 2.) 

The longitudinal lattice modes are the most important 
scattering agent, although transverse modes do contribute and 
may be important in some cases. Intuitively it would be ex- 
pected that phonons with frequencies so high that they are not 
thermally excited to any appreciable extent will not contribute 
to the scattering. This can indeed be verified and the only 
phonons which need be considered are those for which 


ho <kT (6.3) 


Examination of Fig. 6.1 will show that there is a minimum 
value of the lattice wave vector gm for a U-process to take place. 
This value is equal to the minimum distance between the Fermi- 
surfaces in neighbouring zones. For low temperatures the U- 
processes may be ‘frozen out’, but at high temperatures they 
are dominant. Any further discussion of the electron-phonon 
interaction would be out of place in an introductory text and 
the interested reader is referred to Jones’ review article (1956). 

The scattering in metals at normal temperatures is primarily 
caused by the thermal vibrations but at very low temperatures 
the ‘residual resistance’ due to impurities becomes important. 
The collisions of conduction electrons with static imperfection 
is also essentially elastic as the ionic mass is so much greater 
than that of an electron. Often the resistivity of a metal con- 
taining imperfections may be written 


p = p(T) + pe (6.4) 
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where p+ is the resistivity caused by the thermal vibrations and 
ps 1s due to imperfections and is temperature independent. 
Equation (6.4) is an expression of Matthiessen’s rule and may 
be expected to hold if the two sources of scattering are inde- 
pendent of one another. The zero-point thermal vibrations do 
not contribute to the scattering. This is easily understood when 
it is remembered that at absolute zero an electron can not 
absorb phonons as there are no phonons present. Similarly the 
degenerate electrons can not emit phonons as there are no 
lower states for the electrons to fall into. In fact the resistance 
of an ideal metal (pure with no imperfections) would tend to 
zero as [—> 0, 


§ 6.2 The transport equation 


A general theory of electron transport in crystals must not only 
include the effects of external electric fields but also those of 
magnetic fields and applied temperature gradients. The prob- 
lem is analogous to that arising in the theory of transport in the 
kinetic theory of gases and is solved in a similar manner. The 
main differences are: 


(i) the use of quantum statistics (in metals) instead of 
classical Boltzmann statistics; and 

(ii) that the scattering processes arise from electron-lattice 
rather than from electron—-electron interactions. (In kinetic 
theory the gas particle collisions produce the scattering.) 


At any instant the electron system is described by a single- 
particle distribution function f(k, r, t) defined so that 


re Sf (k, ¥, t) dk dr 


is the average number of electrons in the element dk of reci- 
procal space and dr of real space at the time ¢. It is assumed that 
the fields are sufficiently slow varying that dr may be chosen 
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small enough so that the fields and temperature are essentially 
constant within it and yet large enough that it contains many 
electrons. In equilibrium in the absence of external fields (or 
temperature gradients) f(k, r, ¢) is independent of r and t.and 
is just the Fermi distribution fo (2.39). In principle, the trans- 
port properties of a conductor are known when the distribution 
function has been determined under the influence of external 
fields (or temperature gradient). 

Assuming Liouville’s theorem is applicable, then in the 
absence of collisions the distribution function in the neighbour- 
hood of an electron remains constant as the electrons move 
under the applied fields. That is, the electrons which lie within 
the six-dimensional volume of phase space dk dr about k, r at 
time ¢ will occupy an equal volume about k + k di, r + Vdt 
at time ¢ + dé. 


i.e. of (Z) : +#.grad,ftk.grad,f=o (6.5) 


grad, fand grad, f denote gradients in real and reciprocal space 
respectively. In a real crystal the electron-lattice collisions pro- 
duce a net change in f along an electron trajectory. 


(F), = (Z) +f. grad,f+k.grad,f (6.6) 


(Strictly it is only possible to refer to trajectories for classical 
particles.) 

The 1.h.s. denotes the rate of change due to scattering. Equa- 
tion (5.75) may be used to express & in terms of the external 
fields & and B. 


(¥), = (£), + V. grad, f — 5(6 +V A B).grad,f 
(6.7) 


This is the Boltzmann equation. In the steady state (£) =0. 


The collision term is complicated and its evaluation is the 
major difficulty of the theory although it can simply be ex- 
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pressed in terms of a transition probability. As the collisions may 
be assumed to be elastic it is useful to define the probability so 
that P(k, k’) dS’ is the probability per unit time that an electron 
in the Bloch state & makes a transition to some empty state k’ in 
the element dS’ of the constant energy surface including both k 
and k’. As defined P(k, k’) depends on the density of states on 
the energy surface. The number of electrons per unit volume 
leaving the volume element d& for dS’ in time dé depends also 
on the number of electrons in the element dk and the proba- 
bility that the final state k’ is vacant and is given by 


gal, WF (k) dk [x — f(w)] ds" de (6.8) 


The total number of electrons leaving d& in unit time for all 
points on the energy surface is 


sal PU RY UL — F(R] ds" (6.9) 


The integration is over the constant energy surface. Similarly 
the gain of electrons by dk in unit time is 


dk , s ? 
ga) PUR, RP (k Ex —f()] 48" (6.10) 
Clearly the net gain per unit time is 


“(Z) zs a (P(k’, kK) — f(A] 


— Plk, k’)f(k)[t —f(k')]} dS’ (6.11) 


The transition probability between two states is the same for 
either direction. 


P(k’, k) = P(k, k’) 
(F) = [PG WIL) — FH) a8’ (6.12) 


This is precisely the result that would have been obtained if the 
exclusion principle had been ignored. 
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Further progress depends upon the precise form of the prob- 
ability function P(k, k’) and is difficult. The problem is greatly 
simplified when the collision integral (6.12) can be written in 
terms of a relaxation time. A function 7(k) can always be de- 


fined by 
(ZL) = fend Red Gary 


Here fo(k, r) corresponds to thermal equilibrium in the absence 
of applied fields and temperature gradients. As mentioned 
already, this function is the Fermi function and in the absence 
of a temperature gradient does not vary with r. In the general 
case when the temperature varies from point to point then fo 
will depend on r (2.29, 2.41). Of course 7(k) as defined by 
(6.13) is not unique but will depend on the applied fields and 
temperature gradients. However, under certain circumstances 
+(k) to a good approximation is independent of the fields and 
temperature gradients and it is in these cases that 7(x) is a true 
relaxation time. Equation (6.13) states that the ‘rate of relaxa- 
tion’ due to collision processes is proportional to the deviation 
of the function from equilibrium. In the absence of applied 
fields or temperature gradients Boltzmann’s equation is 


f(k, t) — folk) _ (4 
> t(k) - (Fes oe 


The solution is 


f=fo + (fi — foe” 


Ji is the initial value of the distribution function. This result 
shows that under the collision scattering processes alone, f 
approaches the equilibrium distribution exponentially with 
time constant 7. 

The introduction of a relaxation time is an approximation but 
for small deviations of the distribution function from equili- 
brium (i.e. small fields and temperature gradients) it would 
appear to be reasonable. 
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§ 6.3 Electrical conduction and Ohm’s law 


Suppose a uniform electric field @ is applied in the x-direction 
and that no magnetic field or temperature gradient is present. 
All points in the crystal are equivalent and so in the steady 
state, f is a constant in space and time. Assuming a relaxation 
time can be defined, Boltzmann’s equation becomes 


eee (k) = —3é . grad, f(k) (6.15) 
ie. Fk) = folk) + 46r(K)ee(k) (6.16) 


The fields that can be maintained in metals are weak and to 
first order in the field, f on the r.h.s. may be replaced by /o. 


Fk) = folk) + For(RZR(k) +... (6.17) 


The linear term is responsible for Ohm’s law and the higher 
terms produce deviations from it. The internal high fields at 
p-n junctions have enabled the effects of the higher terms to be 
observed. Using Taylor’s theorem (6,17) may be expressed as 


F (ke) = fol ke + 55% hy, he) (6.18) 


If 7 is independent of & the distribution moves bodily a distance 
e&r/h in k-space opposed to the field (Fig. 5.12). In fact it is 
only really necessary to assume that 7 is a constant in the neigh- 
bourhood of the Fermi surface as the second term in the ex- 
pansion (6.17) may be written 


cbr OE os 
h Ok, oF 
and is zero except near the Fermi surface. 
The current density 7 is 


j=— [PV HSH) ak (6.19) 
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The integral is over all partially filled bands. The first term in 
(6.17) represents equilibrium and does not contribute to the 
current. 


2 OE 3 
ja = i Verge De ax (6.20) 


The integration over k-space may be expressed as a combina- 
tion of an integration over the energy combined with an inte- 
gration over surfaces of constant energy. This enables the delta- 


ofo 


function characteristics of (— AB 


With dk = dSd#/| grad,£ | 


aE) to be utilized (2.40-2.47). 


en (a) ds (6.21) 
Je 4meh2 IF. Okz | wadE | ° 


The integral is over the Fermi surface and is correct to first 
order for degenerate electrons in metals although it will not 
apply to all semiconductors. The conductivity tensor o is 
defined by 

j = of (6.22) 


For a material with no symmetry properties o is symmetric with 
six distinct components. Many metals have cubic symmetry 
and for these o is a diagonal tensor with one component. 


Cubic 


7 | grad,£ | dS eee 


oak = onl (6.23) 
Ifr is a true relaxation time then a is independent of the applied 
field and Ohm’s law is verified for weak fields (<108Vm-1). 
Equation (6.23) shows that the conductivity depends only upon 
the properties of the electrons at the Fermi surface. The 
electrons deep within the Fermi surface are prevented by the 
exclusion principle from being scattered. If 7 is taken as con- 
stant over the Fermi surface and the band assumed to be of 
standard form so that 


E(k) = h®k?/2m* 
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then it is easily shown that 
o = e2nr(kp)/m* (6.24) 


where the electron concentration per unit volume is, for highly 
degenerate electrons 


3 kp 
n= zal folk) dk 2 (6.25) 


It is worthwhile observing that n is the total number of con- 
duction electrons per unit volume as all electrons are accelerated 
by the field, but as only those electrons at the Fermi surface are 
scattered then it is this particular relaxation time which is 
important. Also, it is interesting to note that the Fermi statistics 
have a decisive effect upon the electron specific heat (2.5) but 
not upon the electrical conductivity. In the case considered, 
where the band may be approximated by one of standard form 
with a k-independent relaxation time, the expression for the 
conductivity has the same form as that obtained by a classical 
treatment. 
A mean free path / is defined by 


l= V(kr)r(kr) (6.26) 


This simple one-band model may be expected to be satisfactory 
for the monovalent metals (alkalis and noble metals) and then 
l may be calculated from the observed value of the conductivity. 
For lithium / ~ 10-8m. A ‘two-band’ model has been used to 
explain transport phenomena in the transition metals while a 
multi-band model is probably more appropriate for semi- 
conductors. 

When a current flow is set up the electrons gain kinetic 
energy from the field and in the steady state this energy is 
irreversibly transferred to the lattice in the form of Joule heat. 
The electron-lattice collisions are not really elastic (6.1). If the 
theory above was extended to include second order terms in the 
field then Joule heating would be included, but the energy 
dissipation by the lattice must be considered if a steady state is 
to be obtained. 
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The conductivity is not a constant but decreases with rising 
temperature owing to the variation of the relaxation time. 1/7 is 
proportional to the number of phonons present and at high 
temperatures (T' > 9p) this is proportional to T. So in this case 
the resistivity 1/o is directly proportional to the absolute tem- 
perature. At low temperatures (7' << @p) a detailed argument 
shows that the resistivity is proportional to 7°. 


§ 6.4 Thermal effects 


A temperature gradient gives rise to energy transport (or heat 
flow) as well as to a thermoelectric current. Suppose that a 
uniform electric field is applied along the x-axis together with a 


ree: Ree : 
constant temperature gradient ae The distribution function 


is no longer a constant in space and in the steady state Boltz- 
mann’s equation is 
f ¢,4  _f-h 
Neos ho oe (6.27) 


If terms above first order in the field and temperature gradient 
are ignored 


fafot jorge ved (6.28) 


The variation with x arises from the dependence of the tem- 
perature on position (2.39, 2.51). 


af _ AT ae 
Ox dxol 
_ dT &[dEr _ E—Ep 
=—ad zeElar tT (6.29) 


7 OE a dT{dEp , E— Ep 
J=So+ 5 oe ae (ee + wat can g =} 
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The current density is then (6.19) (see Jones’ review) 


T 
a e\[ ef — (= es or ae Ko + ak} (6.30) 


Kn = opa[ (=z) “(kyo ak (6.31) 


This result shows that a temperature gradient acting alone will 
give rise to an electric current. That is, a thermoelectric effect 
is predicted. The electronic contribution to the energy transport 
per unit area is 


eee a) EW) V(k)f(k) d (6.32) 


(In the presence of an electric field the value of E(k) depends 
upon position. See Fig. 5.17.) 
In the case discussed here 
CO, = —|eé — (+ — ae Ki — re & (6.33) 

At first glance this may appear to be the heat flow but this is not 
strictly correct. In the absence of external work the heat 5Q 
added to a system is equal to the increase in internal energy 8U 
minus the increase in the free energy. If the number of electrons 
increases by 8N in a region of the metal where the chemical 
potential (partial potential) is ¢ then 


8Q = 6U — f8N (6.34) 
A free energy current may be defined by 
I 
pal Vs (R) dk (6.35) 


For electrons, the chemical potential is simply the Fermi energy 
Ey. This free energy current is clearly dependent on there being 
a particle current, but the heat flow is not. 

The true heat flow per unit area is 


O = a) {E(h) — EviV(k)f(k) dk (6.36) 
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The electronic thermal conductivity is defined when the 
crystal is in an open electrical circuit so that no current flows. 
The electrons flowing down the temperature gradient are 
balanced by an equal number with slightly lower energies 
flowing the other way. 


«=[-e/S] (6.37) 


Clearly in this case there is no distinction between energy 
transport and heat transport. Equations (6.30) and (6.33) may 
be combined to eliminate the field &. 


a ee fs ee 1dT 
2= Qs = (x ~ Ka) woe (6.38) 
a Ki?\ 1 
as «= (%— 3%) 7 
In metals, the delta-function property of (-#) may be 


employed to write Ky as an expansion in T' (2.47). 


» = Efa(Er) + (kT)? tne (Efa(Er)) +... 
(6.39) 
; OF \2 rk 
with a(Lr) = wR, (xe) oral ds (6.40) 


The integral is over the Fermi surface. If only the first term in 
(6.39) is used the thermal conductivity vanishes. To obtain a 
non-zero thermal conductivity it is necessary to include the 
second order terms and to first order in T 


2 
k= a Ta(Er) (6.41) 

Once again it is demonstrated that the electron transport 

phenomena depend upon the electron properties at the Fermi 


surface. 


198 


ELECTRON TRANSPORT PHENOMENA 


The electric conductivity is defined under isothermal condi- 
tions. From (6.22) and (6.30) 


o= e2Kog (6.42) 


A non-vanishing electrical conductivity is obtained from the 
first term in the expansion of Xo and to first order in T' 


o = ea(Er) (6.43) 
When «(£p) is eliminated from (6.41) and (6.43) 
k= oLT (6.44) 


where L is the Lorenz number (m?k2/3e?) = 2-45 . 1078 watt. 
ohm.deg-?. This is the well-known Wiedemann-Franz law. 
(«/oT is the Wiedemann—Franz ratio.) This law showing that 
« cof can also be derived using Boltzmann statistics al- 
though the Wiedemann-Franz ratio has a different value in 
this case. It is important to realize that equation (6.44) has been 
derived without reference to a specific band model and applies 
to all metals. For most metals the thermal conductivity is 
almost entirely due to the electrons and the observed values are 
in quite good agreement with the Lorenz number for a wide 
range of metals at room temperature. The relation (6.44) be- 
tween the electrical and thermal conductivities is a consequence 
of the ‘unique relaxation time’ assumption. Whenever a single 
relaxation time can not be defined then deviations may be 
expected. This does frequently happen. 

Ifa temperature gradient exists in a conductor in open circuit 
an electric field is produced. The relation between the field and 
a temperature gradient along the x-axis is 


dT 
aaa! ar Fi (6.45) 
pis the Thomson coefficient. 
Equation (6.30) clearly illustrates this effect and may be used 
to evaluate the Thomson coefficient. In the absence of an 
electric current 
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fe I Er dEr I Ki dT 
é.=—|-(7-Gr) + Telde 46) 


The Thomson coefficient is 


_ Ty d (Er 1 Ky 

w= ar( 7) + eR eal 
A more significant parameter is the absolute thermoelectric 
power S which may be defined by 


T 
ie 
sS= i aT (6.48) 
All the thermoelectric effects may be derived from S. 
I K. 


Equation (6.39) may be used to express K, and as in the case 
of thermal conductivity second order terms must be included 
if a non-vanishing result is to be obtained 
wkeT> 1 do(E 
~~ ge Lo(E) aE Ep (6.50) 
where 
e OE \2 = r(k 
(8) =—Sal (Se) panels 659 
o(E) represents the ‘conductivity’ of a metal with Fermi level 
at E (see (6.21)). Equation (6.50) does not apply to semi- 
conductors although it does apply to metals over a wide tempera- 
ture range. Unfortunately it is difficult to use in practice as it is 
not only necessary to know the value of o(£r) but also the 
detailed variation of o with the energy. This implies a know- 
ledge of the energy dependence of the relaxation time in the 
neighbourhood of the Fermi surface. It can be shown that when 
the scattering is produced by lattice vibrations a reasonable 
choice is t o E%/2 if the band is taken to be of standard form. 


Then 


mr2ke2T 
s= aaa) et (6.52) 
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If Ey ~5¢V and J ~g300°K (room temperature) then 
S ~—5.10-6V.deg-1. Thermoelectric powers of this order 
are found for metals at room temperature. It should be noted 
that in contrast to the electric and thermal conductivities, the 
thermoelectric power is linear in the charge (This holds true 
for semiconductors.) and so the sign depends partly upon 
whether the dominant carriers are electrons or holes. The 
alkali metals (except lithium) all have negative thermoelectric 
powers showing that the carriers are electrons as expected. It is 
not necessary to have hole-conduction to produce a positive 
thermoelectric power. If (Z 5B) is itself negative then 
o dk /z, 

electron carriers will produce a positive value for S. This is 
probably the explanation of the positive values observed for 
lithium and the noble metals. Antimony has a positive thermo- 
electric power and this is due to hole conduction. For further 
details the reader is referred to Jones’ review (1956). 


§ 6.5 Magnetic effects 


The two most important phenomena associated with a magnetic 
field are the Hall effect and the magnetoresistance effect. When 
a magnetic field is applied to a current carrying conductor per- 
pendicular to the charge flow a voltage is produced perpendi- 
cular to both the field and the current. This is the Hall effect 
and the voltage is the Hall voltage. Under the same conditions 
the electrical resistance is increased. This is the transverse 
magnetoresistance effect. (There is also a longitudinal magneto- 
resistance effect when the magnetic field is parallel to the 
current.) 

When an electric field & and a magnetic field with induction 
B are applied under isothermal conditions, the steady-state 
Boltzmann equation is 


f=fot+F@+VA B).gradf (6.53) 
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Fic. 6.2 Steady-state distribution in the presence of both 
electric and magnetic fields 


This equation can be solved, but the algebra is slightly more 
involved than the previous cases and instead a more intuitive 
argument will be given. The simplest case of a band of standard 
form is considered. An electric field alone shifts the whole 
distribution a distance e&r/h in k-space against the field. 7 is the 
relaxation time at the Fermi surface. This disturbed distribu- 
tion represents a current parallel to &. A magnetic field rotates 
the distribution about an axis parallel to B. With high fields 
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at very low temperatures several revolutions may be made 
(Cyclotron resonance. Section 5.15) but at higher temperatures 
lattice scattering results in a steady state and the displacement 
at the Fermi surface is 


8k = ZV Br (6.54) 


V , is the velocity component at the Fermi surface perpendicular 
to the field. The combined effect of the electric and magnetic 
fields is to produce a simultaneous translation and rotation of 
the spherical distribution so that the average position of the 
centre of the sphere is off the axis at o’ (Fig. 6.2). (Ziman 1963). 

The new steady state corresponds to a current at an angle 6 to 
the electric field. This is the Hall angle and 6 ~ eV, 7B/hkr. 
For a band of standard form 0 ~ wr where w = eB/m* is the 
cyclotron frequency. In practice the Hall effect is observed for 


@ B out 


X-axis 


Fic. 6.3 


a wire filament (along the x-axis say) so that the current 
direction can not change. The electrons deflected by the mag- 
netic field build up on the sides of the filament so that a trans- 
verse electric field is produced. This Hall field balances B and 
the magnitude is such that the angle between the total electric 
field and the current remains 6 without disturbing the distribu- 
tion. 
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The Hall field is 
oy = —620 8 small 
eB 
~ eet 


The Hall constant R is defined by 


_fy_ 2 ¢ 
— IB i m* 
For a band of standard form 
o= 2+ with o= e2nr/m* (6.55) 
Ox 
; R= —= 
ne 


The sign of the Hall coefficient is the same as that of the 
carriers. In materials where a single standard band is a good 
approximation a measurement of the Hall effect gives 
the number of carriers. The monovalent metals have 
R ~ —10~4cm3/coulomb in fairly good agreement with (6.55). 

The transverse magnetoresistance is measured by the relative 
change in the resistivity. As &, and j are unchanged in this 
approximation (zero-order in kT/Er) the magnetoresistance is 
zero for a band of standard form. 

A more complete treatment including terms second order in 
kT/Ep shows that there is a very small magnetoresistance pro- 
portional to B? for a simple band. Many metals show a 
magnetoresistance proportional to B? but larger by a factor of 
~10%, This may be explained by assuming either two over- 
lapping bands of standard form or a single anisotropic band, 
and in these cases there is a magnetoresistance effect to zero 
order in kT/Er. It is interesting to note that sodium has a very 
small magnetoresistance whereas the effect in lithium is quite 
large. This implies that the Fermi surface in lithium is not 
spherical, a point already suggested by its positive thermo- 
electric power. 
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Problems for Chapter 6 


1. Consider a model in which a crystal is regarded as a 
lattice surrounded by a classical electron gas with n electrons 
per unit volume. In equilibrium in the absence of external 
fields the average electron velocity V = o but when an electric 
field & is applied in the positive x-direction the electrons will be 
accelerated. If a relaxation time + may be used to characterize 
the electron-lattice collisions then it may be assumed that 


dV .m 


Show that the solution of this equation with the given initial 
condition is 
Vz = —— Vy =Vz=0 t>>Tr 
Derive an expression for the current density in terms of 
n, e, T, mand & and verify that the conductivity is 
ao == ne2r/m 


(Compare this result with (6.24).) 
2. Show that for a band of standard form the angle of 
scattering in an N-process is given by 


._ 6 
| q | = 2kp sin 5 


where q is the phonon wave vector and fr is the magnitude of 
the electron wave vector on the Fermi sphere. If the phonon 
Brillouin zone is approximated by a sphere of equal volume 
show that the maximum angle of scattering by an N-process in 
monovalent metals is 


Om = 2 sin-1(2-2/3) ~ 80° 
3. In anisotropic crystals Ohm’s law takes the form 


jt = OlmEm [= 1, 2,3 
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The conductivity is a second order tensor. Show how the 
symmetry of a cubic crystal demands that the conductivity 
takes the diagonal form 


Co .6) oO 
oO Cc oO 
oO oO co 


(Hint — Consider an electric field first along the y-axis and 
then along the negative y-axis and show that ozy = —owy, 


etc.) 
4. Starting from equation (6.19) show that in a crystal with 
no special symmetry properties the elements of the conductivity 


tensor are given by 
ae) as) 
e2 & Ok; 


cu = SR, baradgl| 


Note that the tensor is symmetric about the leading diagonal. 

5. In a non-degenerate semiconductor fo is the classical 
distribution function Ae-£/kT where A is some constant. If the 
conduction band is of standard form with effective mass m.* 
show that the electron contribution to the conductivity is | 


Ce = €2neTe/Me* 


where ne is the total number of conduction band electrons and 
the mean relaxation time is defined by 


te = | 7(E)(E — Ey)e- BI? AE / [(E— By)eesHr dE (a) 
Eg Eg 


The top of the valence band is taken as the zero of energy and 
Fy is the energy gap. 7(£) is assumed not to vary over a surface 
of constant energy. 
The holes also contribute to the current and the total con- 
ductivity is 
NeTe a) 


o = e%{ — + — = 
é& T ana 
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nn is the total number of holes and m,* the hole effective mass. 
Find the expression for 7, corresponding to (a). 


(Hint —— It may be assumed that the conduction and valence 


bands extend to plus and minus infinity respectively. Also 
[Pf de = Skt (* Eafe dE) 
0 0 


6. The motion of an electron state in A-space under the 
action of both electric and magnetic fields is given by 


k= —3(6+V A B) 


Suppose & = (@, 0, 0) and B = (0, o, B). If the band is taken 
to be of standard form show that after time ¢ the centre of the 
spherical distribution is given by 


Em* . 
Ke = nn sin wit 
é&m* 
ky = SR — cos wt) 


kz =0 


where w = eB/m*. Verify that the separation of any two states 
remains constant throughout the motion. This shows that the 
distribution behaves as a ‘hard sphere’. By considering the 
motion of two states initially at equal distances from the origin 
lying on the positive and negative x-axis respectively show 
that the angle through which the distribution is rotated is 
0 = wt. 

7. Suppose in the previous question the electric and magnetic 
fields are parallel. Find the position of the centre of the 
spherical distribution at time ¢. Explain why in this model the 
longitudinal magnetoresistance is zero. 
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Cohesion 


§ 7.1 Solid types 


The problem of solid cohesion has proved itself to be one of the 
most difficult to answer satisfactorily. Qualitatively it is cer- 
tainly possible to understand why solids occur. There are short- 
range attractive forces between atoms and if the temperature is 
sufficiently low so that these forces are more important than the 
effects of thermal agitation the atoms will condense to form a 
solid. One of the fundamental questions a complete solid state 
theory must answer is, ‘Why does an assembly of atoms 
generally form a periodic structure in the lowest energy state?’ 
To date no complete answer has been given and it is necessary 
to accept the periodicity as a fact and then attempt to calculate 
the energy of the lattice as a function of the lattice parameter. 

The cohesive energy of a solid is defined as the energy of the 
crystal in the lowest energy state at absolute zero minus the sum 
of the ground state energies of the isolated atoms and is a nega- 
tive quantity. The magnitude of the cohesive energy is called 
the binding energy (or heat of sublimation) and is equal to the 
work required to separate the crystal into its constituent atoms. 
For a macroscopic crystal for which the surface effects may be 
ignored, the cohesive energy is proportional to the crystal 
volume and may sensibly be expressed in electron volts per 
atom or per molecule. (Sometimes expressed in kcal per mole 
and 23-052 kcal per mole is equivalent to 1 eV per molecule.) 
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In principle the crystal energy may be obtained by finding 
the expectation value of the many-particle Hamiltonian 1.10 
(including both electrons and ions) regarding the ions as fixed. 
In this way the crystal energy is found as a function of the ionic 
separation. Of course, to do this accurately it is necessary to 
express the many particle wave function as a sum of deter- 
minantal wave functions. This procedure has proved most 
difficult in practice and resort has to be made to more approxi- 
mate methods. 

Solids may be roughly classified by the strength and nature 
of their cohesive forces into five major types. They are: 


(2) molecular crystals such as solid nitrogen with binding 
energies of up to o-:1 eV per molecule; 

(b) hydrogen-bonded crystals such as ice with binding 
energies of about 0-5 eV per molecule; 

(c) ionic crystals such as sodium chloride with binding 
energies of about 5-10 eV per molecule; 

(d) valence crystals of which diamond is a good example 
with binding energies of the order of 10 eV per molecule, and 

(¢) metals with binding energies of 1-5 eV per molecule. 


This classification is somewhat artificial and many solids fall 
into intermediate categories. An example is zinc sulphide which 
is both ionic and valence-like in character. The rest of this 
chapter is devoted to a semi-quantitative description of these 
five main types of solid. 


§ 7.2 Molecular crystals 


Molecular crystals are very weakly bound and both the melting 
and boiling points are low. This type of binding is one of the 
easiest to understand and is the source of the cohesion in 
crystals composed of inert molecules or atoms. Interesting 
examples are the rare gases which crystallize in a face-centred 
cubic structure. 
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When two neutral atoms or molecules are a few Angstréms 
apart so that their wave functions essentially do not overlap 
there is a weak attraction between them. This van der Waal 
attraction is caused by a mutual polarization of the atoms and 
is the main cause of binding in molecular crystals. 

If a neutral atom is placed in an electric field & the centre of 
the distribution shifts from the nucleus so producing a dipole. 
This effect is normally quite small as the applied field is small 
compared with the forces inside the atom. A detailed study 
shows that the average induced dipole moment is then pro- 
portional to &. When the field is due to a mutually induced 
dipole moment on another atom distance r from the first then 
€ o 1/r3 and the interaction energy is proportional to 1/r6. A 
more complete quantum mechanical argument shows that the 
potential energy of the two atoms is 


I 3 
Je — AG m8 A positive constant (7.1) 
1 2 


where Vj, V2 are the ionization potentials of the two atoms and 
&1, %2 are the polarizabilities. (7.1) ignores the energy of the 
repulsive forces due to the ‘overlap’ of the atoms but these are 
of very short range and consequently the van der Waal energy 
is much greater than the repulsive energy at the observed 
intermolecular distance. To find the crystal energy, (7.1) may 
be summed over all pairs of atoms but this is not difficult as the 
potential falls off rapidly. (Clearly an approximation is in- 
volved as the energy is not strictly a sum of two-body inter- 
actions.) The results obtained using the observed values of 
V;, o; and r are in quite good agreement with experiment. 


§ 7.3 Hydrogen-bonded crystals 


These crystals are weakly bound. A neutral hydrogen atom has 
only one electron and as expected it normally forms a covalent 
bond with one other atom. However, this is not always so and 
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under certain circumstances the hydrogen atom appears to be 
bonded with two other atoms. This is the hydrogen bond. The 
other atoms involved are electronegative (such as oxygen, 
fluorine, etc.) and it is believed that the hydrogen atom loses 
its electron to one of these. The proton then forms the bond 
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Vic. 7.1 Hydrogen bond, ice 


which is mainly ionic but partly van der Waal in character. A 
notable example of this type of bond occurs in ice where the 
water molecules are bound to each other by hydrogen bonds. 


§ 7.4 Ionic crystals 


Ionic crystals are strongly bound with binding energies of 
about 5-10 eV per molecule. (This is the energy to dissociate 
the lattice into positive and negative ions at infinite separation.) 
In the classical Born—Madelung theory of ionic solids it is 
assumed that the electrons are transferred from electropositive 
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atoms (e.g. Na, K, Ca) to electronegative atoms (e.g, O, F, Cl). 
In the model case each ion is taken to possess a charge which is 
a multiple of the electron charge and has a closed shell structure 
with spherical symmetry. This is probably a reasonable approxi- 
mation to the alkali halides. For example sodium chloride may 
be regarded as composed of positive sodium ions with a neon 
electron configuration together with negative chlorine ions with 
an argon electron structure. The overlap between the ions is 
small and so there is a large gap between the filled and empty 
electron bands. This explains why sodium chloride is an in- 
sulator. 

The ions are arranged in the lattice so that the Coulomb 
attractive force between unlike ions exceeds the repulsion be- 
tween like ions. It is well-known that no stationary configura- 
tion of charges can be in stable equilibrium under their own 
electrostatic forces (Earnshaw’s theorem). The additional force 
(as in the case of molecular crystals) is the quantum mechanical 
short range repulsion between the ions which opposes the over- 
lap of the electron distributions. Of course there are other 
forces as well, such as the van der Waal forces and also the 
van der Waal type forces produced by the polarization of the 
ions by the charges on the other ions. These are not two-body 
interactions and are more difficult to deal with, but fortunately 
their effects are small. In the latter case the symmetry of the 
lattice makes the field due to the neighbours cancel to a large 
extent and in particular the field vanishes at the centre of the 
ions, 

Within the two-body interaction approximation the total 
crystal potential energy is 


ag Uy (7-2) 


where Uj is the interaction energy between the 7 and j ions and 
the sum includes all pairs. The short range repulsion term is due 
to the overlap of the ion structures and may sensibly be repre- 
sented by Ae-u/* where A, p are constants and ry is the distance 
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between the ions. (Originally the interaction energy was taken 
as Ary-® with n~ 10.) 

(—1)e2 
4TEOl i 


Ui = + 


+ Aety/? (7.3) 


é is the magnitude of the charge on each ion and the positive 
and negative signs are taken for unlike and like charges re- 
spectively. If the crystal contains N molecules each with two 
ions then ignoring surface effects the electrostatic contribution 
to the total energy is obtained by summing the first term in 
(7.3) and is 
a, (7.4) 
47r€07 
r is some lattice parameter and may be taken to be the nearest 
neighbour distance. « is a positive number called the Madelung 
constant. 
r 
am 2 de (7-5) 

« does not depend upon the lattice distance but only upon the 
crystal structure. A typical value is « = 1:74.76 as for the NaCl 
structure. The summation (7.5) is only conditionally con- 
vergent and the arrangement of terms is important. This is a 
consequence of the long range characteristic of the Coulomb 
force law. The interested reader is referred to the methods of 
Evjen (1932) and Ewald (1921). 

The exponential term in (7.3) falls off rapidly and is only 
important for nearest neighbours. This non-Coulombic contribu- 
tion to the total energy is 


NM)e- (7.6) 


where M is the co-ordination number (M = 6 for NaCl). With 
this approximation the total lattice potential energy is 


U = N| — : 


2 
47r€0r 


a+ M rerie| (7-7) 
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The quantity in brackets is the cohesive energy per molecule. 
This expression for U contains three parameters 7, A and p. At 


the equilibrium separation es o. This condition may be 


or 
used to eliminate A and writing r = 79 
_ Nea p 
U= Sarl! Pe ro (7-8) 


For a given crystal « and 79 are known and so it remains only 
to find p. This is done indirectly by extrapolating the observed 
compressibility to absolute zero. The compressibility 8 is de- 
fined by 

I dP 

= Vp (7-9) 
where P and V are the pressure and crystal volume respectively. 
From the first law of thermodynamics 


8U = 8Q — PdV 
and at very low temperatures 5Q ~ O and 
sU ~ —PdV (7.10) 
The compressibility at absolute zero is then 


I d2U 
Bo Von (7.11) 


The crystal volume V is equal to the number of ions multiplied 
by the volume associated with each ion. For NaCl V = (2N)r3 
and then from (7.11) and (7.7), after some algebraic mani- 
pulation, it is seen that 


I ery To 
Bs = TBlaweyro—? a - (7-12) 


If Bo is found from experiment then the ratio p/ro can be cal- 
culated and used in (7.8) to calculate the lattice potential 
energy. This value may then be compared with the observed 
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cohesive energy and provides a check on the theory. For the 
alkali halides p/vo ~ 1/10 and in this case 

| Nera 

Ux = o (7.13) 


47r€01%0 “to 


The lattice energy is equal to about 0-9 of the Madelung 
electrostatic energy. For NaCl this is about 8-5 eV per molecule 
and is in good agreement with experiment. 

The Born—Madelung model described above gives the lattice 
energy as composed of two terms. Fréman and Léwdin (1962) 
have pointed out that this is not correct and that the total 
energy involves several terms of the same order as the Madelung 
energy. Briefly their argument is as follows. 

The virial theorem states that in a conservative system in 
equilibrium, if the potential is a homogeneous function of the 
co-ordinates of degree n then the average total potential energy 
V is 

2 
v=2r (7.14) 


where Tis the average total kinetic energy. (This theorem holds 
in both quantum and classical mechanics.) For a system where 


the inter-particle forces are Coulombic n = —1 and 
V = —2T (7.15) 
The total energy is then 
T+V=-—-T 
V 
=> (7.16) 


In a many-body quantum mechanical treatment of ionic 
crystals the potential is Coulombic and the virial theorem is 
satisfied both for the free ions at infinite separation and for the 
crystal. The cohesive energy is 


U=6T+ 8V (7.17) 


where 57 and SV are the changes in the expectation values of 
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the kinetic and potential energies as the crystal is formed. 
From the virial theorem 


bes ~sratt (7.18) 


The kinetic energy is increased by an amount equal in magni- 
tude to half the decrease in the potential energy. Fréman and 
Léwdin have explained this by assuming that the ions are 
‘excited’ in the crystal and possess a greater kinetic energy than 
the free ions. This is related to the contraction of the ionic 
electron distribution on crystal formation. As it is well-known 
that the Madelung energy is close to the cohesive energy the 
true change in the potential energy must be about twice as 
large as the Madelung energy. The extra potential energy is 
due to the ‘overlap of the ions’. (So is the kinetic energy.) The 
final result is that the cohesive energy is essentially equal to the 
Madelung energy but there are other large terms which cancel. 


§ 7-5 Valence crystals 


Valence crystals are very strongly bound with high melting 
points. Diamond is a good example. When two atoms are 
brought together their interaction brings about a deformation 
of the original atomic charge distribution. Bonding between 
similar or identical atoms may occur if negative charge is 
transferred from the atoms to the space between them. The two 
positive ‘atoms’ are held together by the negative electron 
space charge. This is the covalent bond. Of course confinement 
of electron charge in the bond leads to an increase in the kinetic 
energy which tends to reduce the binding energy. If the de- 
crease in potential energy is greater than the increase in kinetic 
energy binding will result. According to the virial theorem, 
when such a stationary state exists the increase in kinetic energy 
is equal to half the decrease in potential energy. (Note that the 
kinetic energy has a maximum value in the lowest energy state.) 
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In a crystal, the total binding energy will not only include the 
energy within a bond but also the interactions between different 
bonds. However in many cases the crystal energy is quite close 
to the sum of the bond energies and this is particularly true for 
diamond. The C—C bond has a strength of about 3-6 eV in 
diamond and in organic molecules. This suggests that a study 
of a single covalent bond would be useful. 


Hydrogen-molecular ion Het 


This is the simplest type of molecular structure and illustrates 
the essential effects of the covalent bond in that two hydrogen 
nuclei are held together by a negative charge between them. 
Within the Born—Oppenheimer approximation this is a single 
particle problem, the two nuclei being in fixed positions. If the 
protons are a distance R apart the Hamiltonian is 


2 2 2 2 
# = —Pys_ € +2) 4 —— (7-19) 


2m Ameo\T1 —r2 A4mregR 


where 71, 72 are the distances of the electron from the two 
protons. 
The electron wave function satisfies the equation 


Hs = E(R) i (7.20) 


+e +e 


Fic. 7.2 Hydrogen-molecular ion 
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and the energy of the system may be found for a given R. 
Equation (7.20) can be solved exactly in elliptic co-ordinates 
and the energy found for different values of R. 


Fic. 7.3 Minimum of E gives stable separation 


The turning point on the E(R) curve defines the equilibrium 
separation of the protons. 

Calculation predicts the separation to be 1-06 Ain agreement 
with experiment. The wave functions may be classified by their 
components of orbital angular momentum about the bond axis 
and it is found that the lowest stable state has an axially sym- 
metric o wave function with zero angular momentum. To find 
the cohesive energy of the molecule it is necessary to subtract 
from the minimum £ the energy of a proton and hydrogen 


219 


SOLID STATE PHYSICS 


atom (in its lowest state) at infinite separation. (This is 
1 Rydberg ~ 13:5 eV.) The binding energy is found to be 
almost 2-8 eV again in good agreement with experiment. 

A less exact but still useful approach is to obtain the mole- 
cular orbital as a linear combination of atomic orbitals 
(L.C.A.O.). The simplest approximation to the ground state 
would be to set the molecular orbital 


ob = Cihr + Carbe (7-21) 


where #1, y2 are the real normalized hydrogen (is) wave 
functions centred on the two protons. The symmetry group of 
the Hamiltonian is the infinite group 0/mm. This group includes 
the inversion through the centre of the bond and for a non- 
degenerate state Ce = +Ci. The two wave functions are 


thg = Cr(hi + de) gerade (even) 
by = Cy’ (da — ge) ungerade (odd) 


C; and C;’ are determined by normalizing the molecular wave 
functions. The corresponding energies are obtained from 


| b*H dr 
k= 
[ure dr 
i.e. Ey = (Hu + #22)/(1 + S$) (7.22) 


Eu = (#11 — KH x2) | (1 —_ S') 
where 


Hy = [bey dr (7-28) 
and § is the’real overlap integral 
ce | da*de dr 


S is considerably less than unity and 


Ey ~ #14 #1 (7-24) 
fy ~ Hun —- Hw 
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Bonding state 


Antibonding state 


Fic. 7.4 
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As 1 and yz are both positive in the region of overlap then it 
can be shown that #12 < o and the lowest energy state is the 
‘gerade’ state with the symmetric wave function. In general it 
is true that the more overlap there is between the atomic 
orbitals the lower the energy will be for the gerade state and the 
stronger the bond will be. The lowest energy state is often called 
the bonding state and the charge density is 


| wo |2 = Cr2(Ya? + Yo? + 2yrye) (7.25) 


The term 2i1%2 represents the bond charge. Compare this with 
the charge density from the ungerade (anti-bonding) state 
where charge is removed from the bond. 


[ au [2 = Cr’2 (yr? + 2? — 2prde) (7.26) 


With the hydrogen (15) states the binding energy is calculated 
to be ~1-8 eV which involves an error of 1 eV. A considerable 
improvement is obtained if some parameter, such as effective 
nuclear charge, is included in the wave functions and varied to 
minimize the energy. 


The hydrogen molecule He 


This is the simplest molecule involving more than one electron 
but even so no exact solution exists. The hydrogen ion bond 
wave function should be able to take up another electron with 
opposite spin and result in a stable state. The energy with the 
electron pair in the bonding orbital will be lower than that 
of two hydrogen atoms at infinite separation. Certainly the 
binding energy of the hydrogen molecule will be greater than 
that of the hydrogen ion but not twice as large because of the 
repulsion between the electrons. These predictions are con- 
firmed by experiment and the binding energy is found to be 
4°7 eV. 

The simple molecular orbital method does not give good 
agreement with experiment. In this model both electrons are 
taken to be in the hydrogen ion bonding state with opposite 
spins and the two-electron wave function is 
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oe soekh Yg(t)a{r) Pg(1)B(r) 
Hale 72) = “7a | fo(a)er(2) bo(2)B(2) (7.27) 


= Ts o(1)%o(2)[a(1)8(2) — «(2)8(1)] 


The spacial part may be expanded to give 


fo(1)g(2) = Cr? [ypi(1)Y1(2) + yo(1)p2(2) 
+ 1(1)y2(2) + e(1)Y1(2)] (7.28) 


The wave function gives equal weight to the ionic states of the 
type y1(1)y¥1(2) where the two electrons are on the same atom 
and to the remaining states in which each atom has one 
electron. Energetically it is clear that the ionic states should not 
be as important as this wave function suggests. The electrons 
tend to keep apart and this may be taken into account by con- 
figuration interaction in which the two-electron wave function 
is expressed as a linear combination of determinants. Group 
theoretical considerations indicate that the lowest deter- 
minantal. wave function which will mix in is that obtained by 
exciting both electrons to the antibonding orbital. This wave 
function is 
— 1 _|pu(t)a(t) u(t) B(1) 
Beles re) = 75 |hu(2)2(2) ul) B(2) (7-29) 


with spacial part 


Pu(t)u(2) = Cy’2[r(1)Y1(2) + yo(1)p2(2) 
— ¥1(1)f2(2) — Ye(t)¥i(1)] (7.30) 


By taking a suitable combination of (7.27) and (7.29) the ionic 
terms may be reduced or eliminated. When the energy is 
minimized in this way, the combination is such that the ionic 
terms are relatively unimportant. 7 
Clearly then, a better first wave function (instead of (7.27)) 
would be to take the electrons on different atoms to begin 
with. This is the basis of the Heitler~London approach where 
the spacial part of the wave function is approximated by 
[y(t )y2(2) + ye(1)y1(2)]. Some electron correlation is built 
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in already in this method and the energy is minimized by 
varying parameters in y and 2. This method is obviously 
correct in the limit as R——> © as the electrons are then indeed 
on different atoms. This Heitler-London wave function and the 
best two-determinant molecular orbital wave function are of 
equivalent accuracy and both involve an error of about 0°5 eV. 


Extension to crystals 


The physics of the covalent bond is illustrated by Het and He 
and the insight gained by the study of these two molecules may 
be used to tackle the much more difficult crystal problem. In 
a covalent crystal there are electron-pair bonds between the 
atoms, one electron being donated by each atom. (The Pauli 
principle does not permit a third electron in a bonding orbital.) 
Essentially the energy of a bond depends upon the overlap 
between the atomic orbitals on the different atoms. It is often 
possible to increase the overlap by mixing the atomic orbitals to 
form hybridized wave functions which are localized along the 
bond. A classic example is given by the carbon atom. 

The free carbon atom in the ground state has the electron 
configuration (15)2(25)2(2p)2. It is well known that the carbon 
atoms in molecules or crystals often have a tetrahedral environ- 
ment (e.g. methane, diamond). This implies that there are four 
unpaired electrons available for bond formation. This quad- 
rivalent property may be explained if it is assumed that in the 
‘yalence state’ one of the (25) electrons is excited to a (2p) 
state to give the configuration (15)?(25)(2p)%. From the four 
valence state electrons it is possible to construct four directed 
hybrid orbital by a unitary transformation. 

These may be chosen to have their maxima in the four direc- 
tions (1 1 1), (7 7 1), (1 7 1) and (7 1 1), Each hybrid has 
the form 


iit = $(s + +/3piin) (7-31) 
with pul = (t= + py + pz) 
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Fic. 7.5 


By analogy with the hydrogen molecule, a bond orbital y(r) 
may be constructed between two carbon atoms A and B. 


. I 
ie. x(r) = Var es) + ta) (7-32) 
where é4 is the hybrid orbital on site A pointing towards the 
nearest neighbour B and S is the overlap integral. The gain in 
energy on bond formation more than compensates for the 
energy required to excite the valence state. 

These intuitive ideas can be formalized by a band theoretic 
approach. In Chapter 5 it was shown how a non-degenerate 
band can be described by a set of orthogonal localized orbitals. 
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Each Bloch function can be expressed as a linear combination 
of these Wannier functions. There is one Wannier function per 
primitive cell. In the same way it can be shown that a set of 
degenerate bands may be described by a set of orthogonal 
localized orbitals which are identical apart from their location 
and/or orientation in space. If the bands are completely 
occupied these orbitals are called equivalent orbitals, otherwise 
they are referred to as site orbitals. For diamond there are four 
filled valence bands corresponding to the atomic states 
(2s)(2p)3. These bands are degenerate at certain points in the 
Brillouin zone and it can be shown that the equivalent orbitals 
may be chosen to be centred on the bonds between the carbon 
atoms, there being four in each primitive cell (Hall, 1952, 1958). 

In a first approximation the many-electron wave function is 
written either as a determinant whose elements include all the 
occupied Bloch functions or whose elements are the equivalent 
orbitals. These are identical, but it is often preferable to work in 
terms of localized functions. Configuration interaction may be 
conveniently introduced by mixing in other determinantal 
wave functions in which one or both electrons in a bond are 
excited. Work by Schmid (1953) on the cohesive energy of 
diamond suggests that a single determinant of bond equivalent 
orbitals is a reasonable approximation. The configuration 
interaction which is so important in the hydrogen molecule does 
not play such a decisive part in diamond. 


§ 7.6 Metals 


The attempts to calculate metal cohesive energies are somewhat 
different in approach than those used for covalent crystals. In 
valence crystals it has proved useful to introduce the equivalent 
orbitals, there being one equivalent orbital for each bond. This 
method is not very appropriate for metals as the number of 
equivalent orbitals is less than the number of bonds. Metallic 
sodium crystallizes in the b.c.c. lattice with eight nearest neigh- 
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bours demanding four bonds per primitive cell. Sodium has one 
‘free’ electron and the valence band is only half-filled in the 
lowest energy state. It is possible to form localized equivalent 
orbitals by taking linear combinations of the Bloch functions of 
the lower half of the band, but there is only one such orbital for 
every two primitive cells. These orbitals are not as localized as 
those obtained for covalent crystals. Clearly it is impossible to 
give these equivalent orbitals the same significance as in the 
case of covalent structures. 

In practice it has proved more fruitful to recognize the 
similarities between metal electrons and free electrons in co- 
hesive energy calculations. This is not to say that the metallic 
bond is of a completely different nature. There is certainly a 
covalent effect in that there is a movement of charge into the 
space between atoms. 


Cellular method 


The first steps in a cohesive energy calculation is to find an 
accurate solution to the single electron problem. Wigner and 
Seitz (1933, 1934) introduced the cellular method to do this 
for the lowest energy state in monovalent alkali metals. This 
method was briefly mentioned in Chapter 5 where it was 
pointed out that the approach is completely different from the 
plane wave and tight-binding method. An attempt is made to 
find as accurate a solution as possible in the region of the ion 
cores and then the lattice periodicity is imposed as a boundary 
condition. 

Real space may be divided into space-filling cells in the 
same way that the reciprocal lattice is divided into Brillouin 
zones. Each lattice point is surrounded by such a Wigner-—Seitz 
cell formed by the planes bisecting the lines joining the point to 
its neighbours. For crystals with b.c.c. direct lattices (e.g. 
sodium) the Wigner-Seitz cell is a truncated regular octahedron 
with fourteen faces. 

The potential at a point inside a cell may be divided into two 
parts, that arising from within and that arising from outside the 
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cell. Again, there are two contributions to the potential from 
within the cell. Firstly, there is the effect of the ion core. This 
has a spherically symmetrical closed shell and the potential may 
be approximated by some analytic representation of the free 
ion potential. Secondly, there is the potential arising from the 
smoothed out charge distribution of all the other valence 
electrons within the cell. This may be neglected in a first 
approximation. The potential from outside the cell is due to the 
other ion cores and valence charge not included in the cell. 
Symmetry requires that each cell be electrically neutral. It is 
consistent to assume that all cells are almost spherically sym- 
metric and so the potential outside an isolated cell is very small 
and may be neglected. The potential V inside any cell may be 
taken to be entirely due to the spherical ion core within. 

In a simple treatment the cell is replaced by a sphere of equal 
volume. The lowest state (k = 0) of the valence band has the 
same symmetry as the atomic functions from which it is derived. 
In sodium, the valence band is derived from the (35) state. The 
wave equation is separable in the co-ordinates and the radial 
part of the wave function R(r) satisfies the equation 


d ,dR 
aa dr) + le) —VIR=0 (733) 


The allowed energies E(o) are determined by the boundary 
conditions imposed which are that the wave function and its 
derivative must be continuous at the cell boundary and also 
satisfy the Bloch condition. In the case considered 


(F),=° (7-34) 


where 7 is the radius of the Wigner-Seitz sphere. In addition 
R(r) must have two spherical nodes as it must resemble the 
(3s) state near the ion. The values of E(o) corresponding to 
different values of ro may be calculated. At the observed lattice 
spacing, the energy obtained from the boundary condition 
(7.34) is considerably lower than that for a free atom and is the 
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main cause of binding in metals. Wigner called this the 
‘boundary correction’. Of course it is not possible for all the 
electrons to be simultaneously in the state k = o. If an effective 
mass m* may be assumed then the energies of the higher states 
are 


E(k) = E(0) += (7.35) 


For the alkali metals it may be shown that the electrons behave 
as if they were almost free and it is reasonable to take m* ~ m. 

When the potential arising from all the other valence elec- 
trons inside the cell is ignored then in effect some electron cor- 
relation is introduced as it implies there is only one electron 
inside the cell at any given moment. The crystal energy is 
essentially the sum ) E(k) of the energies of the electrons in the 
field of their own ion cores. All the other electron-ion, electron— 
electron and ion-~ion interactions in the Hamiltonian 1.10 
cancel in this approximation. It is not difficult to show that the 
average electron energy is 


E=E(0o) + Es (7.36) 
: _ 3 km 
with EE; = 5 om* 


where km is the wave vector for the highest occupied state. The 
binding energy is then 


EA) — [E(o) + Ee] (7-37) 
where E(A) is the free atom energy. For sodium the boundary 
correction F(A) — E(o) is about 3 eV and for a half-filled band 
the mean band energy Ex ~2 eV. The binding energy is 
~1 eV in quite good agreement with experiment. For a full 
band the band energy would just about cancel the boundary 
correction and there would be no binding. 

The above is a greatly simplified treatment of the problem. 
If the smoothed-out electron distribution inside the sphere is 
not neglected (Problem 8) then it can be shown that the energy 
of the crystal is increased and becomes greater than that of the 
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free atoms and so there is no binding. This somewhat surprising 
result is due to the fact that the correlation originally incor- 
porated by ignoring this term has been removed. However, if 
correlation is reintroduced in a more consistent manner it may 
be shown that the correlation energy more or less cancels the 
smoothed electron term, so restoring binding. For this purpose 
Wigner and Seitz considered a gas of free electrons of the same 
density as the metallic electrons but these modifications are 
beyond the scope of this book. 

In the method described above the result depends upon the 
atomic volume and not upon the crystal structure. Refinements 
have been introduced by retaining the Wigner-—Seitz poly- 
hedron and fitting a wave function to satisfy the boundary 
conditions at a selected number of points on the surface. 

Metals are plastic solids as the binding forces are essentially 
non-directional. On the other hand the directional character 
of the bonds in covalent solids resist shear and consequently 
covalent solids are brittle. 


Problems for Chapter 7 


1. Consider an infinite line of ions spaced a distance a apart. 
The ions have charges of equal magnitude but alternate in sign. 
Show that relative to the nearest neighbour distance the Made- 
lung constant has the value 2 In 2. 

2. In the text the short-range repulsion term between the 7 
and 7 ions is taken to be Ae-rwe, The early workers took this 
interaction to be Ary where A, n are constants. Show that a 
value of n ~ Lo is consistent with equation (7.13). 

3. In free space the flux of the electric field & out of a closed 
surface S is 


$é.ds=1 
Ss €0 
where q is the included charge and ¢o is the permittivity of free 
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space. The scalar potential ¢ is defined to within an arbitrary 
constant by & = —grad ¢. Show that ¢ can not have a 
maximum or minimum value except at points where there is a 
positive or negative charge respectively. Deduce that a free 
charge can not be in stable equilibrium at a point unoccupied 
by charge. This result implies that no stationary arrangement of 
charges can be in stable equilibrium under their own influence. 

4. In the discussion of the hydrogen molecular ion #y was 
defined by equation (7.23). Show that 


ps Go een mae 
12 =( H+ re 2) — | pr oe 
and it ae (En + sor) —| hod 
where Ey is the hydrogen atom eigenvalue for the (15) state. 
A similar result holds for #22. Using equation (7.24) derive the 
expression 

He Ha aR ok i oo le Arreqrt oa 
and also the expression for E,. Explain why E, < Ey, and state 
why this implies that #12 < 0. 

5. Suppose a system of classical particles is in motion. 
Newton’s second law states that X; = mi; where X; is the 
force in the x-direction acting on the 7th particle of mass mj. 
Confirm the expression 


My Ma». d A 

kt = — See? + 4 ge (marees) 
Show that provided the particle is confined to a finite region of 
space, the time average as ¢ > ©0 of the last term on the r.h.s. 
is zero. For conservative systems the forces acting on the particles 
may be derived from a scalar potential. 
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SOLID STATE PHYSICS 


If V is a homogeneous function of the co-ordinates of degree n 
show that 


=f 
2 


where V and T are the time averages of the potential and 
kinetic energies of the system as £ ——> 00. (You may use Euler’s 


OV oV OV ae ae: 
theorem ere + ria as) = nV.) This is the virial 


theorem for classical mechanics. 

6. A classical particle of mass m undergoes harmonic oscilla- 
tions with angular frequency w. Show that the potential energy 
at any point in the path is mw?x2/2. Hence show that the time 
average of the potential energy taken over one period is equal 
to the time average of the kinetic energy. Explain how this 
result is in accordance with the virial theorem. 

7. A hybrid orbital % is to be formed from the functions (2s), 


(22), (2py) and (2). 


ie. p=f(r)[a + +/3(dcos 6 + ¢ sin 0 cos d 
+ dsin @ sin ¢)] 


r, 6, ¢ are spherical polar co-ordinates, a, b, c, d are constants, 
J (r) is a radial function. The directions in which ¢ has a maxi- 


mum value may be obtained by putting = ed = 0. Show 


that the value of the orbital in the maximum direction (orbital 
strength) is 


S=f(lat V3 Fete] 


If % is normalized to unity 1 = a? + 5? + ¢2 + d?, Show that 
for a pure (2s) orbital S = f(r) and for a pure (2) orbital 
S = /3f(r). The maximum possible strength of such a hybrid 


is given by F(a) = 0. Confirm that this is obtained when 
tetrahedral orbitals are constructed (equation (7.31)). (Pauling 
(1960) has indicated that the ‘strength’ of a wave function is a 
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measure of its bond forming power. Murrell (1960) has ex- 
tended this principle to produce maximum overlap of the 
hybrids involved in a bond.) 

8. The conduction electron wave functions in an alkali metal 
are sometimes approximated by plane waves yx = ef-r/4/Q 
where Q is the metal volume. Verify that the charge density p 
is a constant. In the Wigner—Seitz method, the Coulomb 
potential at a point depends only upon the charge distribution 
within the sphere that contains the point. 


= PUP) dr 


V is the volume of the Wigner—Seitz sphere. The electronic 
charge in each cell is —e. Show that the contribution to the 
potential from the electronic charge is 


47r€0 2 Ts; 2rs? 


where r; is the radius of the sphere. Explain why this smoothed 
out charge adds a contribution 


+f (—e) dr/V 


to the metallic cohesive energy. (Note the factor 4.) Verify that 
this addition to the energy is 0-6¢2/47eors. (This extra term gives 
the energy in a ‘modified’ Hartree approximation. Correlation 
is ignored entirely.) 
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Solutions to Problems 


Chapter 1 


6. H=£ 44 (p.A+d.p)t+ia—a 
7 om © am \P° BP) om 


x = — yn tha. grad —L ikdiv A+ = a2 — of 
am m "8 2m 2m 
Chapter 3 
5. 
2mm 22 My Mz 
E 22 My Mz 
22 2z E Mz My 
My My Mz E 2z 
Mz Mz My 22 E 
6. 
3m E 32 322 my me m3 
E E 32 322 my me ms 
32 32 322 E ms my, me 
322 322 E 32 mz my my 
my my me ms E 32 323 
me me m3 my 322 E 32 
m3 ms my m2 32 322 E 
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SOLUTIONS TO PROBLEMS 


7. As in 6 with 


(123)~E (231)~3z2 (312) ~ 32 
(132)~m (321)~m, (213) ~ms 


10. Four one-dimensional irreducible representations. 
Two two-dimensional irreducible representations. 
Four three-dimensional irreducible representations. 

11. Two one-dimensional irreducible representations. 
One two-dimensional irreducible representation. 
Two three-dimensional irreducible representations. 

13. f(t) +f(—a), F(#) —f(—2). 

14. Be) xy, E) xz, yz. 

16. Aj, As. 


Chapter 4 
5. Q/ 278 


Chapter 5 
g. Ee(k) = $[Air(1 — sec 6) + Hee(1 + sec 8)] 


ey: 0 
Yo = — tsin > diz — cos = ex 


Chapter 6 


0 0 
5 t= f 7(E)(—E)s eben de] [ (—E)smebier dE 
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Abelian group, 65, 100, 129 
Absorption, infra-red, 92 
of phonons, 105, 186 
Acceleration of electrons, 156, 186 
Acceptor, 176 
Accidental degeneracy, 66 
Acoustic mode, go, 108 
Adams, W. H., to, 17 
Adiabatic approximation, 
Oppenheimer approx. 
Alkali halides, cohesion, 212-217 
Madelung constant, 214 
Alkali metals, cohesion, 226-230 
correlation energy, 230 
paramagnetism, 35-36 
specific heat, 31-33 
Anharmonic terms, 106, 115-116, 118 
Annihilation operator, 105 
Antibonding orbital, 220, 221, 222 
Anti-ferromagnetism, 33 
Antisymmetry of wave functions, gro, 
66 


see Born 


Atomic orbitals, 50, 70, 75 
linear combination of, see L.C.A.O. 
Attenuation coefficient, 83 
Axial rotation group, 67 
Axis, screw, 57, 137 


of symmetry, 43 


Band structure, r26—128, 131, 138, 144, 
150 

Basis vector, 42, 43, 106, 122, 123 

Binding energy, 209 

Bloch functions, 127, 130 

Bloch’s theorem, 124-130 

Blount, E. I., 153, 183 

Body-centred cubic Lacs, 46, 124, 
178, 227 

Boltzmann transport equation, 190, 
192, 193, 196, 201 

hm, D., 14, 17 

Born von Karman, boundary condi- 
tions, see boundary conditions 

Born—Oppenheimer approx., 4-8 

Bose-Einstein statistics, 105, 109 

Bound states, 37, 38 


Boundary conditions, 19, 37, 38 
periodic, 22, 85-86, 128 

Boundary correction, 229 

Bravais lattice, 41-46, 57, 71, 124 

Bravais net, 42, 71 

Brillouin zone, 83, 90, 112, 130—-13% 


Cellular method, 227-230 

Centre of symmetry, 43 

Character, definition, 63 
orthogonality, 74 

Chemical potential of electrons, 27, 
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Class, group, 63-65 
Classical statistics, 24, 28, 185, 206 
Cohesion, covalent, 217-226 
ionic, 212-217 
metallic, 226-230 
Collisions of electrons with phonons, 
186-187 
Conduction bands, 168, 175, 176 
Conductivity, theory, 193-196 
tensor, 194, 205, 206 
Conductor, 121, 161, 162 
Configuration interaction, 14 
Co-ordination number, 214 
Core electrons, 5 
Correlation energy, 14 
Coulomb, hole, 13 
interaction, 6, 9, II, 


216 
Coulson, C. A., 78, 120 
Covalent bonds, 217, 218, 225 
Creation operators, 105 
Critical points, 132-133 
Crystal, momentum, 158-160 
system, 45, 46, 57, 124 
Curie temperature, 33 
Cyclic boundary conditions, see boun- 
dary conditions 
Cyclotron, frequency, 171, 172 
resonance, 169-172 


12, 13, 14, 214, 


Debye, approximation, 111-115 
P., 120 
temperature, 114, 115 


237 


INDEX 


Descneency: and symmetry, 51-53, 65~- 


temperature, 28 
Density of states, 21-22, 30, 31, 132-133 
Destruction operator, 105 
Determinantal wave function, 10, 16, 
223 
Diamagnetism, 33, 34 
Diamond, 224, 226 
Dielectric constant, 175 
Direct lattice, see Bravais lattice 
Directed valence, 224-226 _ 
Distribution function, see Fermi—Dirac 
Donor, 175, 176, 177 
Drift velocity, 164 
Dulong and Petit law, 108-109, 114, 
1X5 


Earnshaw’s theorem, 213, 231 
Effective mass, 167-169, 175, 176, 181 
approximation, 173-177 
Einstein, A., 110, 120 
model, 110 
temperature, 110 
Electric current, 160 
Electron-phonon interaction, 186-187 
Empty lattice, 140, 141 
Energy bands, see band structure 
gaps, 145 
Entropy, 26, 27, 39 
Equivalence of representations, 60 
Equivalent orbitals, 226, 227 
Evjen, H. M., 214, 233 
Ewald, P. P., 214, 233 
Exchange, hole, 12 
energy, 13 
Exciton, 177 
Exclusion principle, see Pauli principle 
Extended zone scheme, 144 


Face-centred cubic lattice, 45, 46, 124, 
178 
F-centre, 177 
Fermi energy, 27, 163 
variation with temperature, 28~31 
Fermi-—Dirac, statistics, 24-28 
distribution law, 27, 28, 31, 160, 192 
Fermions, 24 
Ferimagnetism, 33-34 
Ferromagnetism, 33 
Forbidden region, 83, 127, 128 
Function space, 57-59 
invariant, 58 
irreducible, 58-59 
Free electrons, 18-36 
Free energy, 197 
Frequency distribution, 86, 111, 112, 
117 
Froman, A., 216, 217, 233 
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Glide plane, 57, 137 

Griffith, B. A., 94, 120 

Group, 
Abelian, 65, 100, 129 
and normal modes, 96-98 
class, 63-65 
definition, 53-54 
isomorphic, 73 
permutation, 72-73 
point, 57 

Group velocity, 116 


Hall, G. G., 226, 233 
Hall, angle, 203 
constant, 204 
effect, 201-204 
Hamiltonian, 4, 129 
ae crystal symmetry, 51-53, 65— 
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Harmonic approximation, 92 

Harmonic oscillator, 94, 101-105 

Hartree method, 8, 9, 13 

Hartree-Fock method, 10-14 

Heat capacity, see specific heat 

Heat of sublimation, 209 

Heine, V., 76 

Heitler~-London model, 223, 224 

Hermitian operator, 2, 15 

Herring, C., 179, 183 

Holes, 160, 167~169 

Hybrid orbital, 224, 232 

Hydrogen, molecular ion, 218-222 
molecule, 222~224 


Ice, 212 
Impurities, 173, 175 
Infra-red absorption, 92 
Insulators, 121, 162, 213 
Integrals involving Fermi—Dirac func- 
tion, 28-30 
Inversion symmetry, of crystal, 43 
of bands, 160, 161 
Ionic, solid, 212-217 
states, 223 
Irreducible representations, 60-61, 135 
Isomorphic groups, 73 


Johnston, D. F., 77 
Jones, H., 186, 208 
Joule heat, 195 


Kittel, C., 120 

Koopmans, T., 13, 17 
theorem, 13 

Kramers, H. A., 127, 184 

k-space, see reciprocal lattice 


Lagrange multipliers, 9, 26 
Lagrangian, 15 


Landau, L., 39, 120 
diamagnetism, 36 
levels, 172, 182, 183 
Langevin paramagnetism, see para- 
magnetism 
Larmor precession, 34 
Lattice, translation operator, 49, 100, 
128, 129, 130 
vibrations, 78-116 
L.C.A.O., 151, 179-180, 220 
Lenz’s law, 34 
Liftshitz, E, M., 120 
Liouville’s theorem, 190 
tan crystal, 98-106, 124-128, 1309, 
181 
Linearly independent functions, 10, 
8 


5 
Longitudinal mode, 78 
Lorenz number, rg9g 
Lowdin, P. O., 216, 217, 233 


Madelung constant, 214, 230 
Magnetic, orbits, 170, 171, 182 
susceptibility, 34, 36 
Magnetoresistance, 
longitudinal, 201, 207 
transverse, 201, 204 
Mandl, F., 1, 17 
Many-body problem, 4 
Maradudin, A. A., 120 
Mass, effective, see effective mass 
weighted co-ordinates, 93 
si representation of groups, 59- 


3 

Matthiessen’s rule, 188 
Metals, cohesion, 226-230 
Methane, 75, 224 
Molecular, crystal, 210-211 

orbital, 220, 221, 222, 223, 224 
Momentum, crystal, 158-160 
Murrell, J. N., 234. 


Nearly free-electron model, 139-146 
Nearest neighbour interaction, 79 
Normal, co-ordinates, g2—96 

mode of vibration, 94 
Normalization, 2 
N-process, 187, 205 


Ohm’s law, 193, 194 
One-electron approximation, 8-14 
Optical mode, go, 108 
Orbital, atomic, see atomic orbitals 
strength, 232-233 
Orthogonal transformation, 93 
Orthogonality relation, 74 
Overlap of atomic functions, 148, 151, 
220, 225, 232 
of energy bands, 146, 150 
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Paramagnetism, 33, 34 
of free electrons, 35—36 
Pauli, W. Z., 35, 39 
principle, 10, 19 
Peierls, R. E., 106, 120 
Periodic boundary conditions, see boun- 
dary conditions 
Permutation groups, 72-73 
Perturbation theory, 140, 156 
Phase velocity, 84 
Phonons, 105 
Pines, D., 14, 17 
Pippard, A. B., 172, 184 
Point group, 57 
Primitive cell, 42-43, 45, 46 
Projection operator, 74, 75 


Quadratically integrable, 2 
Quantum, mechanics, 1-4 
numbers, 21, 22 


Reciprocal lattice, 122-124, 129 
Reduced wave vector, 82, 127, 130 
pat ae of a representation, 61~ 
3 
Reflection, 43 
Relaxation time, 192, 193, 195 
Representation, character of, 63 
dimensionality of, 60 
equivalent, 60 
faithful, 63 
irreducible, 61, 135 
matrix, 59-65 
reducible, 63 
Resistivity, 196 
Rhombic Bravais net, 42 
Rotations, 43 


Saddle point, 133 
Schiff, L. I., 1, 17 
Schmidt, L. A., 226, 234 
Screw axis, 57, 137 
Seitz, F., 17, 184 
operator, 49, 51, 134 
Self-consistent field, 13 
Semi-conductor, 121, 162, 185, 194, 206 
n-type, 176 
p-type, 176 
Simple cubic lattice, 148, 149 
Site orbital, 226 
Slater, J. C., 25, 184 
Small groups, 135 
Sommerfeld, A., 18, 40 
theory of metals, 18-23 
Sound, velocity of, 78, 111, 112 
Space groups, 57 
Specific heat, electronic, 31~33 
lattice, 108-115 
Spin paramagnetism, 35-36 
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Splitting of levels, 66-67, 150 
Square, Bravais net, 43 
group, 55-56, 135 
Standing waves, 85, 101 
Star of the wave vector, 135 
Stark ladder, 164, 18r 
Subgroups, 57, 134 
Susceptibility, magnetic, 34, 36 
ayy: of Hamiltonian, 52-53, 65- 
7 


Symmorphic space groups, 134 
Synge, J. L., 94, 120 


Tetrahedral symmetry, 73, 224 
Thermal conductivity, 197-199 
Thermo-electric, effect, 196-201 
power, 200, 201, 204 
Thomson coefficient, 199, 200 
Tight-binding model, 146-151 
Tinkham, M., 77 
Total representation, 97 
Totally symmetric mode, 98 
Translation, group, 100, 128-130 
operator, 49 
symmetry, 41, 42 
Transport equation, 189-193 
Transverse modes, 88 
Tunnelling, 164, 166 


Unit cell, 45, 46 
Unitary operator, 73, 76 
U-process, 187, 188 


Valence, band, 160, 161 

crystals, 217-226 

electrons, 5 
Van der Waal’s force, 211, 213 
Van Hove, L., 106, 120, 132 
Vector, potential, 15, 182 
Velocity of electrons, 155 
Vibrations in solids, 78—116 
Virial theorem, 216-217, 231-232 


Wannier, function, 151-154, 173, 181 
G. H., 151, 184 

Wave vector, 79, 106 

Weak-binding approximation, see 

nearly free-electron model 

Wiedemann—Franz ratio, 199 

Wigner, E., 155, 184 

Wigner-Seitz, cell, 131, 227, 230 
method, 155 
sphere, 228 

Wronskian, 125 


Zener effect, 164-166 
Zero-point energy, 102 
Zone, Brillouin, see Brillouin zone 
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